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Abstract 
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We explore several aspects related to the propagation of high energy cosmic rays (CRs) of 
galactic origin (i.e. in the range E ~ 10^^ — 10^^ eV). In particular, we study in detail 
the diffusion/drift scenario, a rigidity-dependent model in which the knee of the spectrum 
is explained as due to the escape of the light CR component. We show that this scenario 
also explains the second knee (which is produced by the leakage of heavy nuclei), as well as 
the most recent observations of composition and anisotropies (amplitude and first harmonic 
in right ascension) in the range from the knee to the ankle. Moreover, the behavior of 
the diffusion tensor for charged particles propagating in highly turbulent magnetic fields is 
determined numerically. 

We also investigate the high energy neutrino fluxes induced by CRs, and we show that 
a rigidity-dependent scenario for the CR knee has the effect of reducing the background of 
atmospheric neutrinos, hence encouraging the search for high energy astrophysical neutrino 
sources. After exploring this effect in the detection channel of Icpton tracks, wc investigate 
the less studied "shower" channel. As we show here, the shower detection channel is an effec- 
tive means for isolating the prompt atmospheric neutrino component, and it is an important 
complement to the usually-considered track channel. 

Finally, we consider in detail the nuclear photodisintegration model as an alternative 
explanation for the knee in the CR spectrum. In this scenario, CR nuclei emitted by a 
source are disintegrated during their propagation through a surrounding optical/UV photon 
field. Hence, the composition is predicted to turn lighter across the knee, in constrast to the 
results obtained in rigidity-dependent scenarios. 
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Prefacio 



Aiin cuando el fenomeno de los rayos cosmicos viene siendo estudiado en la fisica desde 
hace casi un siglo, sigue siendo en la actualidad un campo de intensa investigacion teorica 
y experimental. Los rayos cosmicos de baja energia, que pucdcn scr dctcctados mcdiantc 
globos, cohetes y satelites, son esencialmente el linico tipo de materia de origen cosmico del 
que puede tenerse evidencia directa, mientras que, en el otro extremo del espectro, los rayos 
cosmicos ultra-energeticos pueden alcanzar energfas de decenas de Joules y constituyen, en 
efecto, las particulas mas energeticas conocidas en el Universo. 

Mas alia del interes intrmseco asociado a cstc fenomeno, la fisica de rayos cosmicos esta 
en intima relacion con otros campos de investigacion (la astronomia, la fisica de particulas, 
la astrofisica, la cosmologia, la fisica nuclear, la fisica de plasmas, etc.), de modo que el 
progreso en esta disciplina tambien contribuye significativamente al desarroUo de otras areas 
de la fisica y la astronomia. 

Pcsc a la notable antigiiedad de este problema, su explicacion no esta aiin bien establecida. 
En verdad, sabemos que el fenomeno de los rayos cosmicos depende del rango de energias 
que sc considcrc, y requiere desarrollar, en cada caso, la explicacion apropiada que de cucnta 
de su naturalcza, de los mecanismos que los originan, y de la manera en que se propagan 
desde las fuentes hasta nuestra ubicacion como observadores. El Capitulo 1 presenta una 
breve introduccion general a la fisica de rayos cosmicos; alii se brinda una revision historica 
de la disciplina, y luego se presentan los aspectos esenciales acerca de la aceleracion, la 
propagacion y la deteccion de los rayos cosmicos. 

Los Capitulos siguientes contienen las contribuciones originales de esta Tesis Doctoral, 
en la que se exploran diversos aspectos relacionados con la propagacion de rayos cosmicos 
galacticos de muy alta energia (es decir, en el rango E ~ 10^^ — 10^* eV, que comprende la 
rodilla y la segunda rodilla del espectro) utilizando diferentes herramientas teoricas, tanto 
analiticas como numericas. 

En el Capitulo 2 se estudia la propagacion de rayos cosmicos en los campos magneticos 
turbulentos de la galaxia; alii se muestra que la rodilla y la segunda rodilla del espectro 
pueden ser explicadas en este escenario como un crossover entre dos rcgfmenes difusivos 
diferentes, a saber, la difusion transversal (que domina a bajas energias) y la antisimetrica 
(o drift) que prevalece a energias altas. Mientras que la supresion de los protones produce la 
primera rodilla, la de la componente pesada (fundamentalmente, niicleos de hierro) da lugar 
a la segunda. La apropiada inclusion de los efectos de drift en las ecuaciones de transporte 
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difusivo no solo explica en forma natural la rodilla y la segunda rodilla del espectro, sino 
tambien el comportamiento observado para la composicion y las anisotropi'as en el rango 
de energias que se extiende desde la rodilla hasta el tobillo del espectro. Por otro lado, 
86 Ueva a cabo un estudio del comportamiento del tensor de difusion en condiciones de 
alta turbulencia. Una aplicacion particular de este estudio es la descripcion detallada de la 
difusion de rayos cosmicos en la galaxia para energias en la region de la rodilla, aunque los 
resultados satisfacen relaciones de escala que permiten ser aplicados en otros problemas de 
interes astrofisico. Los resultados que se presentan en este Capitulo fueron publicados en los 
siguientes articulos: 

• J. Candia, E. Roulet and L.N. Epele, Turbulent diffusion and drift in galactic magnetic 
fields and the explanation of the knee in the cosmic ray spectrum, Journal of High 
Energy Physics JHEP 12 (2002) 033. 

• J. Candia, S. Mollerach and E. Roulet, Cosmic ray drift, the second knee and galactic 
anisotropics, Journal of High Energy Physics JHEP 12 (2002) 032. 

• J. Candia, S. Mollerach and E. Roulet, Cosmic ray spectrum and anisotropics from 
the knee to the second knee. Journal of Cosmology and Astroparticle Physics JCAP 05 
(2003) 003. 

• J. Candia and E. Roulet, Diffusion and drift of cosmic rays in highly turbulent magnetic 
fields. Journal of Cosmology and Astroparticle Physics JCAP 10 (2004) 007. 

El Capitulo 3 investiga el impacto en la prediccion de los fiujos de neutrinos de altas en- 
ergias que rcsulta dc considcrar escenarios en los que la rodilla del espectro de rayos cosmicos 
es un efecto dependiente de la rigidez magnetica, como corresponde, por ejemplo, al esce- 
nario de la difusion turbulenta y drift. Se encuentra que el fondo de neutrinos atmosfericos 
se reduce significativamente, facilitando la biisqueda de fuentes astroffsicas. Asi, se pone en 
evidencia la necesidad de resolver el origen de la rodilla, con el fin de interpretar correcta- 
mente las observaciones en los nuevos telescopios de neutrinos. Este fenomeno se estudia en 
el canal usual de trazas leptonicas (asociado a la deteccion de neutrinos muon), y luego en el 
canal complementario de las cascadas hadronica y electromagnetica (asociado a la deteccion 
de neutrinos electron), que se propone como un nuevo metodo para aislar el flujo de neutrinos 
atmosfericos prompt. Estas investigaciones dieron lugar a las siguientes publicaciones: 

• J. Candia and E. Roulet, Rigidity dependent knee and cosmic ray induced high energy 
neutrino fluxes. Journal of Cosmology and Astroparticle Physics JCAP 09 (2003) 005. 

• J.F. Beacom and J. Candia, Shower power: isolating the prompt atmospheric neutri- 
no flux using electron neutrinos. Journal of Cosmology and Astroparticle Physics (en 
prensa) . 
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El Capitulo 4 explora en detalle el modelo en el que la rodilla del espectro se debe a 
efectos de fotodesintegracion nuclear por la interaccion con fotones opticos y UV blandos 
en la region de las fuentes. Este escenario ofrece una explicacion alternativa a la rodilla 
del espectro, y predice, contrariamente a lo que resulta de los escenarios dependientes de la 
rigidez, una composicion que se torna mas liviana por encima de la rodilla. Este trabajo fue 
publicado en el articulo siguiente: 

• J. Candia, L.N. Epele and E. Roulet, Cosmic ray photodisintegration and the knee of 
the spectrum, Astroparticle Physics 17 (2002) 23. 

El Capitulo 5 contiene las Conclusiones de esta Tesis. Finalmente, se presentan los 
Apendices y la Bibliografia. 

Esta Tesis Doctoral se desarroUo con el apoyo de una beca de postgrado interna otorgada 
por el Consejo Nacional dc Investigaciones Cientificas y Tecnicas (CONICET) en el pen'odo 
04/2000 - 03/2004, y un subsidio otorgado por la Fundacion Antorchas en cl periodo 05/2003 
- 04/2004. Este trabajo fue concluido en Fermilab (Batavia IL, USA) en el marco del 
Programa para doctorandos latinoamericanos, en el periodo 03/2004 - 08/2004. 
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Capitulo 1 

La fisica de rayos cosmicos 



Aqm se ofrece un breve resumen de los aspectos esenciales de la fisica de rayos 
cosmicos, con cl objctivo clc prcscntar un contcxto a las contribuciones originalcs 
dc csta Tcsis. Sc incluyc una lista dc publicaciones de referenda que, desde 
diferentes enfoques, describen los progresos en este campo de investigacion. 

1.1 Introduccion 

Existe una extensa bibliografia que cubre el desarrollo historico de la fisica de rayos cosmicos, 
que expone con gran detalle diversos aspectos de esta disciplina, a saber, el origen, la propa- 
gacion y la deteccion de los rayos cosmicos en los diferentes rangos de energia relevantes, 
y que resume las principales contribuciones que se han desarroUado, desde la teoria o el 
experimento, para comprender este fenomeno. 

Los volumenes de Longair [1] presentan la astroffsica de altas energfas a un nivel de 
pregrado, poniendo especial entasis en el estudio de los rayos cosmicos; Gaisser [2] desar- 
roUa la fisica de rayos cosmicos desde su conexion con la fisica de particulas, y dedica una 
considerable atencion a diversos aspectos de la deteccion (desarrollo de Uuvias atmosfericas, 
parametrizacion de los observables, detectores, tecnicas de simulacion, etc.); Berezinskii et 
al. [3] presentan un estudio abarcativo de toda la astroffsica de rayos cosmicos, que resume 
el estado de la disciplina hasta 1990; mientras que Sokolsky [4] se ocupa, en particular, de los 
rayos cosmicos ultra-energeticos. Un libro mas reciente, de Schlickeiser [5], trata con gran 
detalle los aspectos asociados a la propagacion. 

Con respecto a los trabajos de revision, aqui citaremos a Blandford y Eichler [6] en 
mecanismos de aceleracion, a Cesarsky [7,8] en propagacion, y a Wiebel-Sooth y Biermann [9] 
en la compilacion dc datos experimentalcs. Trabajos de revision mas breves y recientes que 
cubren diversos aspectos de la fisica de rayos cosmicos son, por ejemplo, los de Boratav y 
Watson [10], Kampert [11], Biermann y Sigl [12], Roulet [13] y Horandel [14]; existe ademas 
un gran niimero de publicaciones de revision referidas especificamente a los rayos cosmicos 
ultra-energeticos (por ejemplo, [15-20]), un topico de gran interes actual, pero que escapa al 
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tema especifico de esta Tesis. Finalmente, deben mencionarse los Proceedings of the ICRC 
(International Cosmic Ray Conference)^ las reuniones internacionales bienales que convocan 
a una parte significativa de la comunidad cientifica del area. Las reuniones mas recientes a 
la fecha fueron las de Utah (1999), Hamburg (2001) y Tsukuba (2003). 

1.2 Una breve revision historica 

El cstudio dc los rayos cosmicos se origino a comienzos del siglo XX, vinculado a las in- 
vestigacioncs sobrc radioactividad que se llevaban a cabo con gran exito en esa epoca. Un 
instrumcnto muy utilizado en csos cxperimentos era el elcctroscopio, que pcrmitia medir la 
ionizacion del gas contenido en el aparato (y, en consecuencia, la cantidad de radiacion recibi- 
da) a partir de la velocidad con que las laminas cargadas del instrumento se aproximaban 
entre si. Una peculiar observacion fue que las laminas del electroscopio resultaban descar- 
garse aun en ausencia de fuentes radioactivas; con el fin de explicarla, se realizaron diversos 
cxperimentos que, en principio, buscaban atribuir el fenomeno a la radioactividad natural 
de la superficie terrestre, a la del recipiente, o bien a posibles emanaciones radioactivas en 
el propio gas. 

En 1910, Wulf encontro que la radiacion descendia un 40% al subir a la torre Eiffel, de 
330 m de altura. Estc experimento sugeria, dc todos modos, que la radiacion observada 
no era emitida por la superficie de la tierra; usando el coeficiente de absorcion de los rayos 
7 en el aire, ya conocido en ese momento, se esperaba que la radiacion descendiera a la 
mitad en solo 80 m, de forma que debia ser totalmente despreciable en el punto mas alto 
de la torre. La evidencia definitiva sobre el origen cosmico de la radiacion provino de los 
vuelos en globo realizados por Hess a partir dc 1912. En particular, en su vuelo del 7 de 
agosto de 1912, alcanzo una altitud de 5 km y pudo observar que la tasa de ionizacion 
resultaba varias veces mayor que la medida a nivel del mar ^; estos resultados luego fueron 
confirmados por Kolhorster en vuelos que realize entre 1913 y 1914, alcanzando altitudes de 
hasta 9 km. Durante los anos siguientes, se estudio este fenomeno asumiendo que se trataba 
de algun tipo de radiacion 7 con un poder de penetracion mayor que el obscrvado en la 
radioactividad natural, lo que incidentalmente sugirio a Millikan, en 1925, llamarla con el 
nombre que perdura hasta nuestros dias. En 1927 fue descubierto el efecto geomagnetico (es 
decir, que la ionizacion producida por los rayos cosmicos depende de la latitud), lo que dio 
lugar a establecer definitivamente hacia 1936 que los rayos cosmicos primaries eran particular 
cargadas; mas tarde se demostro que los rayos cosmicos eran fundamentalmente protones, 
aunque tambien contenian niicleos de elementos mas pesados, y una proporcion pequena 
(inferior al 1%) de electrones relativistas. 

^La conclusion de las observaciones de Hess dio lugar al descubrimiento definitivo de los rayos cosmicos: 
"The results of the present observations seem to be most readily explained by the assumption that a radiation 
of very high penetrating power enters our atmosphere from above, and still produces in the lower layers a 
part of the ionisation observed in closed vessels." [21]. 
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El desarroUo de nuevos instrumentos de medicion, como por ejemplo la camara de niebla 
de Skobeltsyn, o el detector Geiger-Miiller, ambos inventados en 1929, hizo posible profun- 
dizar las investigaciones sobre la naturaleza de los rayos cosmicos, y promovio a la vez una 
gran evolucion en la fisica de particulas. En efecto, durante las dos decadas siguientes, los 
rayos cosmicos representaron un acelerador de particulas natural, es decir, una fuente gratui- 
ta e inagotable de particulas de alta energia con la atmosfera como bianco; solo era necesario 
construir el detector para estudiar los productos resultantes de las colisiones. Siguiendo 
esa Imea de investigacion, en ese periodo fueron descubiertos los positrones, muones, piones, 
kaones e hyperones; mas tarde, a partir de los anos '50, la tecnologia de los aceleradores Uego 
a producir haces de particulas de energias comparables a la contenida en los rayos cosmicos, 
con la ventaja adicional de poder manipular y controlar los eventos producidos, y entonces 
la fisica de particulas se traslado definitivamente a los grandes laboratorios. 

Entre las particularidades mas notables del fenomeno de los rayos cosmicos, sin duda la 
mas llamativa es el espectro (ver la fig Jl-H) . que se presenta como una ley de potencias que 
se extiende por mas de 11 decadas en la energia, desde 10^ hasta mas alia de 10^*^ eV, con un 
indice espectral esencialmente constante. El espectro dN/dE oc tiene un indice espectral 
a ~ 2.7 a bajas energias, luego un cambio liacia a ~ 3 que ocurre en i^^ — 3 x 10^^ eV, un 
fenomeno conocido como la rodilla del espectro, un segundo cambio con a ~ 3.3, que aparece 
para Esr — 4 x 10^^ eV, la segunda rodilla, y finalmente un endurecimiento del espectro liacia 
a ~ 2.7 en ^ 5 X 10^^ eV, el tobillo. 

Para energias por debajo de 10^^ eV, los rayos cosmicos primarios pueden ser observados 
en forma directa mediante detectores a bordo de globos, colietes y satelites. Sin embargo, 
a energias mayores esto no es posible en general, debido a que el fiujo cae rapidamente con 
la energia; por ejemplo, el fiujo integral en la region de la rodilla es menor a 1 particula por 
m^ por alio. Para su deteccion, en cambio, debe recurrirse a la reconstruccion de las lluvias 
atmosfericas {extensive air showers, EAS), es decir, observando la cascada de particulas 
secundarias producidas por la incidencia de los rayos cosmicos sobre la atmosfera de la tierra. 
Las lluvias atmosfericas comenzaron a ser exploradas en los afios '30, a partir de los trabajos 
de Rossi ^ y de Auger y sus colaboradores [23,24], midiendo la coincidencia temporal en 
arreglos de contadores Geiger-Miiller. Las tecnicas experimentales fueron progresando con 
el tiempo, permitiendo acceder a la reconstruccion de eventos de mayores energias; el grupo 
del MIT [25] incorporo los centelladores como detectores en un arreglo que, operando entre 
1954 y 1957, permitio obtener el espectro en el rango desde 3 x 10^^ hasta 10^^ eV [26]. 
Este arreglo, llamado Harvard Agassiz Station, sirvio ademas como prototipo para el arreglo 
de gran altitud en Chacaltaya (Bolivia) y el de Volcano Ranch en New Mexico (Estados 
Unidos). En este ultimo, Linsley detecto el primer evento por encima de 10^*^ eV [27]; la 
deteccion de lluvias atmosfericas mediante arreglos de superficie extendidos fue, mas tarde, 

^Rossi parece haber sido el primero en sospechar la existencia de las lluvias atmosfericas, en base a 
sus observaciones en Eritrea. En 1934, escribio: "It would seem ... from time to time there arrive upon 
the equipment very extensive groups of particles which produce coincidences between counters even rather 
distant from each other" [22]. 
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Figura 1.1: El espectro de rayos cosmicos en funcion de la energia. El mdice espectral es 
aproximadamente constante; sin embargo, exhibe cambios, pequeiios pero bien definidos, 
que dan lugar a la rodilla, la segunda rodilla y el tobillo del espectro, senalados en la figura 
(de [16]). 



exitosamente utilizada en otros sitios: Haverah Park (Gran Bretana), Narribri (Australia), 
Yakutsk (Rusia) y AGASA (Japon). 

Tambien fueron explorados metodos alternativos de deteccion; en particular, Suga [28] 
y Ghudakov [29] discutieron en 1962 la posibilidad de utifizar la atmosfera terrestre como 
un centellador gigante, una idea que fue luego investigada en el terreno experimental por 
Greisen en 1965 [30], y por Tanahashi y colaboradores en 1968 [31]. En 1976, se detecto 
en Volcano Ranch la emision fluorescentc de las lluvias en coincidencia con el arreglo de 
detectores de superficie [32], dando lugar al desarroUo de experimentos basados solo en esa 
tecnica, a saber, los exitosos Fly's Eye y HiRes situados en Utah (Estados Unidos). 

Una nueva generacion de experimentos esta por comenzar con el Observatorio Auger en 
Malargiie (Argentina) [33], que consiste en un arreglo hibrido formado por detectores de 
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superficie y de fluorescencia; la operacion combinada de las dos tecnicas es potencialmente 
muy poderosa, y permitiria eliminar las discrepancias sistematicas que se observan al com- 
parar las obscrvacioncs de experimentos que utilizan solo un metodo de deteccion, como es 
el caso, por ejemplo, de AGASA y HiRes. For otra parte, la superficie cubierta por los de- 
tectores sera de 3000 km^, considerablemente mayor que los 100 km^ cubiertos por AGASA. 
La publicacion de los primeros resultados experimentales se espera para mcdiados de 2005, y 
correspondera a una estadistica comparable a la acumulada por AGASA desde su apertura 
en 1990 hasta la fecha. 

Otros experimentos han sido proyectados y se encuentran en diversas fases de planeamien- 
to. Como parte del Proyecto Auger, se planea construir otro Observatorio en Utah (Estados 
Unidos), con el proposito de cubrir toda la esfera celeste. Otro proyecto, basado en la tccnica 
dc la fluorescencia, es cl Telescope Array Project [34], que permitiria la reconstruccion es- 
tereoscopica de todos los eventos ultra-energeticos desde 10^^ eV, con un cxcelente potencial 
para discriminar la especie del rayo cosmico primario. Finalmente, las propuestas mas am- 
biciosas consideran la observacion de la luz fluorescente desde satelites en el espacio, dando 
lugar al proyecto Airwatch/OWL [35]. La observacion desde el espacio ofreceria, por un 
lado, la ventaja dc cubrir un area mucho mayor que desde la tierra, mientras que, adcmas, 
la absorcion de los fotones emitidos seria mucho menor, incrementando la sensibilidad y 
reduciendo los errores sistematicos. La reahzacion de cualquiera de estos grandes proyectos 
quedara, muy probablemente, supeditada a los resultados que se obtengan en los proximos 
aiios en el Observatorio Auger Sur. 

1.3 La aceleracion de rayos cosmicos 

Explicar el origen de los rayos cosmicos rcquiere, por un lado, especiflcar un mecanismo de 
aceleracion que produzca protones y niicleos con espectros de ley de potencias, dN/dE oc 
E~^, con los indices espectrales adecuados a las observaciones. Como veremos mas adelante, 
en el Capitulo 2, la difusion normal domina el transporte de rayos cosmicos a bajas energias 
(debajo de la rodilla), y de ello se espera que (3 = a — 1/3, donde a es el indice del cspcctro 
observado. Siendo que, tipicamente, se obscrva a ~ 2.7 para los protones y las espccics 
nucleares mas abundantes, deberia esperarse {3 ~ 2.4 para el indice de produccion en las 
fuentes ^. Por otra parte, se necesita un mecanismo que permita acelerar rayos cosmicos 
hasta energias extremadamente altas; aun dejando de lado los problemas especificos de los 
rayos cosmicos ultra-energeticos, se requiere un mecanismo que acelere particulas, al menos, 
hasta la region del tobillo del espectro. Y flnalmcnte, asumiendo que la cfcctividad del 
mecanismo sea del ordcn del ~ %, se necesita, adcmas, que la luminosidad total de las 

•^Una observacion notable es que el mecanismo de aceleracion de particulas hasta energias ultrarrelativistas 
parece ser similar para una variedad de diferentes fuentes astrofisicas; por ejemplo, los espectros tfpicos de 
fuentes dc radio corresponden a espectros de elcctroncs con indice /3 « 2.6 ± 0.4, y los espectros continuos 
de cuasares en longitudes de onda opticas y X corresponden a /? ~ 3. 
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fuentes sea mayor, en alrededor de 2 ordenes de magnitud, que la correspondiente a los rayos 
cosmicos. 

Aqui veremos que un mecanismo de aceleracion apropiado es el que resulta de procesos 
difusivos en frentes de onda de choque; en particular, la aceleracion en remanentes de su- 
pernova {supernova remnants, SNRs) satisface los requerimientos que hemos mencionado, 
y constituye un candidate preferencial en la explicacion del origen de los rayos cosmicos. 
La propuesta original de estc mecanismo fuc desarrollada por Fermi cn 1949 [36], y en su 
formulacion actual se conoce con el nombre de "mecanismo de Fermi de segundo orden" . La 
idea basica parte de considerar la transferencia de energia cinetica adquirida por particulas 
cargadas que colisionan con objetos masivos en movimiento; Uamando ^ = /S.E/E a la en- 
ergia relativa ganada por una partfcula en cada colision, al cabo de n colisiones la energia 
de la partfcula es 

E = E,{l+iY, (1.1) 

donde Eq es la energia inicial. Si Tf. es el tiempo caractcristico cntre colisiones sucesivas, y 
Tesc el tiempo caractcristico dc escape de la region de aceleracion, la probabilidad de escape 
luego de una colision es Pesc = Tc/Tesc) entonces, la densidad de particulas aceleradas hasta 
una energia mayor o igual a E resulta 

N{> E)o,Y,{l- Pescr = ^— ^ , (1.2) 



esc 



donde n — log(£'/£'o)/log(l -|- ^). Entonces, la densidad queda dada por 

N{E) oc ^-(W^+i) ; (1.3) 

vcmos que este proceso conduce a un espectro de ley de potencias hasta una energia maxima 
Emax dada por 

Emax ^ Eo{l + ^)-^MTM/T. ^ (14^ 

donde T4 es la vida del acelerador. 

En la conccpcion original dc Fermi, la rcalizacion astrofisica de este mecanismo era la 
reflexion de las particulas cargadas en "espejos magneticos" asociados con irregularidadcs 
en el campo magnetico de la galaxia (es decir, nubes magnetohidrodinamicas no relativistas 
propagandose en el medio intcrcstelar). Como la nube tiene una masa mucho mayor que la 
partfcula, su velocidad no cambia con la colision; en el sistema inercial asociado a la nube, el 
choque es elastico y la partfcula se rcflcja conservando su energfa. Transformando al sistema 
del observador, en el que la velocidad de la nube es V, la de la partfcula es w ~ c, y 6* es el 
angulo formado entre la velocidad de la nube y la velocidad de salida de la partfcula, resulta 

AE 2yvcos^ 



E c2 V c 



+ 2 - . (1.5) 
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Figura 1.2: La aceleracion de particulas en la expansion de frentes de onda de choque (de [1]). 
(a) La propagacion del frente de ondas entre el upstream (a la dereclia) y el downstream 
(a la izquierda). (b) El fiujo de gas interestelar a traves del frente de ondas, en el sistema 
en el que el frente esta en reposo. (c) El fiujo de gas en el sistema del upstream, con el 
frente de onda incidiendo sobre una distribucion de particulas isotropica. (d) El fiujo de gas 
en el sistema del downstream, con el frente de onda incidiendo sobre una distribucion de 
particulas isotropica. 



La ganancia media de energfa se obtiene ahora promediando sobre el angulo 9; debe tenerse 
en cuenta que, de la cinematica de la colision, la probabilidad de un encuentro con un dado 
angulo 9 viene dada por P{9) oc 1 + (V^/c)cos6' [1], de modo que los clioques frontales 
(que corresponden aO < 9 < 7r/2) son mas probables que los clioques desde atras (que 
corresponden a 71/2 < 9 < tt). El resultado es 



E ' 3 _ 

de modo que la ganancia media por colision es de segundo orden en V/c. Siendo que las 
velocidades estocasticas de las nubes interestelares en la galaxia son tfpicamente V/c < 10~'^, 
este resultado dari'a una ganancia de energfa insuficiente; por otra parte, asumiendo ~ 
1/cpa (es decir, considerando las nubes como centros dispersores de densidad p y seccion 
eficaz a), el I'ndice espectral /3 resulta depender de propiedades locales, en lugar de ser 
universal. 

Hacia fines de los anos '70, se encontro [37-40] que estos problemas podian corregirse 
considerando la colision con frentes de onda de choque relativistas, dando lugar a un mecan- 
ismo de aceleracion universal con una ganancia media por colision de primer orden en V/c. 
Como sabemos, una perturbacion en un gas produce ondas que se propagan a la velocidad 
del sonido y transmiten informacion sobre la fuente; si esta se mueve a una velocidad mayor 
a la del sonido, la perturbacion del medio se propaga en la forma de una onda de choque (con 
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la velocidad de la fuente), que representa una discontinuidad entre el estado termodinamico 
del medio perturbado {downstream) y el del medio no perturbado {upstream). Denotare- 
mos la presion, temperatura y densidad del upstream con pi,Ti j pi, respectivamente; a las 
variables termodinamicas del downstream las Uamaremos P2,^2 y P2- La fig ll.2f a) muestra 
el avance del frente de onda de choque con velocidad U = Mcg, donde Cg es la velocidad 
del sonido en el gas y M ^ 1 el numero de Mach. En el sistema de referenda del frente 
de onda, el gas interestelar incide desde el upstream con velocidad fi = |f/| y se aleja del 
frente en el downstream con velocidad V2- En ese sistema de referenda, se pueden rela- 
cionar las propiedades termodinamicas del downstream relativas al upstream, planteando la 
conservacion de la masa (es decir, la ecuacion de continuidad) , de la energfa (la ecuacion 
de Bernoulli) y del impulso, para el flujo de gas desde el upstream hacia el downstream. 
Si consideramos un gas ideal, siendo 7 = Cp/cy el cociente de calores especifi'cos, resulta 
P2/P1 = Vi/v2 = (7 + l)/(7 — 1); en particular, para un gas completamente ionizado, 7 = | 
y resulta f2 = fi/4 (ver fig |1.2f b)). 

Consideremos ahora una distribucion de parti'culas de alta energfa delante del frente 
de ondas. En el sistema en el que el gas esta en reposo, la distribucion de parti'culas es 
isotropica, debido al scattering. Entonces, visto desde el upstream (ver la figlOfc)), el gas 
en el downstream se aproxima con velocidad V = |[/ y la distribucion de velocidades de 
las parti'culas que se encuentran con el frente de ondas es isotropica. Si consideramos ahora 
el proceso opuesto, en el que las parti'culas pasan a traves del frente desde el downstream 
hacia el upstream, y lo observamos desde el sistema del downstream (ver la fig ll.2r d)). se 
encuentra lo mismo: la distribucion de particulas es isotropica y el upstream se aproxima con 
velocidad V = jU. Es decir, una particula que pasa desde el downstream hacia el upstream 
experimenta un incremento AE en su energfa, que resulta igual al del proceso inverso. En 
el mecanismo original de Fermi, el cambio medio de energia en una colision era de orden 
V/c, pero debi'an promediarse colisiones frontales (que producian ganancias de energfa) con 
colisiones desde atras (que producfan perdidas); del promedio, resultaba una ganancia neta 
de orden {V/c)'^. En cambio, aquf el mecanismo asegura una completa simetri'a en el pasaje 
de partfculas de un lado a otro del frente de ondas, ganando energfa en cada colision. Al ir 
y volver, la energfa ganada resulta [1] 



que es de primer orden en V/c; por otra parte, el I'ndice espectral resulta (3 = 2 + 
Este resultado es muy satisfactorio, porque reproduce los valores esperados y, ademas, per- 
mite explicar la universalidad observada en este fenomeno; la presencia de frentes de onda 
supersonicos es plausible en cualquier fuente astrofi'sica de parti'culas de alta energfa, tales 
como remanentes de supernova, micleos galacticos activos y fuentes de radio extendidas. 

Como hemos comentado mas arriba, un escenario ampliamente difundido y aceptado 
en la actualidad es el que atribuye la aceleracion de los rayos cosmicos (exceptuando a los 
ultra-energeticos) a los remanentes de supernova; dependiendo de si los frentes de onda se 
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1 au 1 pc 1 kpc 1 Mpc 

log (size, km) 

Figura 1.3: El plot de Hillas: el campo magnetico y el tamano caracteristicos de diversos 
objetos astrofi'sicos, comparados con la energia maxima que puede adquirir una particula en 
los sitios de aceleracion (de [19]). 

desarroUan en el medio interestelar, o bien en el viento estelar de una estrella predecesora o en 
el de una compaiiera [41], los protones podrian ser acelerados hasta una energia maxima en el 
rango entre Ec = 10^^ y 10^^ eV, mientras que, para micleos de carga Ze, esta energia maxima 
vendria dada por Ecz — ZEc Debe indicarse, sin embargo, que aiin no ha sido posible 
obtener evidencia directa de la aceleracion de protones en remanentes de supernova [42] . En 
el marco del paradigma de las supernovas, el tobillo en ~ 5 x 10^^ eV aparece naturalmente 
como un crossover entre una region dominada por los rayos cosmicos galacticos (para E < Et) 
y otra dominada por los rayos cosmicos extragalacticos (para E > Et). 

La energia maxima alcanzable en un proceso difusivo de aceleracion en frentes de onda 
de choque esta dada por E^ax = kZeUBL [43,44], donde k esta asociado a la eficiencia del 
mecanismo de aceleracion (A; < 1), i? es el campo magnetico en la region de aceleracion, y L 
el tamafio del acelerador. En el caso de aceleracion optima, que corresponde a tomar k = 1 

"^Otros mecanismos permitirian acelerar rayos cosmicos hasta energias aiin mayores, como es el caso, por 
ejemplo, de la aceleracion electrostatica en piilsares jovenes. Sin embargo, estos mecanismos de aceleracion 
alternativos tienen por lo general dificultades para predecir espectros de ley de potencias, y usualmente son 
propuestos solo para explicar el origen de los rayos cosmicos ultra-energeticos. 
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J U = c, se obtiene ^ 



'max 




kpc 



L 



) 



eV. 



(1.8) 



En particular, este limite es muy ilustrativo cuando se lo grafica en un diagrama B vs L 
comparado con los valores tipicos de diferentes objetos astrofisicos, tal como se muestra en la 
fig J1.3l Las rectas corresponden a Emax = 10^° eV para protones y niicleos de Fe, y diferentes 
valores de /? = U/c; solo las fuentes ubicadas en la region superior satisfacen la condicion 
impuesta por la ec. ()1.8|) . Este grafico, conocido como "plot de Hillas", pone en evidencia 
las dificultades que aparecen para encontrar objetos astrofisicos que provean un mecanismo 
apropiado para acelerar rayos cosmicos hasta las mas altas energias observadas. 

1.4 La propagacion de rayos cosmicos 

Como veremos, la abundancia qufmica de los rayos cosmicos de bajas energias, que puede 
obtenerse mediante tecnicas de deteccion directa, provee importante informacion acerca de 
los potenciales sitios de aceleracion y, fundamentalmente, sobre diferentes aspectos referidos 
a la propagacion de los rayos cosmicos en la galaxia. 

La espectroscopia astronomica busca determinar las abundancias qufmicas de los ele- 
mentos en diferentes objetos celestes; en particular, las lineas de absorcion debido a iones, 
atomos y moleculas en la fotosfera solar permiten medir las abundancias en el sol. De la 
teorfa de la evolucion estelar, sabemos que la abundancia cosmica tipica se adquiere mas bien 
rapidamente, de modo que la abundancia solar es representativa de la de toda la galaxia ^. 
Otro medio para determinar las abundancias es el estudio de los meteoritos; si bien se en- 
cuentran diferencias muy significativas en las abundancias en diferentes tipos de meteor ito, 
que probablemente reflejan sus diversas procedencias, existe una clase de meteorito partic- 
ular, los condriticos, que tienen abundancias muy similares a los de la fotosfera solar. Esto 
result a sumamente valioso para extraer informacion sobre las abundancias de algunos de los 
elementos mas raros de la tabla periodica, que no pueden ser bien determinadas mediante 
tecnicas espectroscopicas. 

La fig J1.4l muestra una comparacion entre la abundancia galactica local y la observada en 
los rayos cosmicos; la impresion general es que las abundancias son muy similares, de modo 
que los rayos cosmicos deben ser acelerados a partir de un material de composicion qufmica 
muy similar a la observada en el sistema solar. Sin embargo, se observa en los rayos cosmicos 
un marcado exceso de Li, Be y B, y tambien en el grupo de elementos pesados justo debajo 

^Las unidades astronomicas de longitud usadas habitualmente son el parsec o segundo-parallax (pc), el 
aiio luz (ly) y la unidad astronomica (AU), relacionadas por 1 pc = 3.26 ly — 2.06x10^ AU — 3.09x10^^ m . 

^El sol tiene una edad de 4.6 x 10^ anos y la galaxia, al menos, el doble; las estrellas mas viejas presentan 
un deficit en elementos pesados, pero puede asumirse que la abundancia cosmica no se modifico notablemente 
desde la epoca de la formacion del sol. 
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Figura 1.4: Comparacion entre la abundancia galactica local y la observada en los rayos 
cosmicos (de [45]). 



de Fe. Otras sobreabundancias muy significativas, que no se muestran en la figura, son las 
de los isotopos ^He y ^He; por ejemplo, las abundancias de ^He/^He y de '^He/'^He en los 
rayos cosmicos son unas 5000 veces mayores que las abundancias cosmicas correspondientes. 
Como veremos, estas anomalias pueden explicarse como debidas al proceso de fragmentacion 
nuclear (spallation) de las especies mas pesadas; ademas, esto permite determinar la region 
de la galaxia en la que los rayos cosmicos cstan confinados. 

La ecuacion de transporte para micleos de alta energia propagandose a traves del gas 
interestelar es 

^ = V{DVN, - Vm + ^{b{E)N,) + Q^-^ + ^ + Y.^^^- (1-9) 
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Aqui, Ni = Ni{E,x,t) es la densidad de niicleos de la especie esima (donde el submdice 
ordena las especies de acuerdo con su masa, en forma creciente) , D el tensor de difusion, V 
es una velocidad de conveccion (que puede reprcscntar, por ejemplo, un viento galactico), 
d{b{E)Ni)/dE toma en cuenta ganancias y pcrdidas continuas de energia, donde b{E) = 
—dE/dt (es decir, b{E) < corresponde a una aceleracion sistematica promedio, como la 
debida, por ejemplo, a la aceleracion estocastica de Fermi, mientras que b{E) > representa 
perdidas continuas de energia debidas a la emision de radiacion por bremsstrahlung, a la 
emision sincrotronica, a las perdidas por ionizacion, etc.), Qi representa la tasa de inycccion 
de particulas en las fuentes por unidad de volumen, y rf es la vida media de decaimicnto 
radiactivo para los micleos inestables. Los dos ultimos terminos describen el efecto de ganan- 
cias y perdidas debidas a la fragmentacion nuclear, siendo r/*' la vida media de la especie 
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z— esima asociada a este proceso; a Ni contribuye la fragment acion de todas las especies mas 
pesadas, como indica la suma sobre j > i en el ultimo termino, que acopla las ecuaciones 
correspondientes a diferentes componentes nucleares. Un caso particular de gran interes es 
el correspondiente al estado estacionario, en el que la dependencia explicita con el tiempo 
desaparece. 

Basados en la ecuacion general de transporte dada por la ec. (|1.9|1 . adoptaremos un mod- 
elo muy simplificado (conocido en la literatura como modelo slab [1]) con el fin de ilustrar la 
manera con que el proceso de spallation explica las sobreabundancias observadas. Desprecia- 
remos la difusion, la conveccion, las ganancias y perdidas continuas de energia, la inyeccion 
de nuevas particulas, y el decaimiento radioactivo; solo consideraremos la propagacion rec- 
tilfnea de los micleos, agrupados en el grupo liviano L (formado por Li, Be y B) y el grupo 
medio M (constituido fundamentalmente por C, N y O). Siguiendo la trayectoria de un grupo 
de micleos que se propaga con velocidad v = dx/dt, veremos como es la dependencia espacial 
de la densidad N entre la fuente y el observador, suponiendo que la produccion continua 
en las fuentes da lugar a un estado estacionario. Es conveniente adoptar como variable a 
la columna de densidad ^ = px = pvt, donde p es la densidad promedio en la region de 
propagacion. Entonces, las ecuaciones para la densidad de los grupos L y M se reducen a 



M 



asumiremos, ademas, que iV£,(0) = 0, es decir, que la densidad del grupo liviano producido 
en las fuentes es despreciable. El camino libre medio de spallation viene dado por m = 
nT'p/aL,M] reemplazando las secciones eficaces de spallation, ctl — 200 mb y um — 280 mb, 
resultan — 84 kg m~^ y C,m — 60 kg m~^, respectivamente. Para la probabilidad de 
fragmentacion promedio, puede tomarse el valor Pml — 0.28. Integrando las ecs. p.lO|l y 
()1.11|) . se obtiene 



exp ( - ^ ) - 1 



;i.i2) 



Nm{0 (a -a/) 

De los datos experimentales, sabemos que Nl{^)/Nm{0 = 0.25; reemplazando este valor 
en la ecuacion precedente, resulta ^ = 48 kg m~^. El mismo tipo de calculo puede hacerse 
para la produccion de ^He por spallation de "^He, y se obtiene ^ = 50 kg m~^, es decir, un 
result ado muy similar al obtenido para los micleos del grupo L. 

Cuando se busca explicar la sobreabundacia del grupo pesado (Mn, Cr, V) por spallation 
de los micleos de Fe, aparece una dificultad; la seccion eficaz del Fe es mucho mayor que la de 
elementos mas livianos, y resulta en un camino libre medio mucho menor, = 22 kg m~^, 
de mo do que la abundancia del Fe quedaria excesivamente suprimida luego de atravesar una 
columna de densidad ^ ~ 50 kg m~^. La solucion se obtiene de considerar que los micleos 
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no recorren una linica columna de densidad, como en el modelo slab, sino asumiendo una 
distribucion de longitudes de recorrido, -P(0- 

La ecuacion de difusion y escape es una version simplificada de la ecuacion de transporte 
()1.9|) . en la que se considera que los rayos cosmicos escapan de la galaxia con un tiempo 
caracteristico t^] es decir, 



En el modelo conocido como leaky box, se asume una densidad homogenea N en toda la 
region de propagacion; luego, el termino difusivo en la ec. p.l3|) se anula, y resulta una 
distribucion exponencial de columnas de densidad -P(0 ^ 6xp(— ^/^g)- Alternativamente, 
asumiendo la difusion como el termino preponderante, es decir, tomando Tg — > oo, se obtiene 
una distribucion gaussiana. En cualquier caso, resulta que la ecuacion de transporte conduce 
a distribuciones de caminos de densidad que resuelven el problema que aparecia con la abun- 
dancia del Fe, y que permiten, ademas, reproducir adecuadamente todas las abundancias 
observadas. De deducir la abundancia en las fuentes y compararla con los valores tipicos 
para el sistema solar, se encuentra que los elementos cuyo primer potencial de ionizacion es 
mayor a 10 eV tienen una abundancia suprimida por un factor ~ 5; esto sugiere la existen- 
cia de algiin mecanismo de supresion para la produccion de elementos con potenciales de 
primera ionizacion altos. Por otro lado, en experimentos que resuelven la carga y la masa de 
las particulas se observan otras anomali'as, como por ejemplo la sobreabundancia isotopica 
en ^^Ne/^°Ne, ^^'^^Mg/^^Mg y ^^''^°Si/^^Si, sugiriendo que al menos una parte de los rayos 
cosmicos se origina en medios ricos en neutrones. 

Hasta aqui, hemos visto que la abundancia observada para las especies estables se ex- 
plica a partir de la fragmentacion nuclear; esta explicacion permite, ademas, determinar 
la abundancia en las fuentes y la columna de densidad recorrida por los rayos cosmicos en 
su propagacion. Sin embargo, para determinar el volumen de confinamiento y el tiempo 
de residencia (o de escape), es necesario recurrir a la abundancia relativa de los isotopos 
radioactivos. 

En el disco galactico, la densidad del medio interestelar es p ~ rripUd, donde es la 
masa del proton j rid — lcm~^; entonces, si los rayos cosmicos estan confinados a la region 
del disco, una columna de densidad ^ ~ 50 kg m~^ corresponde a Xg ^ 3 x 10® afios. Notemos 
que este tiempo es mucho mayor a t = c Lq ~ 3 x (10^ — 10"^) afios, el tiempo que tardan'a 
una particula relativista en atravesar la galaxia, tomando Lq ^ 1-10 kpc como valores 
indicativos de la dimension de la galaxia. Esto implica que los rayos cosmicos de bajas 
energias estan efectivamente confinados por los campos magneticos galacticos. 

Una determinacion mas precisa del tiempo de escape puede obtenerse a partir de la 
abundancia del isotopo radioactive del berilio, ^°Be, cuya vida media es ^ 3.9 x 10® afios, 
es decir, similar a la primera estimacion de Tg que hemos hecho. El berilio es producido por 
la fragmentacion de carbono y oxfgeno; el 10% del total producido es la especie radioactiva 
^°Be, que luego decae a por decaimiento /?, mientras que el resto corresponde a las 
especies estables ''Be y ^Be. En el modelo estacionario leaky box, las densidades vienen 
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dadas por 

------ + a = o, (1.14) 

r- T- r- 

2 I I 

donde, para simplificar la notacion, definimos como la tasa de produccion de la especie 
i— esima: 

j>i i 

Asumiendo que los isotopos de berilio son producidos por spallation de los elementos del 
grupo M (es decir, C, N y O), las densidades resultan 

n 

Ni = —, ; (1.16) 

para los isotopos estables, esta expresion requiere tomar el Ifmite rf ^ oo. El valor observado 
de la abundancia del isotopo inestable ^°Be relativa a la abundancia total del Be producido 
es 0.028 [47]; reemplazando este valor en la solucion dada por la ec. ()1.16p . resulta que 
el tiempo de escape es tipicamente Te ~ 10^ aiios, correspondiendo a una densidad media 
n ^ 0.2 cm~^. Esto implica que las particulas pasan mucho tiempo en el halo que rodea la 
galaxia, cuya densidad (en el caso de un halo extendido de ~ 10 kpc de altura) es tipicamente 
n ~ 10~^ cm~^. El tiempo de escape decrece con la energi'a de las particulas, dado que el 
confinamiento debido a los campos magneticos galacticos se torna menos efectivo a energias 
mayores. De la ec. ()1.13|) estacionaria, resulta que Tg ~ L"^ /D, donde L es una dimension 
caracten'stica de la region de propagacion; como veremos en el proximo Capitulo, la difusion 
a energias por debajo de la rodilla tiene una dependencia con la energia D oc E^^^, de donde 
resulta que Xg oc E~^^^. 



1.5 La deteccion de rayos cosmicos 

La deteccion de los rayos cosmicos de alta energfa {E > 10^^ eV) es indirecta, y requiere 
reconstruir la cascada de particulas secundarias que produce en la atmosfera la incidencia 
de la particula primaria. 

Consideremos, en primer lugar, la incidencia de un proton de alta energia sobre un 
nucleo atmosferico de masa Amp. Como el radio nuclear, dado por R ^ 1.2 A^/^ fm (1 fm 
= 10~^^ m), es mucho mayor que el proton relativista (por ejemplo, la longitud de onda 
de de Broglie es ~ 0.01 fm para un proton de 10 GeV), el proton realiza un scattering 
multiple con los nucleones individuales en el interior del nucleo. En cada interaccion, se 
producen principalmente piones y kaones, emitidos fundamentalmente en la direccion de 
incidencia; a la vez, estas particulas secundarias pueden tambien iniciar otras colisiones 
dentro del nucleo. Cada nucleon que colisiona con el proton incidente o con las particulas 
secundarias es, por lo general, removido del nucleo; si el nucleo resultante es de una especie 
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inestable, se fragmentara en niicleos estables mas livianos. De cada fragmento excitado, 
ademas, pueden emitirse neutrones. Una colision entre un proton de energia E > 1 GeV y 
un nucleo atmosferico produce tipicamente 2{E/GeVy^^ particulas relativistas cargadas. 

En el estudio de Uuvias atmosfericas, suele utilizarse come variable al slant depth o 
grammage, que representa la coordenada longitudinal X que mide la columna de densidad 
atravesada desde la parte superior de la atmosfera. La profundidad atmosferica vertical X^ 
viene dada por 

poo 

X^= p{h')dh' ^X COS 9 , (1.17) 

Jh 

donde p{h') es el perfil de densidad de la atmosfera, y ^ es el angulo de inclinacion de la 
Uuvia, medido desde la vertical. El camino libre medio de un proton de alta energia en la 
atmosfera es Xp_atm ~ 80 g cm~^, valor mucho menor que la columna de densidad total de 
la atmosfera en la superficie al nivel del mar, Xg^p ^ 1000 g cm~^. 

Si la particula primaria es un nucleo, en lugar de un proton, la colision con los niicleos 
atmosfericos puede generar fragmentos relativistas que desarroUan luego cascadas por sep- 
arado; en efecto, para E > 10^^ eV se observan EAS que tienen mas de un eje. Por otro 
lado, la interaccion del primario ocurre mas arriba en la atmosfera, debido a que la scccion 
eficaz aumenta con la masa del nucleo. Una regla sencilla, conocida como "principio de su- 
perposicion" , considera cada nucleo relativista de masa Arrip como un grupo de A nucleones 
independientes, ya que la energia de ligadura entre los nucleones puede despreciarse; en ese 
caso, la seccion eficaz de interaccion con los niicleos atmosfericos depende linealmente de la 
masa. 

Los nucleones secundarios, piones cargados y kaones continiian propagandose en la atmosfera, 
generando nuevas colisiones nucleares hasta llegar a la energia por particula umbral para la 
produccion de piones multiples, alrededor de 1 GeV. A cstc proccso sc lo llama cascada 
nucleonica o hadronica y, en el, la energia inicial del primario incidente se degrada en un 
gran niimero de hadrones de baja energia. Sin embargo, los piones y kaones tambien pueden 
decaer antes dc intcractuar. En particular, los piones neutros tienen una vida media muy 
corta (r = 1.78 x 10^"*^^ s) y dccacn rapidamcntc via 7r° 27. Los fotoncs resultantcs 
producen pares electron-positron, y estos producen, a la vez, fotones via bremsstrahlung; el 
camino libre medio de ambos procesos en la atmosfera es similar, Xem ~ 280 m. Esta cas- 
cada de fotones, electrones y positrones que se reconvierten sucesivamente entre si recibe el 
nombre de lluvia electromagnetica; tipicamente, exhibe un crecimiento inicial exponencial, 
desarrolla un mimero maximo de particulas proporcional a la energia de la particula que 
inicia la cascada y, mas alia del maximo, se atenua rapidamente. Cuando la energia por 
particula cae por debajo de la energia critica Ec (que es aquella para la cual las perdidas por 
bremsstrahlung igualan a las perdidas por ionizacion), los electrones y positrones pierden 
energia en la produccion de iones y electrones libres; para electrones en aire, E^ = 83 MeV. 
Con rcspccto a los fotones, por otra parte, la seccion eficaz de produccion dc pares tambien 
decrece y se hace del orden de las correspondientes al scattering Compton y absorcion fo- 
toelectrica. Los piones cargados tambien pueden decaer, via las reacciones ir'^ ^ /i'^ + y 
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7r~ ^ /x~ + P^, cuya vida media es r = 2.55 x 10~^ s. Los muoncs, a su vez, pueden decaer 
via /i"*" — i> e"*" + + i^/i y — ^ e~ + Pg + ^/x (con r = 2.2 x 10~^ s) e iniciar otras Uuvias 
electromagneticas. 

Como caracteristicas generales, diremos que una Uuvia atmosferica se desarroUa, alcanza 
un maximo de Nm particulas a una profundidad Xm J luego se ateniia, debido fundamen- 
talmente a perdidas por ionizacion. Cuanto mayor cs Eq, la cncrgia dc la particula primaria, 
mayores resultan y una estimacion (prccisa dcntro del ~ 25%) cs Eq ^ 1.4 Nm GeV. 
El frente de una lluvia de particulas ticnc un cspcsor tipicamente no mayor a algunos metros, 
y una extension que depende de Eq, por ejemplo, las Uuvias verticales producidas por rayos 
cosmicos ultra-energeticos tienen su maximo desarroUo a nivel del mar, involucran ~ 10^^ 
particulas y su extension lateral es de varios km. Hasta cl maximo, la mayor parte de la 
energia de la lluvia esta contcnida cn los clectrones (y positroncs); luego, csta componentc se 
atemia rapidamente, y la encrgia pasa a estar contenida fundamentalmente en los muones. 

La tarea de extraer la energfa y masa de la particula primaria a partir de los parametros 
observados en la lluvia atmosferica (el numero total de electrones, muones y hadrones, la 
forma de sus distribuciones latcralcs dc densidad, la reconstruccion de la altura del maximo 
de la lluvia, etc.) requiere adoptar modelos que describan el desarrollo de las cascadas; 
con ese fin, se ban desarroUado varios programas de simulacion Monte Carlo como, por 
ejemplo, CORSIKA y AIRES. Los modelos QCD fenomenologicos que se incorporan a los 
programas de simulacion (por ejemplo, QJSJET, VENUS y SYBILL) quedan constreiiidos 
por el gran cuerpo de datos cxpcrimentales disponibles a partir de los experimentos de 
colisiones de iones pesados (en el CERN y Brookhaven) y de colisiones proton-antiproton 
(en el CERN SPS y el Tevatron del Fermilab). Sin embargo, las interacciones que resultan 
de rayos cosmicos con energias por encima de la rodilla involucran energias de centro de 
masa mayores a las disponibles en el Tevatron; en consecuencia, los modelos de interaccion 
hadronica que se incorporan a los codigos dc simulacion dcbcn cxtrapolar los rcsultados 
experiment ales conocidos por varios ordenes de magnitud, y esto introduce grandes inccrtezas 
en los resultados. En este contexto, los experimentos de rayos cosmicos pueden ofrcccr la 
posibilidad de testear las simulaciones de Uuvias atmosfericas y, en consecuencia, los modelos 
de interaccion hadronica en regiones no accesibles a los laboratorios. 

Como hemos comentado mas arriba, en la fisica de rayos cosmicos se ha desarroUado, a 
lo largo de casi un siglo de investigaciones, un numeroso y variado conjunto de experimentos; 
a modo de ejemplo, aqui dcscribircmos brcvcmcntc cl cxpcrimcnto KASCADE {Karlsruhe 
Shower Core and Array Detector), en operacion desde 1996 en Karlsruhe, Alemania. Un 
arreglo de centelladores de 200 x 200 m^ mide las componentes electromagnetica y muonica. 
Un sistema central de detectores, de 320 m^, combina un gran calorimetro hadronico con 
varios detectores de muones, mientras que, adicionalmente, las trazas de muones de alta 
energia son medidas por un sistema subterraneo de 128 m^. Estc cxperimcnto combinado 
permite medir simultaneamente las cascadas electromagnetica, muonica y hadronica; ademas 
de ser, probablemente, el experimento que provee las medidas mas precisas y confiables en 
el rango £■ ~ 10^^ — 10^^ eV, la informacion "redundante" que es capaz de proveer presenta 
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la posibilidad de poner a prueba los modelos fenomenologicos de interaccion hadronica en 
procesos de muy alta energi'a. Como parte de las perspectivas futuras del experimento, se 
proyecta su extension para permitir la deteccion de Uuvias atmosfericas para energias del 
primario de hasta 10^^ eV, en un proyecto conocido como KASCADE-Grande [48,49]. 



Capitulo 2 



La difusion de rayos cosmicos en la 
galaxia 

En este Capitulo, cstudiamos cn dctalle la propagacion de rayos cosmicos en 
los campos magneticos turbulentos de la galaxia, teniendo en cuenta su caracter 
esencialmente difusivo. En este escenario, la apropiada inclusion de los efectos 
de drift macroscopico provee una explicacion natural de las observaciones del 
espectro, la coinposicion y las anisotropias (amplitud y fase en ascencion recta) 
de los rayos cosmicos en el rango E ^ 10^^ — 10^^ eV. En efecto, la rodilla y la 
segunda rodilla del espectro resultan del crossover entre el regimen dominado 
por la difusion transversal (a bajas energfas) y el dominado por los drifts (a en- 
ergias altas); veremos que la rodilla aparece como consecuencia de la supresion 
de los protones, mientras que la segunda rodilla resulta de la supresion de los 
nucleos de hierro. Por otra parte, la tendencia hacia componentes mas pesadas 
en la composicion, y la dcpcndcncia con la cncrgia de las anisotropias, repro- 
ducen adecuadamcntc las observaciones mas recientes. Finalmente, se presenta 
un estudio del tensor de difusion en condiciones de alta turbulencia, relevante 
en el contexto del transporte de rayos cosmicos galacticos de muy alta energia, 
aunque tambien util en otras aplicaciones de interes astrofisico. 

2.1 El escenario de la difusion turbulenta y drift 

2.1.1 Las componentes regular y turbulenta 
del campo magnetico galactico 

Como otras galaxias espirales, la Via Lactea tiene un campo magnetico regular de gran 
escala Bo(x) en el que las Imeas de campo tienden a seguir la estructura de los brazos 
espirales. Este campo se observa haciendo uso del efecto Faraday, es decir, midiendo la 
rotacion en el piano de polarizacion de la radiacion producida en piilsares y fuentes de radio 
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Figura 2.1: Esquema del campo regular del disco de la Via Lactea. El sol, ubicado sobre 
el piano galactico, a una distancia de 8.5 kpc del centro, aparece indicado en coordenadas 
galacticas (de [52]). 

extragalacticas [50,51]. Diferentes estructuras del campo magnetico regular de galaxias 
espirales son descriptas por una variedad de modelos. Una primera clasificacion distingue 
entre campos bisimetricos (BSS) y axisimetricos (ASS), definidos de modo que los campos 
BSS (ASS) son antisimetricos (simetricos) con respecto a una rotacion de tt radianes en torno 
de un eje normal al piano galactico que pasa por el centro. Por ejemplo, en una galaxia 
espiral con dos brazos, el campo invierte su signo entre ellos en un modelo BSS, mientras 
que mantiene su orientacion en un modelo ASS. En el caso particular de la Via Lactea, el 
campo magnetico regular en el disco galactico es predominantemente toroidal, alineado con 
la estructura espiral y con inversiones de orientacion entre brazos contiguos, de modo que 
parece corresponder a un modelo BSS. La fig J2.1l muestra esquematicamente una estructura 
BSS que podria corresponder al campo regular del disco de la Via Lactea, e indica tambien la 
posicion del sol en coordenadas galacticas (ubicado sobre el piano galactico, a una distancia 
de 8.5 kpc del centro). Adicionalmente, los modelos de campo magnetico pueden clasificarse 
como simetricos (S) o antisimetricos (A) con respecto al piano galactico; en particular, el 
campo magnetico regular en el disco galactico de la Via Lactea es usualmente considerado 
como simetrico con respecto al piano galactico, de modo que corresponderia a un modelo 
BSS-S. 

La deteccion de una significativa emision de ondas de radio no termicas a grandes lat- 
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itudes galacticas indicaria la existencia de una corona con campos magneticos formando 
un halo, que podria extenderse posiblemente hasta una distancia de 10 kpc desde el piano 
galactico. Aunque en este halo magnetico la orientacion del campo regular no es conocida, 
parece ser dominante la componente toroidal. La estructura del campo magnetico regular 
del halo podria ser similar a la del campo en el disco galactico (por ejemplo, si su origen esta 
relacionado con un viento galactico que de algiin modo extiende hacia el halo las propiedades 
del campo del disco), pero tambien podria ser diferente (por ejemplo, considerando que es 
producida en forma independiente mediante un mecanismo de dinamo ^). 

Considerando las principales caracteristicas del campo magnetico regular que hemos men- 
cionado anteriormente, aqui asumiremos que la region de propagacion de los rayos cosmicos 
en la galaxia es un cilindro de radio R = 20 kpc y altura 2H, con un halo de tamano 
H = 10 kpc. El angulo que forman los brazos espirales respecto de la direccion azimutal 
es pequeno (menor a 10°), de modo que es una buena aproximacion suponer que el cam- 
po magnetico regular esta orientado segun la direccion azimutal. En consecuencia, resulta 
natural adoptar coordenadas cilfndricas (r, 0,2;), tomando Bq = -Bo0 y asumiendo por sim- 
plicidad que el sistema posee simetria azimutal. Puede notarse, ademas, que un campo 
azimutal independiente de la coordenada tiene divergencia nula, de manera que la con- 
sistencia con las ecuaciones de Maxwell esta asegurada (equivalentemente, desde un punto 
de vista mas fisico, las Imeas de campo son circunferencias cerradas, de modo que no hay 
fuentes ni sumideros de Imeas de campo). 

Para reproducir las inversiones de orientacion del campo propias del modelo BSS-S, 
adoptaremos para la componente regular del disco galactico la expresion 

Bt-(r, .) ^ B.^^j^.n th= (^) -±^^ . (2.1) 

donde Tc = 4 kpc representa un radio caracteristico que suaviza el comportamiento radial 
global 1/r en la region del centro galactico, es una altura caracterfstica del disco galactico, 
y Bd es una constante que fija el valor del campo en nuestro punto de observacion {vobs = 
8.5 kpc, Zobs = 0), que es negativo (ya que el campo magnetico galactico local esta dirigido 
aproximadamente hacia la longitud / ~ 90°, es decir, en la direccion de —0). El factor th^(r) 
ha sido introducido solo para asegurar que no haya singularidades en el eje r = 0, dado que, 
de lo contrario, apareceria una corriente de deriva (drift) singular espuria a lo largo de ese eje 
(ver la ec. ()2.15|) mas adelante). Ademas, puede notarse que el factor sinusoidal introduce 
inversiones de orientacion que ocurren para r = 4, 8, 12 y 16 kpc. 

La componente regular correspondiente al halo puede describirse mediante una estructura 
de campo con inversiones radiales y simetrica con respecto al piano galactico (que Uamare- 

"'^El mecanismo de dmamo, sugerido originalmente por Larmer en 1919 [53], explica el origen de estructuras 
de campo magnetico a gran escala como resultado de una realimentacion constructiva en el acoplamiento entre 
el campo magnetico y el campo de velocidades en un flui'do magnetohidrodinamico turbulento. Bajo ciertas 
geometrias particulares, las corrientes inducen campos magneticos que, a su vez, realimentan y amplifican 
las corrientes iniciales, hasta alcanzar un estado de saturacion debido a mecanismos disipativos. 
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mos R — S), o bien por un modelo sin inversiones radiales (que llamaremos NR); en este 
ultimo caso, cuatro diferentes configuraciones resultan de considerar la orientacion relativa 
del campo en el halo con respecto a la del disco, y de la simetria con respecto al piano 
galactico, que denotaremos con NR — S± y NR — A± (donde el signo ± es el del producto 
j^haio . j^disco gyaluado en r = Tobs Y z > 0). Debe notarse que, al considerar un halo con 
inversiones radiales, estamos suponiendo que el halo tiene su origen en un fenomeno asoci- 
ado al disco galactico (por ejemplo mediante un viento galactico, como hemos mencionado 
anteriormente) , de modo que la estructura del campo magnetico deberia ser similar. Esto 
justifica que, en este caso, sea suficiente considerar solo la configuracion simetrica R — S, 
mientras que un halo sin inversiones radiales podri'a tener un origen independiente del campo 
del disco y requiere el estudio de varias configuraciones diferentes. Recientes estudios que 
utilizan el efecto Faraday podrian indicar que la estructura regular del halo magnetico es 
la de un campo azimutal con orientacion invertida arriba y debajo del piano galactico, que 
podria corresponder a una configuracion de un dmamo AO [54]. Estos resultados podrian 
sugerir que el modelo NR — A es el mas apropiado, aunque debe tenerse en cuenta que la 
evidencia observacional no esta aiin suficientemente bien establecida, de modo que todas las 
configuraciones que hemos descripto resultan plausibles. 

Si la componente regular del halo tiene inversiones radiales, consideraremos que estara 
dada por 

B^'^'^ir, z) = B. ^ , ^ sin ( th^ ( , (2.2) 

mientras que si su estructura no tiene inversiones radiales, asumiremos 

5„^«'°(r,z) = Bh , ^ = , ^ =th' ( ^— ^ 7^ , (2.3) 

donde Zh es una altura caracteristica del halo, y Bh es una constante que fija el valor y 
la orientacion local del campo del halo. La funcion TZ es simplemente la unidad para el 
caso del halo simetrico (con respecto al piano galactico), mientras que es tomada como 
TZ{z) = ta.nh.{z / Za) en el caso de los modelos antisimetricos. 

Las ecs. (j2.1|) - ()2.3p proveen una descripcion simple pero realista para la estructura 
regular a gran escala de los campos magneticos de la galaxia, y contienen las caracteristicas 
esenciales (como, por ejemplo, las inversiones radiales de orientacion) relevantes para la 
determinacion de los efectos de drift en la propagacion de los rayos cosmicos. 

Superpuesto a las estructuras regulares del campo magnetico del disco galactico y el 
halo extendi do, existe un campo magnetico irregular Br(x) estocastico [random) produci- 
do por la turbulencia del plasma interestelar. La escala maxima de las irregularidades es 
Lmax — 100 pc, y localmente su intensidad es comparable a la de la componente galactica 
regular, (i?j.)/(|Bo|) ~ 1 - 3. Los datos observacionales [55] muestran que el espectro de 
inhomogeneidades en la densidad del gas interestelar podria ser el mismo que el correspon- 
diente a la componente irregular del campo magnetico, y esto corresponde a un espectro 
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de Kolmogorov [56]. En consecuencia, las componentes de Fourier del campo magnetico 
random contienen una densidad de energi'a magnetica dEr/dk oc k~^^^ (con k > 27i/Lmax)- 
Asumiremos que la intensidad de la componente irregular posee una dependencia espacial 
suave, analoga a la envolvente asumida para los campos regulares, dada por 



donde Zr es la escala caracteristica del perfil vertical y B* la amplitud cuadratica media 
(rms) local. 

En el contexto de este trabajo, estamos interesados en describir la propagacion de parti'culas 
ultrarrelativistas, cuya velocidad f ~ c es mucho mayor que la velocidad de Alfven {va — 
5 X 10^ m/s) caracteristica de la propagacion de las ondas magneticas en el medio inter- 
estelar ^. En consecuencia, esta bien justificada la suposicion de considerar perturbaciones 
magneticas estaticas, despreciando los campos electricos. 

2.1.2 La ecuacion de difusion 

Una particula de carga Ze, propagandose en un campo magnetico regular uniforme Bq, 
describe una trayectoria lielicoidal, caracterizada por el angulo de pitch 6 (definido como el 
angulo entre la velocidad instantanea y la direccion del campo) y el radio de Larmor, dado 
por 



donde la ultima expresion es valida en el li'mite relativista. La componente de la velocidad 
paralela a Bq es v\\ = c cos9, mientras que el radio de la helice es ^ a = rz, sen6'. La rigidez 
magnetica de la particula se define como 



para una particula relativista, R = E[eV]/Z Volts. 

^En magnetohidrodinamica, se considera la propagacion de tres tipos de movimiento ondulatorio: las 
ondas de sonido, que son longitudinales y se propagan en la direccion del campo con velocidad Vg = y^kT/rrii, 
donde T es la temperatura de los electrones termicos y rrii la masa de los iones en el plasma; las ondas 
de Alfven, que son transversales y se propagan en la direccion del campo con velocidad va — B/^/Anp, 
donde p es la densidad del plasma; y las ondas magnetosonicas, que son longitudinales y se propagan 
perpendicularmente al campo con velocidad y^i)f + iij. En las condiciones del medio interestelar, es habitual 
asumir que Vs <SC va [3], de modo que la velocidad de Alfven es caracteristica de la propagacion de las 
perturbaciones magnetohidrodinamicas. 

•^Otra definicion del radio de Larmor resulta de reemplazar p p± = p sen0 en la ec. H2.5(l . de modo 
que a = r^; esta eleccion es la usual en los libros de texto sobre electromagnetismo (por ej., Jackson [57], p. 




(2.4) 




(2.5) 
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En presencia de un campo magnetico random Br, las particulas se dispersan en las 
irregularidades del campo estocastico y cambian el angulo de pitch; en este proceso, conocido 
como "dispersion del angulo de pitch" {pitch angle scattering), la velocidad, la rigidez y el 
radio de Larmor ()2.5p permanecen constantes. La dispersion del angulo de pitch ocurre 
con probabilidad maxima en condiciones de resonancia (es decir, cuando la escala de las 
irregularidades es del orden de r^), de modo que es un mecanismo de isotropizacion eficiente 
siempre que < L^ax- Por ejemplo, en el campo magnetico galactico, caracterizado por 
una amplitud |Bo| de algunos /xG y una escala maxima de turbulencia L^ax — 100 pc, 
la dispersion del angulo de pitch es eficiente para protones hasta energias por encima de 
10^^ eV. 

Por consiguiente, para energias compatibles con la condicion de resonancia, ri < Lmax, 
el transporte de rayos cosmicos a traves de los campos magneticos galacticos turbulentos 
queda bien descripto por medio de una ecuacion de difusion, en la cual las componentes 
del tensor de difusion dependen tanto del campo regular como de la turbulencia estocastica. 
Despreciando la fragmentacion nuclear, las perdidas de energfa de cualquier tipo (por ejem- 
plo, por ionizacion, bremsstrahlung, emision sincrotronica, etc.), efectos convectivos (que 
podrian deberse, por ejemplo, a un posible viento galactico), procesos de reaceleracion, etc., 
la ecuacion estacionaria de transporte p.9p resulta 

V ■ J = g(x) , (2.7) 

donde Q es el termino que describe a las fuentes, y donde la corriente de rayos cosmicos, 
dada en terminos de la densidad A^(x), es 

Ji = -Ai(x)VjA^(x) . (2.8) 

En general, el tensor de difusion Dij puede escribirse como 

Ai = {Dii - Dj_) b^bj + D^6ij + DAe^jkbk , (2.9) 

donde b = Bo/|Bo| es un vector unitario a lo largo del campo magnetico regular, 6ij es la 
delta de Kronecker, y eijk es el tensor antisimetrico de Levi-Civita. Los terminos simetricos 
contienen los coeficientes de difusion que corresponden a las direcciones paralela y perpendic- 
ular al campo regular {D\\ y D±) y describen la difusion debida a las fluctuaciones turbulentas 
de pequena escala. El termino antisimetrico contiene el coeficiente de difusion de Hall (-Da) 
y es el responsable de la corriente macroscopica de drift ^ de las particulas en propagacion. 
La ecuacion de difusion puede ser reescrita como 

-V^Ds^JVjN + UiViN = Q , (2.10) 

^La estructura del tensor de difusion, ec. H2.9|l . esta explicada en el Apendice A. Aqui asumimos Da > 
para la difusion de particulas con carga positiva, mientras que las de carga negativa tienen asociado el signo 
opuesto. 

^El drift macroscopico resulta del comportamiento colectivo de una distribucion de particulas, y no debe 
ser confundido con el drift microscopico, asociado a la curvatura y los gradientes en Bq. 
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donde Dsij = {D\\ — D±^ bibj + D±Sij contiene la parte simetrica del tensor de difusion, 
mientras que Ui = —eijkVk (DAbj) es la velocidad de drift. Por otro lado, la corriente de 
rayos cosmicos dada por la ec. ()2.8|) puede tambien reescribirse como 



-D^_\/±N - Z^ii ViiA^ + DAh X ViV 



(2.11) 



donde Vy = b(b ■ V) y V± = V — Vy. En esta expresion se observa explicitamente que la 
corriente macroscopica de drift (dada por el ultimo termino de la ec. 1)2.111) ) es ortogonal 
tanto a la direccion del campo regular como al gradiente de densidad de rayos cosmicos. 

Bajo la suposicion de que el campo magnetico regular esta dirigido en la direccion azimu- 
tal (es decir, br = = 0, b^i, = ±1) y de que el sistema posee simetria azimutal {d/d(j) = 0), 
result a V±N = ViV y la corriente macroscopica de rayos cosmicos queda dada por 



-D^VN + DAh X VN 



(2.12) 



es decir, desaparece el termino que involucra a D\\ y la corriente resulta perpendicular al 
campo magnetico regular. Dicho en otros terminos, las simetrias asumidas excluyen la 
difusion en la direccion del campo (o sea, en la direccion que corresponde a los brazos 
espirales). Reemplazando ()2.12|) en ()2.7|) . la ecuacion de difusion en coordenadas cilmdricas 
resulta 
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En estas expresiones, vemos que un cambio de signo en b^ (como resulta de una inversion en 
la orientacion de los campos regulares) afecta las velocidades de drift y, por consiguiente, la 
difusion en esa region. 

Sobre la superficie exterior S correspondiente a los bordes del halo, la densidad de rayos 
cosmicos debe anularse, dado que las partfculas pueden escapar libremente hacia el espacio 
intergalactico. Por lo tanto, las ecs. ()2.13|) - ()2.15p deben considerarse con la condicion de 
borde A^|s = 0. 

Con respecto a la distribucion espacial de las fuentes, asumiremos que estan contenidas 
en un disco delgado de altura 2hs = 400 pc, de modo que el termino de las fuentes puede 
escribirse como Q{r,z) = q{r)6{hs — \z\), donde 6 es la funcion de Heaviside y q{r) es una 
funcion de distribucion radial. Como se describe mas adelante, en la Seccion 2.1.4, la ecuacion 
de difusion puede resolverse numericamente mediante la discretizacion de Crank-Nicolson y 
el metodo de Liebmann [58]. 
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2.1.3 Calculo de los coeficientes de difusion 

La difusion en la direccion del campo magnetico regular, debida a la dispersion del angulo 
de pitch, conduce a un coeficiente de difusion dado por 

^11 = |A|| , (2.16) 

dondc A II cs el camino libre medio en la direccion paralela al campo [59]. En esta expresion, 
A|| depende de la densidad de energia asociada a los modos del campo magnetico random en 
escalas del orden de r^, es decir. 



All oc 



dEr/d Ink 



(2.17) 

fc=27r/r^ 



donde dEr/dk es el espectro de Fourier de la densidad de energia contenida en el campo 
magnetico random. 

En general, una dada componente del tensor de difusion esta asociada a ciertos procesos 
fisicos particulares que pueden no contribuir a las otras componentes. La difusion transversal 
al campo magnetico regular es debida tanto a la dispersion del angulo de pitch (es decir, 
el mccanismo que prcvalccc cn cl case dc la difusion paralela) como a la propia difusion de 
las Imeas de campo, que arrastran con cUas a las particulas cn la direccion perpendicular 
a Bo (un proceso conocido como field line random walk). Cuando el nivel de turbulencia 
es pequeiio, la difusion en la direccion transversal es mucho mas lenta que en la direccion 
paralela; en cambio, en cl limite de turbulencia muy alta, la anisotropia debida al campo 
regular ticndc a dcsaparecer, y los movimientos paralelo y perpendicular resultan similares. 

Sin embargo, no existe un tratamiento teorico general de este problema que permita, una 
vez definidos los campos, la determinacion de los coeficientes de difusion. Considerando el 
limite de baja turbulencia, que permite un tratamiento perturbativo del campo magnetico 
random, los coeficientes de difusion resultan 

D, = ^^uir, , (2.18) 



-3 1 + {u;t±) 



Da^-^T^^^^.^ (2.20) 



donde ~ c es la velocidad dc las particulas, = f /r^^ cs la frecuencia angular de Larmor, 
y T\\ = T± = ta = T representa un tiempo caracteristico de descorrelacion de las trayectorias. 
Debido a que en este tipo de analisis se asume un linico mecanismo de scattering como 
responsable de todas las descorrelaciones, interviene una linica escala temporal r para los 
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tres coeficientes de difusion. Esto ocurre, por ejemplo, cuando se considera la evolucion de 
una distribucion de parti'culas casi isotropica bajo un proceso de dispersion homogeneo, en 
el que r representa el tiempo caracterfstico de isotropizacion [60]. Una estructura similar 
se obtiene tambien para los coeficientes de difusion asociados a procesos de colision en un 
plasma en equilibrio termico [61], o bien, considerando r como un tiempo caracterfstico de 
dispersion y tratando las perturbaciones magneticas como centros dispersores constitufdos 
por esferas rfgidas [62], en lo que se conoce como el "resultado de dispersion clasico" [63-65]. 

En condiciones de alta turbulencia, en cambio, el tratamiento perturbativo no puede 
aplicarse, y debemos recurrir a un enfoque analftico de otra naturaleza [66], o bien a simu- 
laciones numericas [67,68]. Debe notarse que el problema de la difusion de rayos cosmicos 
en la galaxia requiere tratar con niveles de alta turbulencia debido a que, por un lado, la 
intensidad del campo magnetico turbulento es del mismo orden que la del campo regular 
y, por otra parte, a que en las regiones en las que el campo regular tiene inversiones de 
orientacion su amplitud es pequefia, de modo que allf prevalece largamente la componente 
estocastica. 

Utilizando el formalismo de Kubo [66,69,99], los coeficientes de difusion pueden 
obtenerse, alternativamente, de la descorrelacion entre diferentes componentes de la veloci- 
dad de una partfcula, promediada sobre una distribucion isotropica de muchas partfculas y 
sobre muchas configuraciones de campo random diferentes. En efecto, resulta^ 

poo 

Dij = dt Rij{t) , (2.21) 
Jo 

con Rijit) = {voiVj(t)) (donde (...) corresponde al promedio sobre distribuciones y configura- 
ciones de campo, y con Vi{t = 0) = Voi). En [66], las funciones de descorrelacion se obtienen 
a partir de las expresiones de Rij correspondientes a trayectorias helicoidales, agregando 
factores ad hoc de modulacion exponencial para representar el efecto de la turbulencia. Es 
decir, asumiendo que el campo regular esta dirigido segiin z, se adoptan las expresiones 

Rxxit) = Ryyit) = - coscjt e-*/"^ , (2.22) 



-Ry^{t) = — sinujt e"*/^^ (2.23) 

R.M = y e-*/-ii , (2.24) 

donde las tres escalas temporales de descorrelacion, t\\,t± y ta, son en principio diferentes. 
Si bien este enfoque pretende ser general, y valido aun en condiciones de alta turbulen- 
cia [66], debe notarse que las expresiones ()2.22|1 - ()2.24|1 estan fuertemente inspiradas en el 
comportamiento que se observa a baja turbulencia; el rango de validez de este enfoque no 

^La formula de Kubo aparece discutida en el Apendice B. 



Rxy(t) 

y 
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es evidente de inmediato, y sera reconsiderado mas adelante, en la Seccion 2.4. Sin em- 
bargo, resulta interesante observar que, reemplazando las ecs. fl2.22|) - ()2.24|) para Rij en la 
integral ()2.2H) . se reobtienen los coeficientes de difusion dados por las ecs. ()2.18|) - ()2.20|) . 
donde D\\ = Dzz,D± = D^x = Dyy y Da = D^y = —Dy^. Es decir que este enfoque, que 
incorpora diferentes escalas de descorrelacion, puede ser visto como una generalizacion de 
los resultados perturbativos. 

Para poder aplicar las ecs. ()2.18|) - ()2.20p al calculo de los coeficientes de difusion en un 
problema concreto, todavia debemos especificar de algiin modo los tiempos de descorrelacion. 
Una alternativa es recurrir a hipotesis adicionales que simplifiquen el problema, vinculando 
las escalas temporales de descorrelacion. Por ejemplo, en [66] se asume que el tiempo de 
descorrelacion asociado a la difusion perpendicular es el mismo que el correspondiente a 
la difusion antisimetrica. Asumiendo, ademas, que la contribucion principal a la difusion 
perpendicular es la debida a la caminata aleatoria de las Imeas de campo, se Uega a la 
expresion 

2 rr 

UJT^ = --i- , (2.25) 

<J ^ flrw 

que relaciona la escala temporal de descorrelacion t±_ con el coeficiente de difusion Djirw 
asociado a la caminata aleatoria de las Imeas de campo. Luego, estimando Dfirw analitica o 
numericamente en funcion del nivel de turbulencia y de la forma del espectro de turbulencia 
asumido, se pueden calcular los coeficientes D± y Da- Sin embargo, la ec. ()2.25|) conduce a 
discrepancias evidentes con los resultados de [67], obtenidos mediante simulaciones extensivas 
Monte Carlo, sugiriendo que las hipotesis adicionales asumidas en [66] no son completamente 
apropiadas. 

Recientemente, Casse et al. [68] extendieron los resultados numericos de [67] y determi- 
naron la difusion paralela y transversal de particulas de prueba relativistas propagandose en 
campos estocasticos de Kolmogorov de gran intensidad. En este estudio se observa que, por 
lo general, la extrapolacion de los resultados perturbativos es inapropiada para describir la di- 
fusion de particulas en campos altamente turbulentos. Por ejemplo, aunque la ec. (j2.18|) para 
la difusion paralela parece concordar bien con los resultados obtenidos en [68], la ec. ()2.19p 
para la difusion transversal (con ry = t±) es claramente inapropiada, como es de esperar del 
hecho de que, en general, la caminata aleatoria de las Imeas de campo no puede despreciarse. 

En este contexto, lo que haremos es extraer directamente el coeficiente de difusion 
transversal de los resultados numericos de [68] . Estos resultados son aplicables en el regimen 
de interes para el estudio de la region de la rodilla, ya que el parametro de energia relevante, 
definido como p = VL/Lmax, toma el valor p ~ 0.02 para protones cerca de la rodilla, que 
esta contenido en el rango 0.001 < 27rp < 10 que se explora en [68]. Definiendo el nivel de 
turbulencia como i] = B^/^Bq + B^), se encuentra que una expresion aproximada para D±, 
valida en el rango 0.1 < rj < 1, es 

D^ ~ 8.7 X 102^e3-24'7 (1 _ ^)i/6 ( !±] ''^ ^ . (2.26) 

Vpcy s 
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Aunque esta expresion se obtiene de simulaciones numericas, la dependencia con la rigidez a 
traves de una potencia 1/3 aparece tambien en expresiones teoricas, deducidas en el li'mite 
de baja turbulencia, bajo la suposicion de un espectro de Kolmogorov para la componente 
irregular del campo magnetico [59]. 

Para caracterizar los efectos de drift partiremos de la ec. (|2.2(Jj) . que permite obtener 
el valor de en casos extremos, tanto de como de muy baja turbulencia. En 

el caso de turbulencia muy alta {ujta 0) resulta Da — * 0, y los efectos de drift tienden 
a ser despreciables. Este Ifmite es de particular interes en las regiones con la inversion de 
orientacion del campo regular, en las que la turbulencia tiende a suprimir los efectos de drift. 
Por otra parte, en el li'mite de turbulencia muy baja {ujta ^ 1) la expresion obtenida tiende 
al resultado conocido de la aproximacion cuasi-lineal, es decir. Da = vri/S [59,66]. Con el 
fin de obtener una transicion suave entre los dos comportamientos li'mite, consideraremos las 
ecs. fl2.18|) - ()2.20|) con una escala caracterfstica temporal comiin para todos los procesos de 
descorrelacion ^. En particular, resulta que 



Este cociente, calculado en [68] en funcion de la rigidez y del nivel de turbulencia, es esencial- 
mente independiente de la energfa en el rango estudiado; se encuentra que la simple relacion 
aproximada ut ~ 4.5 ( | Bo| /-Br) ^'^j sustitui'da en la ec. ()2.27|1 . reproduce razonablemente los 
resultados numericos. Este resultado para lot puede ahora ser reemplazado en la ec. ()2.2Up 
(cambiando ta ^ r), j conduce a 



De este modo, se obtiene una expresion para el coeficiente de difusion de Hall que exhibe el 
comportamiento apropiado en los Ifmites de muy alta y de muy baja turbulencia. Tal como 
se muestra en el estudio numerico de la Seccion 2.4, la prescripcion ()2.28|) es adecuada en 
este contexto, y los resultados no cambian sustancialmente al adoptar una parametrizacion 
mas precisa. 

2.1.4 Resolucion numerica de la ecuacion de difusion 

Adoptando un modelo de campo magnetico galactico determinado, eligiendo valores partic- 
ulares para el conjunto de parametros que interviene en cada modelo, y especificando una 
funcion de distribucion radial para las fuentes, la difusion de rayos cosmicos se obtiene de 
resolver la ecuacion que resulta de las ecs. (j2.13|) - (j2.15j) con la condicion de borde A^ls = 0. 

^Mas arriba hemos hecho notar que r|| ^ t±. Sin embargo, con el fin de obtener una prescripcion que 
regule el comportamiento de Da (conectando el regimen de baja turbulencia con el de muy alta turbulencia), 
asumiremos, en lo que sigue, que todas las escalas temporales de descorrelacion son iguales. 




(2.27) 




(2.28) 
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La ecuacion diferencial queda definida una vez calculados los coeficientes de difusion y sus 
derivadas, que son determinados en funcion de la rigidez y del nivel de turbulencia mediante 
las ecs. (fTI^ y ^TI^ . 

Una solucion analitica de este problema queda naturalmente excluida, debido a la com- 
plejidad de las ecuaciones que deben resolverse. La resolucion numerica se llevo a cabo 
usando la discretizacion de Crank- Nicholson y aplicando el metodo de Liebmann [58], luego 
implementados en un programa FORTRAN. A continuacion, describiremos el esquema de 
discretizacion empleado en este trabajo. 

Comencemos por suponer que se busca resolver una ecuacion diferencial de la forma 

Cu = , (2.29) 

donde C es un operador diferencial que actiia sobre la funcion u{x), definida en el rango xa < 
X < xb, con condiciones de borde de Diriclilet u{xa) = ua, u{xb) = us- Si incorporamos a u 
una dependencia temporal, u{x) u{x,t), y planteamos la ecuacion diferencial dependiente 
del tiempo, 

^ = £m , (2.30) 

con una condicion inicial arbitraria (compatible con las condiciones de borde, que deben 
satisfacerse para todo t), podemos encontrar la solucion de la ec. ()2.29|1 como la solucion de 
equilibrio a la que relaja ()2.30|) para t oo. Este procedimiento es conocido como metodo 
de relajacion. 

Implementando la tecnica de diferencias finitas, discretizamos la variable espacial en un 
lattice regular de + 1 sitios, con = (x^ — xa)/A, mientras que al tiempo se lo discretiza 
con un mdice entero n > 0. En este punto, debe adoptarse un esquema de discretizacion 
adecuado para representar al operador diferencial £; la estabilidad del calculo numerico y 
la rapidez con que los resultados relajan a la solucion asintotica dependen, en general, del 
esquema elegido. En particular, para £ = d^/dx^, un esquema de discretizacion estable de 
la ec. (|2.3(J|) es el de Crank-Nicholson, dado por 

- < = % ■ (2.31) 

Resolver la ec. ()2.31|) implica la resolucion de un conjunto de A^ — 1 ecuaciones lineales 
simultaneas que involucran los uf^ con 0<i<NYm = n,n+l. Agrupando los terminos 
en (|2.31|) apropiadamente, resulta 

- (A2 + 2)<+i + = -A V (1 < ^ < AT - 1) , (2.32) 

un sistema tridiagonal que puede ser resuelto facilmente para cada paso temporal [58] . Pidi- 
endo que las diferencias entre dos configuraciones temporales sucesivas sean menores a una 
cota predeterminada, se contimia el esquema de calculo hasta alcanzar la convergencia. 
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Si C es un operador diferencial que puede escribirse como una suma de operadores, £ = 
Li (tales que, para cada uno por separado, se conoce un esquema de discretizacion estable 
que resuelve la ecuacion diferencial), puede seguirse el metodo de division de operadores, o 
division temporal, que indica actuar en forma sucesiva con cada operador £j. En particular, 
consideremos una ecuacion diferencial parcial lineal de segundo orden inhomogenea, 

d^u du ,du , , 

ox^ oy^ ox oy 

con condiciones de horde de Diriclilet o de Neumann. Esta ecuacion puede escribirse como 

Cu + ri = , (2.34) 

donde C = L\x + l^\y + l^2x + l^2y, con = d^/ da;^, = x d/ dx, y expresiones analogas 
para x ^ y. Generando el primer paso temporal mediante Cx^i resulta 



A2 A2 

^■i "'""2r'"" + V-., "^'""^"^^-" + -to , (2.35) 

que conduce a un sistema tridiagonal de ecuaciones. En los pasos temporales siguientes de 
la serie, se actualizan del mismo modo los terminos que corresponden a L^x Y l^2y 

La ec. ()2.33|) corresponde, en efecto, a la ecuacion de difusion ()2.13|) . identificando 
u[x, y) N{r, z), (j) ^ 1/r + {dD^/dr - Ur) /D^, if) [dD jdz - u ^) j t] ^ -Q, 
donde Ur y Uz son las velocidades de drift dadas por las ecs. ()2.14p y ()2.15|) . El pasaje de 
coordenadas cilmdricas a rectangulares presenta la dificultad de que el mapeo no es biyectivo 
en r = 0; la condicion de borde A^|s = que debe aplicarse sobre la superficie exterior S 
(correspondiente a los bordes del halo, de radio R y altura 2H) se traduce en condiciones 
de borde liomogeneas de Dirichlet en r = i? y 2; = ±if. Segiin el esquema de calculo que 
describimos, faltaria especificar la condicion de borde sobre r = 0, uno de los bordes del 
lattice en el piano r — z. Una solucion se obtiene mediante la extension a la region no fi'sica 
r < 0; extendiendo 0, y 77 al rango —R < r < R como funciones pares de la variable r, 
podemos imponer la condicion = a los bordes r = ±_R y el esquema de calculo queda 
bien definido. 



2.2 El espectro y la composicion de rayos cosmicos 
en la region de la rodilla 

En el contexto del modelo de difusion turbulenta que hemos descripto en la Seccion anterior, 
una explicacion para la rodilla del espectro de rayos cosmicos surge de un modo natural 
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Figura 2.2: Distribucion espacial de la densidad de rayos cosmicos en la galaxia en el modelo 
A^-R — S+ (halo sin inversiones radiales y simetrico con respecto al piano galactico), para 
una fuente en = 6 kpc: (a) a bajas energias (debajo de la rodilla), el regimen dominado 
por la difusion transversal; (b) a energias altas (arriba de la rodilla), el regimen dominado 
por los drifts. Los picos muy agudos que aparecen alrededor de la ubicacion de la fuente han 
sido truncados. 

como consecuencia del cambio de regimen entre la difusion transversal, dominante a bajas 
energias, y la difusion de Hall, que prevalece a energias altas ®. Como se observa de las 
ecs. (I2.26P y (j2.28j) . la dependencia de los coeficientes de difusion con la energi'a viene dada 
por D±{x.) ~ (x)(E/Eo)^/3 y L'a(x) ~ D\{-k)E/Eq. A bajas energias {E ^ Z Er), la 
difusion transversal es el mecanismo dominante que gobierna el transporte de rayos cosmicos, 
y conduce a dN/dE oc {D^/ Dj_)dQ/dE oc E~^~^/^, donde P es el indice espectral de 
la fuente, mientras que a energias altas {E ^ Z Er) se tiene en cambio que dN/dE cx 
{D\/DA)dQ/dE oc E~^~^. De este modo vemos que, si el indice espectral de los flujos 
producidos en las fuentes es /3 ~ 2.4, se obtiene la pendiente espectral correcta debajo de 
la rodilla {dN/dE oc i?"", con a ^ 2.7). Cada componente de rayos cosmicos de carga 
Z comienza a verse afectada por los drifts a la energia E Z Er, j su espectro se torna 
progresivamente mas empinado, con el indice espectral cambiando finalmente en Aa ^2/3 
en alrededor de una decada de energi'a. Veremos que la envolvente del espectro total, obtenido 
de la suma de las diferentes componentes nucleares, reproduce muy adecuadamente el cambio 
de un espectro oc E~'^-'^ debajo de la rodilla, a uno oc E~^ arriba de ella. Adicionalmente, 
como todas las componentes livianas quedan fuertemente suprimidas arriba de energias del 

^En el Apendice C se describe un modelo simplificado, cuya solucion anali'tica reproduce las caracteristicas 
esenciales de la rodilla del espectro. 
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Figura 2.3: Espectro diferencial normalizado versus E/Z para los diferentes modelos de 
campo magnetico regular investigados. R — S denota un halo con inversiones radiales y 
simetrico con respecto al piano galactico, mientras que NR — S{A) se refiere al halo sin 
inversiones radiales y simetrico (antisimetrico) con respecto a este piano. Ademas, en estos 
casos la orientacion del halo relativa a la componente regular del disco galactico se distingue 
mediante el signo de ■ B^, como se indica. 

orden de 10^^ eV, la componente pesada dominante (correspondiente a nucleos del grupo del 
hierro) cambia gradualmente su espectro hasta que el espectro total se torna oc £'~2-7-2/3 
encima de 10^^ eV. Asi, este escenario reproduce tambien la segunda rodilla observada en 
Egr — 4 X 10^^ eV y el espectro oc E~^'^ que se observa hasta el tobillo, tal como se describira 
en detalle mas adelante, en la Seccion 2.3. 

Con el objetivo de ilustrar el rol de los drifts en la difusion de rayos cosmicos, consideremos 
en primer lugar una fuente localizada en = 6 kpc, tomando para el campo magnetico 
random una intensidad rms local B* = 2 fiG y una escala vertical Zj. = 3 kpc. Las figs. 
I2.2r a) y (b) muestran la distribucion espacial de la densidad de rayos cosmicos en la galaxia, 
segiin resulta para una estructura de campo regular NR — S+, para energi'as debajo y arriba 
de la rodilla, respectivamente. En concordancia con las discusiones precedentes, en las 
figuras se observa que a bajas energias la difusion de rayos cosmicos se debe principalmente 
a fluctuaciones turbulentas de pequena escala (gobernadas por D^), mientras que a mayores 
energi'as el bulk de los rayos cosmicos se desplaza de acuerdo a los drifts que resultan del 
campo magnetico regular de gran escala, dando lugar a evidentes asimetrfas de gran escala 
en la distribucion de densidad. 

La infiuencia de los efectos de drift en el espectro de rayos cosmicos dependera, en general, 
del modelo de campo galactico, del conjunto particular de parametros especificados para el 
modelo (es decir, las amplitudes y escalas verticales asumidas para las componentes regular 
y estocastica del campo magnetico) y de la distribucion de las fuentes. Debe notarse que un 
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E/Z [eV] 

Figura 2.4: Espectro diferencial normalizado versus E/Z para el caso R — S j diferentes 
distribuciones radiales de fuentes: homogeneamente extendida (g(r) =cte.) y localizada 
(g(r) oc 6{r — r^), con = 3 y 6 kpc). 

parametro relevante en los coeficientes de difusion es el radio de Larmor de las particulas; 
esto implica que el impacto de los efectos del drift dependeran de la rigidez de los rayos 
cosmicos, es decir, del valor de E/Z. 

La fig ]2.3l muestra el flujo de rayos cosmicos normalizado como funcion de E/Z para todos 
los modelos de campo magnetico galactico descriptos anteriormente (ver la Seccion 2.1.1), 
con la fuente localizada en = 6 kpc, y con B* = 2 /iG, Zr = 5 kpc para el campo magnetico 
irregular. Puede notarse que en esta y las proximas figuras, el flujo de rayos cosmicos aparece 
multiplicado por E^^~^^^^, donde Pz se refiere al mdice espectral de la componente de carga 
Z en la fuente. Como la difusion transversal depende de la rigidez a traves de una potencia 
de 1/3, ec. ()2.26|) . resulta dNz/dE oc EJ-t^+^Z^ a bajas energfas, debajo de la rodilla. De 
la figura se ve que la mayon'a de los modelos considerados muestran la tendencia apropiada 
que da cuenta de la rodilla, ya que el cambio de mdice espectral aparece, por lo general, 
en el valor apropiado de E/Z (es decir, E ^ E^. para protones, y a mayores energfas para 
nucleos mas pesados). Los efectos de drift son particularmente notables para el caso del 
halo antisimetrico, como podria anticiparse de la ec. ()2.14|) teniendo en cuenta que, para 

~ ±signo(2;), aparece una velocidad de drift singular Ur — =f2-Da<^(-2) sobre el piano 
galactico. Siendo que alli es donde estan localizadas las fuentes, es razonable esperar que 
esta contribucion adicional a los drifts afecte los resultados significativamente. Por otra 
parte, en los modelos antisimetricos el halo contribuye con una velocidad de drift vertical Uz 
que converge o diverge desde el piano galactico (dependiendo de la orientacion del campo 
del halo relativa a la del disco), a diferencia de lo que ocurre en los modelos simetricos, en 
los que el drift vertical tiene la misma orientacion a ambos lados del piano. En modelos 
antisimetricos, se obtiene una rodilla mas pronunciada para ■ B^ > 0, mientras que, por el 
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Figura 2.5: Espectro diferencial normalizado versus E/Z para el caso R — S j una fuente 
localizada en = 6 kpc. Los parametros asociados al campo irregular toman los valores 
B* = 2 6 A /xG, y = 3 6 5 kpc, tal como se indica. 

contrario, una protuberancia {hump) justo antes de la rodilla se observa para Bd ■ Bh < 0. El 
pronunciado bump que muestra la fig l2.3l para este caso esta claramente en desacuerdo con los 
espectros observados (comparar con los datos experimentales de la figEZI^a)), de modo que, 
en principio, este modelo particular estaria desfavorecido por las observaciones. Sin embargo, 
Kempa et al. [71] hicieron notar que los datos experimentales podrian sugerir la existencia de 
un bump menos pronunciado justo antes de la rodilla. Desde esta observacion, realizada hace 
tres decadas, la existencia del bump sigue dando lugar a controversias; de cualquier forma, 
es interesante comentar que caracteristicas de este tipo podrian originarse en efectos de drift 
que aparecen en algunas configuraciones particulares de los campos magneticos regulares en 
la galaxia. 

La fig |2. 41 muestra el caso R — S para una distribucion de fuentes extendida g(r) = cte., 
y para fuentes localizadas en = 3 kpc y = 6 kpc, con el mismo campo irregular que se 
uso en la fig l2.2l (es decir, B* = 2 fiG j Zr = 3 kpc). Como es de esperar, los efectos de drift 
son mas intensos cuanto mayor es la distancia a la que esta situada la fuente, debido a que 
los drifts pueden remover del piano galactico a los rayos cosmicos que se propagan a lo largo 
de todo el camino entre la fuente y el observador. El mismo comportamiento cualitativo ha 
sido tambien observado en los otros casos estudiados. 

La fig ]2.5l exhibe el efecto de variar los parametros que describen al campo random. Los 
datos corresponden al caso i? — S* y los parametros asociados al campo irregular toman los 
valores B* = 2 6 4 /iG, j Zr = 3 6 5 kpc. Como se espera, los efectos del drift son mas 
intensos cuanto mayor es la supresion del campo random, ya sea asumiendo una amplitud 
menor, o bien disminuyendo su escala vertical caracterfstica. Un comportamiento analogo 
se observa en las otras configuraciones de campo relevantes, debido, en todos los casos, a la 
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Tabla 2.1: Abundancias fraccionales (para E = 1 TeV) e indices espectrales a bajas energfas 
(debajo de la rodilla) correspondientes a las componentes nucleares de rayos cosmicos desde 
liidrogeno liasta mquel (de [9,72]). 



mayor supresion de Da bajo condiciones de mayor turbulencia. 

Como estamos tratando con un escenario dependiente de la rigidez, los resultados obtenidos 
dependen de E/Z; para comparar nuestros resultados con los datos experimentales, debe 
calcularse la contribucion de cada componente de carga Z convolucionando el espectro pro- 
ducido por las fuentes con los efectos de modulacion que resultan de la propagacion en la 
galaxia. Finalmente, sumando sobre todas las componentes nucleares desde liidrogeno liasta 
niquel ^ (es decir, para 1 < Z < 28), obtenemos la prediccion del modelo para el espec- 
tro total de rayos cosmicos, que luego puede ser comparado con los espectros medidos por 
diferentes experimentos (KASCADE, CASABLANCA, DICE, PROTON y Akeno) [73-77]. 

A bajas energfas (debajo de la rodilla) el espectro observado localmente para cada com- 
ponente de carga Z viene dado por 



fz^J^] , (2.36) 



-az 

dE """"V^o/ 

donde (po es el flujo total de rayos cosmicos correspondiente a la energfa de referenda Eq, 
y fz es la abundancia fraccional de la componente nuclear de carga Z a la misma energia 
(Ylizfz = !)• Aquf adoptaremos Eq = 1 TeV, las abundancias e indices espectrales que 
resultan de los datos experimentales compilados en [9,72] (reproducidos en la Tabla EHj), 
y la normalizacion 0o = 3.5 x 10~^^m~^s~^sr~^eV~^. Asumiremos que los espectros pro- 
ducidos por las fuentes tienen indices espectrales {[3z) constantes (es decir, independientes 
de la energia), consistentes con los indices observados {az) a bajas energfas. Ademas, en 
la extrapolacion de los fiujos no tomamos en cuenta posibles energfas de corte asociadas al 
mecanismo de aceleracion en las fuentes. 



^La inclusion de elementos mas pesados es generalmente considerada irrelevante en este contexto, dado que 
son elementos mucho menos abundantes. De todos modos, en [72] se sugiere que estos micleos 'superpesados' 
(con Z > 28) pueden ser dominantes entre la rodilla y la segunda rodilla. En la Seccion 2.3 consideraremos en 
detalle el espectro de rayos cosmicos en esa region de energias, y discutiremos el posible rol de la componente 
superpesada. 
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Figura 2.6: Principales contribuciones al flujo total de rayos cosmicos, provenientes de pro- 
tones y nucleos de helio, oxigeno y hierro, obtenidas en el modelo R — S para una fuente 
localizada en = 6 kpc. 

La fig ]2.6l muestra las principales contribuciones al flujo total de rayos cosmicos, prove- 
nientes de protones y nucleos de helio, oxigeno y hierro, segiin se obtiene en el modelo 
R — S para una fuente localizada en = 6 kpc. Los parametros para el campo random son 
B; = 2 fiGj Zr = 3 kpc. 

La fig ]2.7f a) muestra el flujo total correspondiente a este caso en comparacion con algunas 
observaciones; como puede verse claramente, el acuerdo entre las predicciones y los datos 
experimentales es muy satisfactorio. 

La fig l2.7r b) muestra la dependencia de la composicion masica media (In A) con la 
energi'a, correspondiente al mismo caso de la figura anterior, comparada con las mediciones 
de diferentes experimentos [73-78] . Siendo que el escape de partfculas hacia el exterior de la 
galaxia depende de E/Z j que, por lo tanto, es mas efectivo para las partfculas mas livianas, 
una prediccion natural de este escenario es que la composicion se torne mas pesada arriba de 
la rodilla. Esta tendencia hacia componentes mas pesadas es la que manifiestan los datos de 
algunos experimentos, como por ejemplo las mas recientes observaciones de KASCADE [77] 
y EAS-TOP/MACRO [78] 

Aqui, sin embargo, debemos hacer notar que los resultados reportados por experimentos 
que utilizan diferentes tecnicas de deteccion no son compatibles en todos los casos, eviden- 
ciando discrepancias incluso en el comportamiento cualitativo de las observaciones. Si bien 
los datos mas recientes de varios experimentos parecen confirmar la tendencia hacia el pre- 

^°Es particularmente notable el acuerdo con los ultimos resultados de EAS-TOP/MACRO, dados a conocer 
luego de la publicacion de nuestro trabajo: "The present data explain therefore the observed knee in the 
cosmic ray primary spectrum as due to the steepening of the spectrum of a light component (p, He) at 
£^0 « 4 X 10'^^ eV, of A7 = 0.7 ± 0.4" [78]. 




37 




Figura 2.7: Comparacion entre las observaciones experiment ales y las predicciones teoricas 
obtenidas en el modelo R — S para una fuente localizada en = 6 kpc: (a) espectro total 
de rayos cosmicos; (b) composicion masica media (In A). 

dominio de las componentes mas pesadas por encima de la rodilla, deben tenerse en cuenta 
tambien otros resultados recientes [73] que sostienen la tendencia opuesta y que obligan a la 
cautela. En la fig J2.7r b) se puede ver que, segun algunos conjuntos de datos, la composicion 
masica media se liace mas liviana en ~ 10^^ eV, y comienza a crecer hacia componentes 
mas pesadas solo a partir de ~ 3 x 10^^ eV. Una depresion (dip) de esta naturaleza podria 
result ar de un bump en el espectro de la componente de protones en ese mismo rango de 
energias. 

De las discusiones y figuras precedentes, es claro que puede obtenerse una variedad de 
diferentes espectros modificando el valor de los parametros que intervienen en el modelo, o 
bien considerando otros modelos de campo galactico. Sin embargo, un analisis cuantitativo 
detallado de los diferentes resultados posibles no tiene mayor sentido en el marco de este tra- 
bajo. Por un lado, los datos experimentales muestran aiin errores estadisticos y sistematicos 
considerables, mientras que, por otra parte, las predicciones teoricas tambien se apoyan en 
un cierto mimero de consideraciones simplificadoras. Por ejemplo, la ecuacion de transporte 
fue considerada como puramente difusiva, despreciando la fragmentacion de los micleos en 
el medio interestelar, la reaceleracion, o la conveccion debido a la posible existencia de un 
viento galactico. Ademas, la dependencia espacial de las diferentes componentes del campo 
magnetico galactico fue tratada en terminos muy simples; aunque los modelos que hemos 
adoptado son plausibles y capturan las caracteristicas esenciales de los campos magneticos 
en la galaxia, una descripcion mas realista exigin'a abandonar las simetrfas del sistema, in- 
corporar parametros adicionales, etc. La tendencia de los resultados obtenidos, sin embargo, 
muestra un acuerdo notable con las observaciones experimentales. 
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2.3 El espectro y las anisotropias de rayos cosmicos 
desde la rodilla hasta la segunda rodilla 

En la Seccion precedente, hemos visto que el modelo de difusion turbulenta y drift provee una 
explicacion natural para la rodilla del espectro de rayos cosmicos, basada en propiedades bien 
establecidas sobre la propagacion de parti'culas cargadas en campos magneticos regulares y 
turbulentos. El espectro de una componente de carga Z comienza a verse afectado por los 
drifts a partir de E c:^ Z Er, con el mdice espectral cambiando finalmente en Aa ~ 2/3 en 
alrededor de una decada de energi'a. De este modo, la rodilla en el espectro total (es decir, 
el cambio de mdice 2.7 — 3 a partir de £^r. ~ 3 x 10^^ eV) aparece debido a la supresion de 
las componentes livianas, fundamentalmente protones y niicleos de helio. En esta Seccion, 
veremos que este escenario permite explicar tambien la segunda rodilla del espectro total (el 
cambio de mdice 3 3.3 a partir de Esr — 4 x 10^'' eV) como resultado de la supresion 
de las componentes pesadas (correspondientes al grupo del liierro); incorporando, como es 
usual, una componente de rayos cosmicos extragalacticos que domina a partir del tobillo 
{Et ^ 5 X 10^*^ eV), mostraremos que este escenario es capaz de reproducir el espectro 
observado hasta la region de ultra-alta energi'a. 

Con respecto a la composicion, en la region de la segunda rodilla no se dispone de medidas 
precisas (la magnitud de los errores en las medidas por encima de 10^^ eV, que pueden verse en 
la fig l2.7f b). da una buena nocion acerca de las limitaciones experiment ales en este aspecto). 
Sin embargo, existen algunas observaciones de interes que permiten, hasta cierto punto, un 
contraste entre la teoria y el experimento; en particular, la determinacion de Haverah Park 
de las abundancias nucleares usando la distribucion lateral de Uuvias atmosfericas [79] (que 
sugiere que, en el rango 2 x 10^ ''-10^^ eV, solo un ~ 30% de los rayos cosmicos es liviano), y 
la observacion de HIRES-MIA de un cambio de composicion hacia una mezcla mas liviana 
a partir de la segunda rodilla [80]. 

En este contexto, las observaciones referidas a la anisotropia del flujo de rayos cosmicos 
representan una herramienta adicional de gran importancia. En particular, la colaboracion 
AGASA [81, 82] reporto una anisotropia del ~ 4% para E ~ 0.8-2 x 10^^ eV, con un 
exceso observado en una direccion cercana al centro galactico y un deficit cercano a la 
direccion del anticentro galactico. Adicionalmente, ha sido observado algun exceso de eventos 
en la direccion de la region de Cygnus (a lo largo del brazo espiral de Orion). En los 
datos de SUGAR se ha encontrado tambien un exceso desde una direccion cercana al centro 
galactico [83], y con una direccion separada apenas unos grados del exceso reportado por 
AGASA, aunque de una extension espacial menor (consistente con una fuente puntual). 
Debe notarse que estas observaciones ocurren justamente en la region de energi'as en la que la 
componente galactica de rayos cosmicos comienza a atenuarse para dar lugar al predominio de 
la componente extragalactica. Esto sugiere que la explicacion de estas observaciones podn'a 
estar directamente vinculada con el mismo proceso que es responsable de la atenuacion de 
la componente galactica, y veremos que, en efecto, el escenario de difusion turbulenta y drift 
las puede reproducir adecuadamente. 
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Como un ultimo examen de este escenario, estudiaremos (para diferentes modelos de 
campo magnetico) la dependencia con la energi'a de las predicciones sobre el flujo, la amplitud 
de anisotropia y la fase del primer armonico en ascencion recta, comparando los resultados 
con las observaciones experimentales en todo el rango que se extiende desde la rodilla hasta 
el tobillo del espectro. Encontraremos que las predicciones muestran la tendencia adecuada 
para dar cuenta de las observaciones, dando robustez al modelo de difusion turbulenta y 
drift como posible explicacion de los fenomenos observados en los rayos cosmicos galacticos 
de muy alta energia. 

Siguiendo los pasos de la Seccion anterior, consideraremos que la componente galactica de 
rayos cosmicos se obtiene de resolver la ecuacion de difusion correspondiente, luego de asumir 
un dado modelo de campo magnetico regular, una realizacion particular de los parametros 
que describen los campos, y una distribucion radial para las fuentes. Adicionalmente, aliora 
incorporaremos una componente extragalactica dada por 

La amplitud de esta componente se ajusto para reproducir la normalizacion correcta de los 
flujos observados para E ~ 10^^ eV. Si consideramos que el flujo extragalactico es isotropico 
en el borde de la galaxia, entonces, en ausencia de efectos de reaceleracion, tambien sera 
isotropico en el interior (es decir, la densidad de rayos cosmicos no puede ser incrementada 
localmente, dentro de la galaxia, solo por procesos difusivos ^^). En consecuencia, la depen- 
dencia con la energi'a de la componente extragalactica debe ser similar a la que se infiere 
para la produccion en las fuentes (es decir, dada por un mdice espectral /3 ~ 2.4). 

La fig ]2.8l muestra los espectros que resultan de considerar una componente galactica (cal- 
culada, como en la Seccion anterior, para diferentes modelos de campo magnetico y diferentes 
distribuciones de fuentes) y la extragalactica, ec. (j2.37p La eleccion de los parametros 
para cada modelo estuvo guiada solo por el requerimiento de reproducir correctamente la 
primera rodilla, dentro del rango de valores plausibles. Tal como muestra la figura, el acuerdo 
con los datos experimentales es, en general, muy bueno. 

Las abundancias relativas de las diferentes componentes, agrupadas en livianas {1 < Z < 
5), intermedias {Q < Z < 19) y pesadas (20 < Z < 28), se muestran en la fig l2.9l como funcion 
de la energia, asumiendo que el flujo extragalactico esta constituido por protones. De alli 
resulta que en nuestro modelo la composicion de rayos cosmicos para ~ 5 x 10^'' eV consiste 

"'^"'^Esto puede entenderse como una consecuencia del teorema de Liouville [84,85], o bien puede verse 
directamente de la ecuacion de difusion, notando que la solucion en ausencia de fuentes y con condiciones de 
borde isotropicas es la de una densidad constante, de modo que los efectos difusivos en la galaxia no pueden 
producir un incremento local de densidad en la componente extragalactica. 

^^Con el fin de preservar la presentacion de las figuras en su contexto original (es decir, segiin la version 
en la que han aparecido publicadas), se opto por sacrificar la uniformidad de presentacion a lo largo de esta 
Tesis. Asf, por ejemplo, los espectros de la figl^Tfa) se presentan en terminos de (j>E'^-^ [MeV"'^'^cm~^s~"'^sr~^], 
y en la fig l2.8l en cambio, como <j)E^ [eV^m^^s^^sr^^]. 
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Figura 2.8: Espectros de rayos cosmicos obtenidos para diferentes modelos de campo 
magnetico y distribuciones de fuentes, en comparacion con los datos experimentales. La 
recta punteada representa el flujo extragalactico, ec. ()2.37p . 

en un ~ 50% de nucleos del grupo del hierro (fundamentalmente Fe y Mn), en alrededor de 
un ~ 20% de nucleos de masa intermedia, mientras que el resto corresponde esencialmente 
a nucleos galacticos de helio y a protones extragalacticos. Estos valores son consistentes con 
la determinacion de Haverah Park de las abundancias nucleares [79], que observa solo un 
~ 30% de componentes livianas en el rango 2 x 10^^-10^^ eV. Por otra parte, la observacion 
de HIRES-MIA de un cambio de composicion liacia una mezcla mas liviana a partir de la 
segunda rodilla [80] se explica en este escenario como debida a la supresion del grupo de 
nucleos galacticos pesados, que da lugar al predominio de los protones extragalacticos hacia 
la region del tobillo. 

Consideremos ahora las anisotropfas que resultan en este escenario. Como ya hemos 
comentado, estas seran producidas solo por la componente galactica, ya que, si suponemos 
que la componente extragalactica es originalmente isotropica, conservara su isotropfa luego 
de que los efectos de propagacion en la galaxia sean tenidos en cuenta. La componente 
galactica de rayos cosmicos de carga Zi contribuye a la anisotropia total una cantidad dada 
por [3] 

^ • 

donde Ji es la corriente macroscopica correspondiente a esa componente, dada por una 
expresion analoga a la ec. ()2.12j) . pero involucrando a la densidad Ni que corresponde a los 
rayos cosmicos de carga Zj. La anisotropia total viene dada por 6 = J2i fi^ii donde /* = 
Ni/N son las abundancias relativas de todas las especies nucleares, y donde N es la densidad 
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Figura 2.9: Abundancias relativas de las diferentes componentes de rayos cosmicos, asumien- 
do que el fiujo extragalactico esta constitmdo por protones. Los elementos estan agrupados 
en componentes livianas (1 < 2' < 5), intermedias {Q < Z < 19) y pesadas (20 < Z < 28). 
Los datos corresponden al mismo modelo de campo magnetico y a la misma distribucion de 
fuentes que en la fig l2.10l 

de rayos cosmicos total, obtenida de sumar las contribuciones galactica y extragalactica 
La presencia de la componente extragalactica isotropica tiende a reducir las fracciones /* 
de las componentes galacticas, con el efecto de suprimir el crecimiento de la amplitud de la 
anisotropfa total hacia la region del tobillo. 

Como hemos mencionado antes, a bajas energias la difusion esta dominada por D±] de 
la ec. ()2.12|) se ve que, en ese caso, la contribucion a la anisotropfa de los rayos cosmicos 
de carga Zi es Si oc D_iVNi/Ni, de forma que 5i oc E^^^. Sin embargo, a mayores energi'as 
hay un crossover hacia el regimen de dominio del drift, conduciendo al comportamiento 
6i oc E. Aiin a mayores energias, a partir de ~ 10^^ eV, la densidad de rayos cosmicos 
comienza gradualmente a quedar dominada por la anisotropfa intrfnseca de la componente 
extragalactica (que, por simplicidad, fue considerada despreciable en este trabajo). 

De la ec. ()2.12j) se ve que la contribucion a la corriente macroscopica que resulta de 
la difusion transversal es en la direccion de VA^ (y ortogonal a las curvas de nivel de la 
densidad), mientras que la contribucion del drift es normal a VA^ (y, por lo tanto, paralela 
a las curvas de nivel). 

En lo que sigue, estudiaremos la anisotropfa asociada a una linica componente galactica^^; 
mas adelante, analizaremos la anisotropfa total que resulta de considerar el conjunto de 

""^•^Notese aqui la diferencia con las abundancias galacticas fz que se introdujeron en la Seccion anterior 
(ec. (jOHl V que no tienen en cuenta la componente extragalactica. 

^"^Para simplificar la notacion, omitiremos el subindice que distingue entre las diferentes contribuciones 
galacticas. 
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Figura 2.10: Curvas de nivel de densidad para rayos cosmicos de diferentes energi'as, corre- 
spondientes a: a) E/Z = 10^^ eV, b) E/Z = E^, y c) E/Z = 1.5 x 10^^ eV/26. Cada dos 
curvas de nivel, la densidad cambia en un orden de magnitud. Las fleclias representan la 
direccion de los vectores de anisotropfa, mientras que los sombreados indican la amplitud 
de las anisotropfas. El modelo corresponde a una fuente constante contenida en un disco de 
radio 3 kpc en el piano galactico, y la estructura de campo magnetico descripta en el texto. 



todas las componentes galacticas y la extragalactica. Como la corriente macroscopica es 
perpendicular al campo magnetico regular, esta contenida en el piano r — z; en consecuencia, 
el vector de anisotropfa tambien esta contenido en este piano, con sus componentes dadas 
por 

3 f ^ ON , ^ dN\ 



y 

dN , ^ dN\ 



Las figs ]2.10l y 12.111 muestran, para diferentes modelos de campo magnetico y distintas 
distribuciones de fuentes, las curvas de nivel asociadas a la densidad de rayos cosmicos. 
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asi como tambien la direccion de los vectores de anisotropia y sus amplitudes (indicadas 
por medio de los sombreados). Los resultados dependen del valor de EjZ^ es decir que 
corresponden a diferentes energfas de acuerdo a la especie nuclear particular que se considere. 
La importancia relativa de la difusion transversal y el drift para las distintas regiones de la 
galaxia puede ser inferida directamente a partir de la direccion entre las curvas de nivel y 
las flechas. 

La fig l2.10l resulta de considerar, para diferentes energfas, un modelo NR — 5*+ con 
parametros {Bl^i'^° = -0.75,5*'^"'° = -0.75,5; = 1.5, = 0.5, z/, = 5,Zr = 5) (los 
campos medidos en /xG y las escalas verticales en kpc), y con una fuente constante contenida 
en un disco de radio 3 kpc en el piano galactico. La fig l2.10f a) corresponde a energfas por 
debajo de la primera rodilla {E/Z = 10^^ eV), donde la difusion esta casi completamente 
dominada por D±. La fuente corresponde a un disco uniforme de radio r = 3 kpc contenido en 
el piano galactico. Puede verse que las curvas de nivel resultantes son simetricas con respecto 
al piano galactico (comparar con la fig ]2.2r a)). mientras que las anisotropfas son de amplitud 
pequefia (~ 10^^) y perpendiculares a las curvas de nivel. La fig ]2.10f b) corresponde a 
E/Z = Er, donde los efectos de drift comienzan a ser mas notables y la anisotropia local es del 
orden de ~ 10~^. Finalmente, la fig l2.10f c) corresponde a la rigidez E/Z = 1.5 x 10^® eV/26 
(es decir, muestra el comportamiento que corresponde, por ejemplo, a micleos de hierro para 
ii^ = 1.5 X 10^^ eV). En este regimen, las densidades estan en gran medida determinadas por 
los drifts, de forma que reflejan las asimetrfas del campo magnetico regular a gran escala 
(comparar con la fig l2.2f b)). La amplitud de anisotropia local es, en este caso, del orden de 
10-2. 

En la fig |2.11l se considera el regimen dominado por los drifts (con E/Z = 1.5 x 
10^^ eV/26, la rigidez correspondiente a la fig l2.10r c)). para fuentes localizadas en diferentes 
posiciones radiales r^, y para diferentes parametros y estructuras de campo magnetico. Si 
bien todas estas realizaciones ajustan los datos del espectro local de rayos cosmicos adecuada- 
mente (ver la fig ]2.8|) . se observa que las densidades globales, las corrientes macroscopicas y 
las anisotropfas varfan significativamente, por lo general, de un modelo a otro. En conse- 
cuencia, las anisotropfas locales observadas pueden proveer una herramienta muy litil para 
establecer las propiedades generales del campo magnetico. 

Resulta importante seiialar que, tanto el exceso de eventos hacia el centro galactico como 
el deficit hacia el anticentro galactico observados por AGASA [81], podrfan ser reproducidos 
mediante un flujo difusivo de rayos cosmicos que este localmente dirigido hacia el anticentro 
galactico, con una amplitud que de lugar a la anisotropia observada del ~ 4%. En el rango 
de energfas relevante a estas observaciones, la difusion de Hall domina sobre la difusion 
transversal; localmente, = —1, de modo que 

^^^^'^ (2.41) 



c N dz c N dr 

Entonces, vemos que la anisotropia local estara dirigida aproximadamente hacia la direccion 



iSDefiniendo B^'^^'^" = B^''^''{rob.. Zobs) Y B^''^'" = B^'''°{robs, Zobs), donde {vobs^Zobs) = (8.5 kpc,0). 
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r=0 NR-S (-1,-1,2,0.5,5,4) 

s 



r^=0 NR-A+ (-1,-1,3,0.5,2,3) 





r^=6 NR-S (-1,-1,1,1,4,4) 




Figura 2.11: Densidades, corrientes y anisotropias para diferentes modelos de campo 
magnetico y distintas distribuciones de fuentes, en el regimen dominado por los drifts. Entre 
parentesis se indican los valores de {B^'^'"^'^", Bq^"''", B*, Zd, z^, Zr) en /iG y kpc, respectiva- 
mente. 

radial si, localmente, \dN/dz\ ^ \dN/dr\. Por otra parte, volviendo sobre la ec. ()2.20|) . 
que describe al coeficiente de difusion antisimetrico, podemos estimar la amplitud de la 
anisotropia: 

15 I ~ ^ 04 ^ 26kpc^ [utaY 

' ^' - hl + (cur^)2 - ^-^^ low Z h |Bo| 1 + (cuTAy ' ^ ^ ^ 

donde h = \d\nN/dz\~^ es la escala vertical caracteristica de la variacion local de la densidad 
de rayos cosmicos. Al considerar la mezcla de todas las especies nucleares, debe evaluarse 
S{E) = fi{E)^i{E)'i como hemos comentado, el drift tiende a remover primero las com- 
ponentes mas livianas, pero a mayores energi'as (hacia la region del tobillo) la presencia de 
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la componente de protones extragalactica tiende a reducir las abundancias relativas de todas 
las especies nucleares galacticas (ver la fig l2.9|) . Si bien en la componente galactica dominan 
los niicleos del grupo del hierro, debe notarse que la anisotropia, para una dada energia, es 
mayor para las componentes mas livianas (de la ec. ()2.42p . se ve que 5z ~ 1/^)- Podemos 
escribir 



donde la suma es sobre las componentes galacticas. Considerando, como valores indicativos 
para 1.5 x 10^^ eV, abundancias del orden del ~ 1% para protones (galacticos), ~ 3% para 
niicleos de lielio, ~ 5% para micleos de oxfgeno y ~ 20% para nucleos de hierro, se obtiene 
que 5 ~ dpe- De este modo, la anisotropia de los nucleos de hierro a esta energia resulta un 
buen estimador de la anisotropia total, y sera lo que discutiremos en lo que sigue (antes de 
tratar con el caso detallado en el que sumamos la contribucion de todas las componentes). 
En particular, las figs l2.1Uf c) v l2.11l muestran dpe^ esta energia para diferentes realizaciones 
particulares de campo magnetico y fuentes que reproducen localmente las observaciones de 
anisotropia de AG AS A que ya hemos mencionado. 

De la ec. ()2.42|) . se ve que las anisotropias crecen linealmente con la energia y pueden 
alcanzar los valores observados, del orden de 5 ~ 0.04, si la escala vertical caracteristica para 
las variaciones de la densidad de rayos cosmicos es del orden del kpc. En los modelos de 
campos magneticos simetricos aparecen, de un modo natural, gradientes verticales significa- 
tivos en la densidad de rayos cosmicos sobre el piano galactico, ya que en esos modelos los 
drifts verticales estan igualmente orientados en ambos hemisferios. 

Para que el vector de anisotropia apunte en la direccion radial hacia afuera, el valor 
local de dN/dz debe ser negativo (ver la ec. (j2.41|) ). Para esto, las corrientes macroscopicas 
verticales deben estar esencialmente orientadas en la direccion negativa, algo que es mas 
pronunciado, por ejemplo, en modelos donde un halo simetrico esta dirigido en la direccion de 
—(f). Por otra parte, para tener localmente \dN/dr\ < \dN/dz\, es conveniente situar la fuente 
en una posicion no muy cercana a nuestra ubicacion. El modelo NR—A+ ffig r2.111 arriba a la 
derecha) muestra, ademas, la interesante caracteristica de que los drifts tienden a converger 
hacia el piano galactico, resultando en un drift dirigido hacia afuera para 2; ~ 0. Por el 
contrario, un modelo de halo con inversiones radiales en la orientacion del campo (o bien, un 
modelo en el que se considera solamente el campo regular del disco galactico) exhibe drifts 
verticales muy pronunciados ffig ]2.11l abajo a la izquierda). El factor {u!Ta)'^/{1 + (ujta)'^) 
en las ecs. (I2.20j) y (|2.42j) refleja la supresion de en presencia de alta turbulencia. Si bien 
tiende a la unidad para Br <^ |Bo|, puede hacerse pequeno si la turbulencia es alta, con el 
efecto de suprimir las anisotropias. 

Otras propuestas, elaboradas con el fin de explicar el exceso de eventos en la direccion del 
centro galactico, asocian la anisotropia en esta direccion a rayos cosmicos que se propagan en 
forma casi rectilinea desde una fuente central [86-89]. Sin embargo, aiin en el caso de que esos 
rayos cosmicos fueran protones, la deflexion producida por los campos magneticos galacticos 
(regulares y random) excederia largamente el desplazamiento de 10° que se observa entre 




(2.43) 



z 
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la direccion del exceso de eventos y la del centro galactico; para producir una deflexion tan 
pequena, la fuente deben'a estar a una distancia no mayor de ~ 2 kpc [87]. La liipotesis que 
atribuye el exceso a neutrones producidos a partir de protones o nucleos acelerados en fuentes 
en el centro galactico (mediante la produccion de fotopiones o por procesos de fragmentacion 

ha sido explorada como un escenario alternativo [81,87,89]. La longitud caracteristica 
de decaimiento de los neutrones es ~ 10 kpc eV), de modo que, para E ~ 10^^ eV, 

es justamente del orden de la distancia galactocentrica. Ademas, esta propuesta tambien 
explicarfa porque no se observa ningiin exceso significativo a menores energfas. Sin embargo, 
disponer de una fuente de neutrones tan luminosa no es particularmente natural; en cambio, 
la energfa maxima de los protones galacticos en el escenario de la difusion turbulenta y 
drift no necesita exceder ~ 10^'' eV, un requerimiento muclio mas plausible. Por otra parte, 
tampoco resulta evidente, en el escenario que propone el decaimiento de neutrones, como 
puede explicarse el deficit observado en la direccion del anticentro. En lo que respecta al 
exceso de eventos observados en la direccion de Cygnus, aqui solo liaremos notar que, en 
principio, una explicacion similar a la que hemos venido desarroUando en esta Seccion puede 
encontrarse en terminos de la difusion de rayos cosmicos a lo largo del brazo espiral. Como ya 
hemos comentado, esto requerirfa remover la condicion de simetrfa cilmdrica para permitir 
la difusion paralela de rayos cosmicos en la direccion azimutal, gobernada por D\\ (ver las 
ecs. ()2.1H) y ()2.12p ). En cualquier caso, una prediccion general que resulta en este escenario 
es que deberia observarse un deficit en la direccion opuesta, algo que resulta de potencial 
interes para las futuras observaciones en el hemisferio sur con el observatorio AUGER [33]. 

En la discusion precedente, hemos analizado en detalle la relacion entre el vector de 
anisotropfa y las corrientes macroscopicas de difusion transversal y de drift; ademas, estu- 
diamos la dependencia de las anisotropias con diferentes estructuras de campo regular y, 
finalmente, tomando la anisotropfa del hierro como estimativa de la anisotropfa total a en- 
ergfas del orden de ~ 10^^ eV, mostramos que el escenario de difusion turbulenta es capaz de 
reproducir las observaciones de anisotropfa de AGASA. Desde un enfoque mas cuantitativo, 
en lo que sigue consideraremos, para diferentes realizaciones de campo magnetico y distin- 
tas distribuciones de fuentes, la dependencia con la energfa de las medidas de amplitud de 
anisotropfa y fase del primer armonico en ascencion recta, para todo el rango de energfas que 
comprende las dos rodillas (es decir, aproximadamente desde ~ 10^^ eV hasta ~ 10^^ eV). 

En la fig J2.12l los espectros y las amplitudes de anisotropfa obtenidas para los modelos 
de campo magnetico NR — A+ y NR — S+ se comparan con las observaciones de diferentes 
experimentos [16,43,90]. Debe notarse que, a partir de aquf, estamos tratando con las 

^^En el Capi'tulo 4 se estudiara en detalle este tipo de procesos, en el contexto de una explicacion alternativa 
para la rodilla del espectro de rayos cosmicos. 

^^Los puntos experimentales de anisotropfa se obtienen de la amplitud del primer armonico en ascencion 
recta A, que se relaciona con la anisotropfa S a traves de A = S cos d cos A, donde d es la declinacion a la cual 
las observaciones son realizadas, mientras que A es la latitud de la direccion donde el flujo es maximo. Los 
puntos que aparecen en el grafico corresponden a, A/ cos d para los diferentes experimentos, de forma que, en 
verdad, representan una cota inferior para S, dado que cos A no puede ser recuperado mediante un analisis 
armonico en ascencion recta. 
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NR-A+ (2,1,4,5) NR-A+ (2,1,4,5) 




E [cV] E [eV] 




Figura 2.12: Espectros y amplitudes de anisotropia calculados para diferentes modelos de 
campo magnetico galactico, con Bq'^^^'^° = Bq'^"-'-° = —1, y los parametros {B*,Zd,Zh,Zr) en 
/iG y kpc, respectivamente, indicados sobre cada figura. Las rectas punteadas corresponden 
al fiujo extragalactico dado por la ec. ()2.37|) . Tambien se muestran los datos experimentales 
relevantes. 

anisotropias totales, calculadas con la inclusion de todas las componentes nucleares galacticas 
y la extragalactica de protones. En esta figura, se consideran tanto fuentes localizadas (en 
= 3 y 6 kpc) como extendidas en un anillo de densidad uniforme (para 3 < < 6 kpc). 
Como ya hemos observado anteriormente, los efectos de drift son mas notables para fuentes 
mas lejanas, ya que las corrientes de drift pueden remover del piano galactico a los rayos 
cosmicos que se propagan a lo largo de todo el camino entre la fuente y el observador. Por 
otra parte, las distribuciones de fuentes extendidas tienden a producir espectros mas pianos 
(es decir, muestran un cambio de fndice espectral menos abrupto en la region de la rodilla), 
debido a que los gradientes de densidad tienden a ser menores. Finalmente, los espectros 
de la fig l2.12l tambien muestran que, cuando se considera un anillo liomogeneo extendido de 
fuentes, el flujo observado dominante es el producido cerca del borde mas cercano al punto 
de observacion. 
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De la fig J2.12l se observa que, asumiendo valores plausibles para los parametros de los 
modelos de campo, el escenario de difusion turbulenta y drift es capaz de reproducir no solo 
las observaciones del espectro de rayos cosmicos, sino tambien las medidas de anisotropia. 
De las ecs. fj2.12|) y ()2.38|) puede verse que |ViV| juega un papel crucial en el incremento 
de la amplitud de la anisotropia, de modo que fuentes cercanas tienden a producir mayores 
anisotropias que otras mas lejanas. El efecto de considerar distribuciones de fuentes extendi- 
das puede ser comprendido como el resultado de una suma de los resultados producidos por 
fuentes localizadas, pesada por la contribucion relativa de cada una a la densidad total de 
rayos cosmicos observada. De esta forma, puede suceder que solo una pequeiia region dentro 
de la distribucion extendida sea responsable de las contribuciones dominantes al espectro 
y la anisotropia observadas. O bien, puede ocurrir que diversas regiones de la distribucion 
de fuentes contribuyan con corrientes macroscopicas comparables en magnitud pero orien- 
tadas segiin direcciones muy diferentes, dando como resultado que la anisotropia aparezca 
suprimida en algunas regiones, e incrementada en otras. En cualquier caso, los resultados 
que aparecen en la fig l2. 121 para fuentes localizadas deberfan ser indicativos de la tendencia 
de los que se obtendn'an para otras distribuciones de fuentes. 

Aunque algunos resultados para el caso NR — S+ exhiben anisotropias algo bajas, debe 
notarse que, aun en el caso de fuentes localizadas, estamos en realidad tratando con dis- 
tribuciones anulares alrededor del centro galactico, ya que nuestro enfoque asume desde el 
comienzo la simetria azimutal del sistema. Removiendo esta suposicion simplificadora y 
resolviendo la ecuacion de difusion tridimensional para una fuente puntual cercana, podri'a 
obtenerse un incremento adicional de la anisotropia. En efecto, ya ha sido reconocido en la 
literatura [3, 91] que algunas fuentes cercanas discretas (por ejemplo, piilsares y remanentes 
de supernova tales como Vela, Loop III y Geminga) podri'an proveer la contribucion domi- 
nante a las anisotropias observadas a bajas energias. Una estimacion precisa de este efecto 
requiere, sin embargo, indicar la distribucion espacial, la luminosidad, la edad y la evolucion 
de las fuentes; en consecuencia, no parece necesario buscar un ajuste preciso de nuestros 
resultados a los datos experimentales debajo de la rodilla, aunque es interesante observar 
que las predicciones en este escenario tienden en general a exhibir un ajuste muy aceptable 
a las observaciones, especialmente en la region de altas energias. 

Al sumar sobre todas las componentes galacticas para obtener el espectro total de rayos 
cosmicos, hemos estado considerando la contribucion de todas las especies nucleares desde 
hidrogeno hasta niquel. De todos modos, siguiendo las sugerencias de [72], tambien explo- 
ramos la posible contribucion de los elementos 'superpesados' hasta el uranio (es decir, con 
28 < Z < 92). El mdice espectral az de un elemento con Z <28 puede medirse a partir del 
espectro observado a bajas energias (debajo de ~ 10^^ eV), y de el puede inferirse el indice Pz 
correspondiente a la produccion en las fuentes. En cambio, para los elementos superpesados 
no hay informacion observacional disponible acerca de los indices espectrales. En [72], se 
propuso obtenerlos mediante la extrapolacion de la posible dependencia con Z de los indices 
espectrales asociados a los elementos mas livianos, usando una parametrizacion dada por 
—az = A + BZ^ para ajustar los datos. Mientras que esta expresion (en la que intervienen 
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tres parametros a ajustar) favorece una extrapolacion no lineal (con C = 1.51 ± 0.13), tam- 
bien el caso de una extrapolacion lineal (C = 1) fue considerada, encontrandose un buen 
ajuste con los indices espectrales medidos [72]. Asumiendo la extrapolacion lineal de los 
indices espectrales, nosotros encontramos que la contribucion de las especies superpesadas 
es completamente despreciable en el escenario de la difusion turbulenta y drift. Estos elemen- 
tos superpesados solamente tendrian algiin efecto notable con las extremas extrapolaciones 
no lineales sugeridas en [72] (segiin las cuales el mdice espectral del uranio, por ejemplo, 
resultan'a au — 1.9), pero en nuestro modelo esto conduciria a un marcado exceso en el 
espectro arriba de ~ 10^^ eV, excepto que se considerara un cutoff dependiente de la rigidez 
en el espectro producido en las fuentes, para energfas cercanas a ^ Z x 10^'' eV. En 
cambio, en [72] se asume la hipotesis ad hoc de que el mdice espectral en la region de la 
rodilla cambia en Aa ~ 2, y la supresion del espectro resultante termina siendo claramente 
excesiva por encima de la segunda rodilla Por otra parte, seria natural que el espectro 
producido en las fuentes fuera universal (en terminos de la rigidez), de modo que las difer- 
encias observadas en los indices espectrales a bajas energias solo reflejasen los efectos del 
proceso de fragment acion que afecta a la propagacion de los micleos pesados. Entonces, no 
deben'a esperarse en este caso ningun incremento de micleos superpesados por encima de la 
rodilla. 

Las amplitudes de anisotropfa para diferentes componentes galacticas de rayos cosmicos 
se muestran en la fig l2.13l como funcion de la energi'a, para una fuente anular extendida 
{rg = 3 — 6 kpc), el caso NR—A+ y los mismos parametros que se usaron en la fig l2.12l Como 
ya hemos discutido antes, el regimen dominado por los drifts es un fenomeno dependiente de 
la rigidez que incrementa la anisotropfa de todas las componentes galacticas. Hemos visto 
tambien que, aunque los elementos mas livianos (particularmente, los protones y los micleos 
de helio) estan fuertemente suprimidos a altas energias debido al escape de la galaxia, aiin 
asi contribuyen significativamente a la anisotropia total. Por ejemplo, los micleos de helio 
contribuyen, para el caso de la fig J2.13l alrededor del 30 — 40% de la anisotropfa total de 
los rayos cosmicos galacticos en todo el rango de energfa entre 10^^ y 10^^ eV, aunque la 
abundancia cae debajo del 10% para 10^^ eV. 

La fig J2.14l muestra estos mismos resultados en comparacion con los Ifmites experimentales 
de KASCADE [92-94], obtenidos del reciente analisis de ~ 10^ eventos en el rango de energfa 
(0.7—6) X 10^^ eV. Si bien la cota de KASCADE es compatible con observaciones anteriores, se 
espera que este experimento (o bien su sucesor, KASCADE-Grande) incremente la estadfstica 
y pueda poner a prueba las predicciones de este escenario [92-94]. Por otra parte, serfa de 
gran interes que se pudieran resolver las anisotropfas de las diferentes componentes por 
separado, ya que eso permitirfa discriminar entre algunos de los modelos propuestos hasta el 

-"^^En el Capi'tulo 3 se discutira, en otro contexto, el ajuste al espectro de rayos cosmicos en escenarios 
dependientes de la rigidez. Teniendo en cuenta solo componentes pesadas hasta Z = 28, se observa que 
un cambio de fndice espectral de Aa ~ 2 conduce a una supresion excesiva en el espectro total a partir de 
~ 10^^ eV (ver la fief l3.1|l . Con la extrapolacion no lineal de los indices espectrales de elementos superpesados, 
se evidencia un problema similar por encima de la segunda rodilla [72]. 
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Figura 2.13: Amplitudes de anisotropia correspondientes a las principales componentes de 
rayos cosmicos galacticos: protones y niicleos de helio, oxigeno y hierro. Para comparar, 
tambien se muestra la amplitud de anisotropia total, que involucra a las componentes 
galactica y extragalactica. Los datos corresponden a un anillo homogeneo extendido de 
fuentes (3 < < 6 kpc) para el caso NR — A+, con los mismos parametros que en la 
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Figura 2.14: Amplitudes de anisotropia en el escenario de la difusion turbulenta y drift, en 
comparacion con los recientes limites experiment ales de KASCADE (de [92-94]). 



momento. KASCADE, en particular, podn'a ser capaz de realizar este tipo de observaciones 
en un futuro proximo [92,95]. 



51 




Figura 2.15: Fase de anisotropi'a del primer armonico en ascencion recta como funcion de la 
energi'a, correspondiente a los mismos casos de la fig ]2.12l Los datos medidos por diferentes 
experimentos tambien se muestran para comparar. A la derecha, se indican direcciones 
de referenda en la galaxia (ACG: Anticentre galactico, PNG: Polo norte galactico, etc.), 
tomando en cuenta que S esta contenida en el piano r — z. 

En la fig J2.15l se grafica la fase de la anisotropfa del primer armonico en ascencion recta 
en funcion de la energia, comparada con las observaciones compiladas en [43,90], para los 
mismos casos que se trataron en la fig l2.12[ A bajas energfas, la difusion transversal domina 
y la direccion de maxima intensidad corresponde al centro galactico, como ya hemos visto 
antes (ver la fig l2.10f a)). A energfas mayores, la fase de la anisotropfa pasa a depender 
de la geometrfa detallada del campo regular adoptado; en los modelos considerados en esta 
figura, las diferencias no son grandes (y tienden a indicar un exceso mas bien hacia el polo 
norte galactico), pero, tal como hemos visto antes para otros modelos de campo magnetico, 
el maximo de la anisotropfa puede estar en la direccion del centro galactico para energfas 
cercanas a 10^^ eV (ver figs ]CTir c)- ITTT|l . 

2.4 Una reconsideracion de los coeficientes de difusion 

Como hemos visto en las secciones precedentes, la apropiada inclusion de los efectos de drift 
en la propagacion de rayos cosmicos en los campos magneticos turbulentos de la galaxia 
es un elemento crucial que permite explicar el espectro, la composicion y las anisotropfas 
observadas en la region que comprende las dos rodillas. La difusion turbulenta y el drift 
son caracterfsticos del transporte de rayos cosmicos tambien en otros escenarios de interes, 
como por ejemplo en la propagacion de rayos cosmicos en los medios interplanetario o in- 
tergalactico, o en el mecanismo de Fermi de aceleracion en diversos objetos astroffsicos. 



52 



Sin embargo, la informacion disponible sobre el comportamiento del tensor de difusion es 
muy limitada; recientes investigaciones numericas [67, 68] muestran que la extrapolacion de 
los resultados analiticos, obtenidos perturbativamente para pequefios niveles de turbulen- 
cia, no provee una descripcion apropiada de la difusion en condiciones de alta turbulencia. 
Complementando esas investigaciones, aquf calcularemos los coeficientes de difusion paralela 
y transversal, parametrizando los resultados mediante relaciones de escala que puedan ser 
luego utilizadas en diferentes aplicaciones. Siguiendo un procedimiento analogo, tambien 
evaluaremos numericamente el coeficiente antisimetrico, que es responsable, como ya hemos 
visto, de los efectos de drift, y que no ha sido estudiado previamente en condiciones de alta 
turbulencia. Mientras que en [67, 68] solo ha sido estudiado el caso de fluctuaciones en el 
campo magnetico random dadas por un espectro de Kolmogorov, aqui tambien consider- 
aremos otros tipos de turbulencia (a saber, los espectros de turbulencia de Kraichnan y de 
Bykov-Toptygin). Por otra parte, analizaremos la propuesta analitica de [66] (en donde se 
representa el efecto de la turbulencia mediante factores ad hoc de modulacion exponencial, 
ecs. ()2.22p - ()2.24p ) y discutiremos su rango de validez. 

Siguiendo el enfoque de [67], consideraremos particulas cargadas propagandose en un 
campo magnetico B(r) = SoZ + Br(r), donde el primer termino representa un campo regular 
uniforme dirigido segiin la direccion z, mientras que el segundo corresponde a la componente 
random. Un campo turbulento isotropico y espacialmente homogeneo puede ser aproximado 
numericamente mediante la suma de un numero grande (Nm) de ondas planas con direcciones 
del vector de ondas, polarizaciones y fases elegidas al azar [67,96], es decir, 

A^m 2 

Br(r) = E E ^(^-)^n cos(k„ -r + rn), (2-44) 

n=l a=l 

donde las dos direcciones de polarizacion ortogonales (a = 1,2) estan contenidas en 
el piano perpendicular a la direccion del vector de ondas (o sea, ± k„, con el fin de 
asegurar que V ■ B = 0). La distribucion de numeros de onda se toma de acuerdo a un 
espaciamiento logarftmico constante entre kmin = '^'^/Lmax J k^ax = "^T^l^min-, donde Lmm 
y ^max son las escalas de turbulencia minima y maxima, respectivamente. La densidad de 
energi'a de la componente turbulenta se toma como d£'r/d/c oc ^ donde el mdice espectral 
7 queda determinado por el tipo de mecanismo que genera la turbulencia. La amphtud de 
las ondas planas, en consecuencia, esta dada por J^ikn) = ■^{Br)k~'^{kn — kn-i), donde M 
es una constante de normalizacion tal que ^^^^(fcn) = {B^)- Como hemos mencionado 
mas arriba, aquf consideraremos tres tipos de espectro de turbulencia, a saber: el espectro 
de Kolmogorov con 7 = 5/3 (un caso particularmente atractivo, debido a que, de acuerdo 
con las observaciones [50,55], las fluctuaciones en la densidad del medio interestelar siguen 
este espectro de turbulencia), el espectro hidromagnetico de Kraichnan con 7 = 3/2 [97], y 
el espectro de Bykov-Toptygin con 7 = 2 [98]. 

El metodo de calculo consiste, basicamente, en generar una configuracion particular de 
campo magnetico random (eligiendo al azar la direccion de propagacion, polarizacion y 
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fase de las Nm ondas planas en la ec. ()2.44|1 ). y seguir la trayectoria de una parti'cula que 
se propaga desde el origen, con una direccion inicial arbitraria, en el campo magnetico 
total (es decir, en el campo que resulta de considerar la suma de componentes regular y 
random). Luego, los resultados deben ser promediados sobre un gran mimero de diferentes 
configuraciones de campo para calcular los coeficientes de difusion. Por ejemplo, en [67] 
los resultados fueron obtenidos integrando primero las trayectorias de 2500 particulas para 
una dada configuracion de campo, y luego promediando sobre 50 realizaciones de campo 
diferentes. Nosotros encontramos, sin embargo, que los resultados muestran una mayor 
dependencia con la eleccion particular de la configuracion de campo que con la eleccion de la 
velocidad inicial de la particula. En consecuencia, resulta mas conveniente promediar sobre 
un gran mimero de configuraciones de campo; en nuestro trabajo, generamos tipicamente 
~ 10^ configuraciones de campo, siguiendo la trayectoria de una linica parti'cula en cada 
configuracion, y usamos A^^, = 100 modos en todas las simulaciones. 

Para la integracion numerica de las trayectorias, utilizamos una rutina de Runge-Kutta 
con un paso temporal At = 0.1 r^/c. En [68], los coeficientes de difusion paralelo y perpen- 
dicular se calcularon a partir de la tasa de incremento asintotica (a tiempos grandes) de los 
desplazamientos cuadraticos medios en cada direccion, es decir. 



Sin embargo, este metodo requiere seguir las trayectorias de particula para tiempos con- 
siderablemente largos con el fin de alcanzar la region asintotica, tipicamente mas alia de 
t = 1000 tl/c. En cambio, nosotros encontramos mas conveniente calcular las funciones de 
descorrelacion de velocidad Rij(t), luego integrandolas en la formula de Kubo (|2.21|) . dado 
que este procedimiento requiere seguir las trayectorias solo hasta tiempos tipicamente no 
mayores a t = 100 r^/c. Por otro lado, este metodo permite calcular el coeficiente de di- 
fusion antisimetrico por medio de Ry^, mientras que promedios tales como {Ax Ay) se anulan 
a tiempos largos y, en consecuencia, no resultan utiles para calcular Da- 

La fig l2.16l muestra la dependencia temporal de las funciones de descorrelacion asociadas 
a diferentes componentes de la velocidad, para el espectro de Kolmogorov y diferentes nive- 
les de turbulencia Para baja turbulencia, se observa en la figEHEfa) que, luego de un 
perfodo transitorio, las funciones de descorrelacion se comportan segun las ecs. ()2.22|) - ()2.24|) . 
Esto permite calcular las escalas temporales de descorrelacion; por ejemplo, siguiendo los 
maximos locales de las funciones sinusoidales, se encuentra t±_ = r^, tal como muestra la 
fig ]2.17l En [66], estas escalas temporales de descorrelacion se consideraban iguales sola- 
mente como una hipotesis simplificadora, pero alii se sugen'a que podrian ocurrir nuevos 

Aqui definimos el nivel de turbulencia como = (i?^) /Bq [67] . Alternativamente, el nivel dc turbulencia 
puede definirse como rj = {B^)/{{B^) + B^) = a^/{l + a^) [68]. 




(2.45) 



y 




(2.46) 
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Figura 2.16: Dependencia temporal de las funciones de descorrelacion asociadas a diferentes 
componentes de la velocidad, para el espectro de turbulencia de Kolmogorov y dos niveles de 
turbulencia diferentes: (a) cr^ = 0.3 y (b) = 5, ambos correspondientes a = 0.1 Lmax- 



efectos en el caso general en que t± ^ ta- Aqui encontramos, en cambio, que efectiva- 
mente estas escalas temporales son iguales, como puede esperarse del heclio de que el origen 
fisico de la descorrelacion es el mismo en ambos casos. Para altos niveles de turbulencia, 
la descorrelacion es muclio mas abrupta, tal como muestra la fig ]2.16r b). Las funciones de 
descorrelacion se tornan despreciables mientras aun prevalece el periodo transitorio, y las 
ecs. ^^-^2M> dejan de ser validas. 

Definiendo el parametro adimensional p = ri/Lmax, los resultados, expresados en terminos 
de D / [cLmaxP) vs. p, son universales y satisfacen relaciones de escala que permiten calcular 
el tensor de difusion para diferentes conjuntos de valores para la amplitud del campo regular, 
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Figura 2.17: Maximos locales de \Rxx\ J \Ryx \ en funcion de t, graficados en una escala lineal- 
log, segun corresponde al caso de la fig J2.1fjf a). De aqm, resulta t± = ta (ver ecs. fl2.22|l - 
(im). 

las escalas de longitud caracteristicas asociadas a los modos del campo random, y la energia 
de las particulas que se propagan. For ejemplo, el rango 0.01 < p < 1 que estudiamos 
aqui puede considerarse como correspondiente a protones de energia 10^^ < E/eV < 10^^ 
propagandose en la galaxia (con Bq = 1 fiG and Lmax = 100 pc). Alternativamente, con- 
siderando protones propagandose en el campo interplanetario, de intensidad Bq = 50 /iG y 
Lmax = 0.01 AU, el li'mite al regimen difusivo para p = 1 corresponde a la energia cinetica 
Ek = 1.8 GeV. Dado que la dispersion del angulo de pitch ocurre fundament almente en 
condiciones de resonancia, la propagacion es esencialmente independiente de la escala minima 
de la turbulencia Lmin, en tanto sea Lmm ^l- Aqui tomamos ti'picamente Lmin = 0.1 r^. 

La fig l2.18l muestra los resultados numericos obtenidos para la difusion paralela, cor- 
respondientes al espectro de Kraichnan, para diferentes niveles de turbulencia en el rango 
0.1 < 0"^ < 100. Como resulta de las ecs. ()2.16|) - ()2.17|) . en el regimen de dispersion resonante 
se espera que D\\ oc p^"'*', mientras que a mayor rigidez (fuera del regimen de dispersion res- 
onante) es de esperar un decrecimiento en la efectividad de la dispersion segun E'"^ [3,59], 
conduciendo asi a D\\ oc p^, independientemente del espectro de turbulencia considerado. 
Los resultados que muestra la fig |2.18l concuerdan muy bien con el comportamiento espera- 
do; la region de transicion entre el regimen de dispersion resonante y el no resonante ocurre 
alrededor de p ~ 0.2. En efecto, la longitud de escala que separa ambos regimenes esta 

^°Sin embargo, en [67] el espectro de ley de potencias de la componente random del campo magnetico 
interplanetario se continiia con amplitud constante hasta una escala maxima de 1 AU, extendiendose con- 
sider ablemente el regimen difusivo. 
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Figura 2.18: Resultados correspondientes a la difusion paralela, para el espectro de Kraichnan 
y diferentes niveles de turbulencia en el rango 0.1 < < 100. 

determinada por la longitud de correlacion del espectro de turbulencia L^, definida como 

/oo 
dL(Br(0)-Br(r(L)))=L,(i?2) , (2.47) 
-oo 

donde el punto r(L) esta desplazado una distancia L con respecto al origen, a lo largo de 
una direccion fija. Entonces, considerando un espectro de fluctuaciones de mdice 7 que se 
extiende entre las longitudes de escala y Lmax, la longitud de correlacion queda dada 
por [100] 



-T ■ (2.48) 

2 ^ ^ — ( r IT 

I ^ y-'-' mm I ■'-'max j 

De aqui puede verse que Lc/Lmax — 0.2 resulta un valor representativo para la longitud 
de correlacion de los espectros de campo random que consideramos en este trabajo, y esto 
permite explicar el cambio de regimen observado en p ~ 0.2. 

De acuerdo con estas consideraciones, una manera apropiada de ajustar los resultados 
es interpolando D\\ entre las leyes de potencia que caracterizan los regimenes de baja y alta 
rigidez. Esto puede Uevarse a cabo mediante la expresion 



cL 

r 




(2.49) 



donde los parametros N\\ y py estan dados en la Tabla 12.21 Esta expresion ajusta bien 
nuestros resultados numericos hasta cr^ ~ 10, que corresponde a un nivel de turbulencia 
sobradamente alto para la mayor parte de las aplicaciones astrofisicas de interes. Para 
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Espectro 


7 


iV|j 


Pll 




asint 
P\\ 




a± 


Na 


Kraichnan 


3/2 


2.0 


0.22 


3.5 


0.65 


0.019 


1.37 


17.6 


Kolmogorov 


5/3 


1.7 


0.20 


3.1 


0.55 


0.025 


1.36 


14.9 


Bykov- Toptygin 


2 


1.4 


0.16 


2.6 


0.45 


0.020 


1.38 


14.2 



Tabla 2.2: Parametros de las formulas de ajuste correspondientes a los coeficientes de difusion 
paralelo, transversal y antisimetrico, para diferentes tipos de espectros de turbulencia. Ver 
mas detalles en el texto. 

niveles de turbulencia aiin mayores, los resultados quedan mejor descriptos asintoticamente 
por los parametros A^|^**"* y tambien dados en la Tabla 1?!^ 
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Figura 2.19: Comparacion entre el ajuste a D\\ dado en este trabajo (valido formalmente 
solo en el rango 0.01 < p < 1) y los resultados de [68], ambos correspondientes a un espectro 
de Kolmogorov, y para los mismos niveles de turbulencia. 

Con el fin de chequear la consistencia con los resultados numericos previos para la difusion 
paralela en condiciones de alta turbulencia (que fue estudiada en [68] asumiendo un espectro 
de fluctuaciones de Kolmogorov) , la fig J2.19l muestra una comparacion entre nuestra formula 
de ajuste y los resultados de [68], ambos correspondientes a un espectro de Kolmogorov, 
y para los mismos niveles de turbulencia. Debe notarse que en [68] el radio de Larmor^^ 
r2 se define reemplazando Bq + en la ec. fl2.5j) . es decir, acoplando en la 

dependencia con la rigidez y con el nivel de turbulencia, mientras que, por otro lado, el 
parametro adimensional p* se define como p* = 27rr^/ Lmax- En consecuencia, los resultados 
deben ser reescaleados de acuerdo con las relaciones D/{cLmaxP) = -D/ (r^c)/ a/I + o^ y 

^^Para evitar confusiones, adoptamos aqui un asterisco para referirnos a las cantidades definidas en [68]. 
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Figura 2.20: Ajuste a D±/ D\\ como funcion de la rigidez para el espectro de Kolmogorov en 
el regimen de baja rigidez (es decir, p < 0.2), para diferentes niveles de turbulencia. Para 
mayor rigidez, D±/D\\ oc p^^ (ver la ec. ()2.50|) ). Para comparar, se muestran tambien los 
resultados que se obtienen de ajustes a los datos dados en [67] y [68]. 



p = p*vT+o^/27r. Como puede verse en la fig ]2.191 ambos conjuntos de datos estan en 
buen acuerdo. 

La manera mas simple de parametrizar el coeficiente de difusion perpendicular es a traves 
de los resultados para el cociente D±/D\\, dado que este cociente exhibe una dependencia 
sencilla con la rigidez. Con el fin de evitar el regimen subdifusivo (que discutimos brevemente 
mas adelante), nos restringimos a la region delimitada por p > 0.03 y cr^ > 1. Como 
en las investigaciones numericas previas [67,68], encontramos que el cociente D±/D\^ es 
independiente de la rigidez en la region de baja rigidez (es decir, para p < 0.2), mientras que 
exhibe una dependencia segiin en la region de alta rigidez. En consecuencia, nuestros 
resultados pueden ser parametrizados convenientemente por medio de la expresion 



Dii 



N± X {ay^ X 



1 (p < 0.2) 

ip/0.2)-' (p > 0.2) 



(2.50) 



donde los parametros A^_l y a± estan dados en la Tabla lT^ Hemos visto que los datos de [6^ 
(para el regimen de baja turbulencia) pueden representarse mediante 



1 + 4.52/ ((J 



2^1-5 



(2.51) 



(ver la ec. ()2.27|) y el parrafo que le sigue); analogamente, los resultados que se presentan 
en [67], correspondientes al rango de turbulencia 0.03 < ct^ < 1, pueden ajustarse por medio 
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Figura 2.21: {{Ax^f) /2At/ (cLmaxP) en funcion del tiempo, para diferentes valores de rigidez 
y niveles de turbulencia, en el caso de un espectro de fluctuaciones de Kolmogorov. El 
regimen subdifusivo aparece en la region de rigidez y nivel de turbulencia pequenos. 

de la expresion 

^ = 0.025 X (^2)1-835 _ (2.52) 

La fig l2.20l muestra el ajuste dado por la ec. ()2.5()j) para el caso de Kolmogorov y diferentes 
niveles de turbulencia. Para comparar, tambien se muestran los resultados de [68] (segun 
se obtienen de la ec. ()2.51|) ) y de [67] (dados por la ec. ()2.52|) ). Para niveles de turbulencia 
moderados (a^ = 1) nuestros resultados muestran un acuerdo excelente con los de [67], 
mientras que el acuerdo con [68] tiende a mejorar para turbulencias mayores. 

Para rigidez y nivel de turbulencia pequeiios (por ejemplo, p ~ 0.01 y o"^ < 1), hemos 
encontrado evidencia del fenomeno de subdifusion, en concordancia con lo reportado en [68]. 
Una particula de rigidez pequefia, propagandose en condiciones de baja turbulencia, tiende 
a permanecer siempre sujeta a una misma Ifnea de campo; entonces, la difusion transversal 
ocurre fundamentalmente debido a la caminata aleatoria transversal de las propias Imeas de 
campo. En consecuencia, se espera que el desplazamiento perpendicular medio evolucione 
con el tiempo mas lentamente, en comparacion con el caso de la difusion normal (es decir, 
{AxD oc t™, con m < 1). La fig HTCl muestra {{Ax±) )/2At/{cLmaxp) como funcion del 
tiempo, para diferentes valores de rigidez y nivel de turbulencia, en el caso de un espectro de 
fluctuaciones de Kolmogorov. Un plateau se obtiene a tiempos grandes para p = 0.1 y o"^ = 1 
(correspondiendo, por lo tanto, a la relacion de difusion usual (Ax^) oc t), mientras que, para 
valores mas pequenos de rigidez y nivel de turbulencia, aparece el fenomeno de la subdifusion. 
Como una primera aproximacion a la subdifusion, se ha propuesto e investigado el caso 
li'mite conocido como "difusion compuesta", en donde las particulas estan estrictamente 
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Figura 2.22: Comparacion entre el ajuste a D^/p obtenido aqui (para el espectro de Kol- 
mogorov y los valores de rigidez indicados) y la prescripcion independiente de p adoptada 
previamente (ecs. ()2.28|) y ()2.55|) ). 



restringidas a moverse a lo largo de las Imeas de campo [99, 101-103]. Para el caso de 
la difusion compuesta, se encuentra que m = 1/2, mientras que, para el transporte de 
particulas en tres dimensiones, se espera que m tenga una dependencia suave con la rigidez 
y la turbulencia, tal que 1/2 < m{p, cx^) < 1. 

Con respecto al coeficiente de difusion antisimetrico, los resultados numericos pueden 
ajustarse adecuadamente mediante la expresion 



Da 1 



cLmaxP 3 + (^2/^2)2 

donde 



(2.53) 



y donde los valores para el parametro Na estan dados en la Tabla 12.21 En el limite de 
turbulencia muy baja, esta expresion tiende apropiadamente al valor Da — ctl/S [59] (ver 
tambien la ec. ()2.20|) en el limite de ta grande), mientras que, como es de esperar, se anula 
en el limite de turbulencia muy alta. Anteriormente, habiamos adoptado una simple pre- 
scripcion para el coeficiente antisimetrico (ver la ec. ()2.28|) ): en terminos de cr^, la expresion 
que habfamos usado resulta 

Da 1 1 , , 

(2.55) 



cL^axP 3 1 + ((t2)^-V4.52 
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Puede verse que, segun esta prescripcion, Da/ cLmaxP es independiente de la rigidez. En la 
fig ]2.22l se compara esta prescripcion con el ajuste dado por las ecs. fj2.53p - ()2.54p . calculado 
para el caso de Kolmogorov y diferentes valores para la rigidez. Se observa un excelente 
acuerdo entre la prescripcion ()2.55|) y el ajuste correspondiente a p = 0.01; la dependen- 
cia de Da/ p con p es, en cualquier caso, solamente significativa para niveles de turbulencia 
muy altos (a^ ^1). Finalmente, resulta importante destacar que la adopcion de las nuevas 
expresiones que presentamos aqui para los coeficientes de difusion en condiciones de alta tur- 
bulencia no modifica sustancialmente los resultados obtenidos en las Secciones precedentes. 



Capitulo 3 



Neutrinos de muy alta energia 
inducidos por rayos cosmicos en 
escenarios dependientes de la rigidez 

En este Capitulo, estudiainos el flujo de neutrinos de altas energias en escenarios 
en los que la rodilla del espectro de rayos cosmicos es un efecto dependiente de la 
rigidez magnetica; en esta clase de escenarios, el fondo de neutrinos atmosfericos 
se reduce significativamente, facilitando la husqueda de fuentes astrofisicas. Asi, 
se pone en evidencia la necesidad de resolver el origen de la rodilla, con el fin de 
interpretar correctaniente las observaciones en los nuevos telescopios de neutri- 
nos. Luego de investigar este fenomeno en el canal usual de trazas leptonicas, 
asociado a la deteccion de neutrinos muon, proponemos un nuevo metodo para 
aislar el flujo de neutrinos atmosfericos prompt a traves del canal de las cascadas 
hadronica y electromagnetica, asociado a la deteccion de neutrinos electron. 

3.1 Introduccion 

La existcncia de campos magneticos en la region de propagacion de los rayos cosmicos tiene 
como efecto la alteracion de sus trayectorias, dificultando, en consecuencia, la identificacion 
de fuentes puntuales. Mas alia del regimen difusivo, que hemos estudiado en detalle en 
el Capitulo anterior, a mayores energias se hacen presentes otros fenomenos, asociados al 
lensing magnetico (la magnificacion o atenuacion de las fuentes observadas, el desplazamiento 
de sus posiciones aparentes, la aparicion de pares de imagcnes adicionales, etc. [100, 104- 
106]); eventualmente, la astronomia de rayos cosmicos podria resultar dificil de ser llevada 
a la practica. Las mismas fuentes astrofisicas de altas energias podrfan producir tambien 
rayos 7, que conservan la direccion de propagacion, pero que, para altas energias y grandes 
distancias, serian absorbidos por la reaccion 77 — > e'^e~ con el fondo infrarrojo de radiacion 
(por ejemplo, el camino libre medio es de 100 Mpc para fotones con energias alrededor de 
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10^ GeV [107,108]). En muchos modelos de fuentes de altas energi'as, tambien los neutrinos 
son producidos en abundancia. Los neutrinos tienen la ventaja de no ser deflectados ni 
absorbidos aun recorriendo vastas distancias; ademas, son capaces de escapar de fuentes de 
gran densidad. Una desventaja evidente es que, precisamente debido a su elusivo caracter, 
resultan tambien muy dificiles de detectar, ya que solo interactuan debilmente. 

Sin embargo, por primera vez, los tclcscopios de neutrinos actualmente en operacion o 
en construccion tcndran la sensibilidad rcquerida para poner a prueba modelos rcalistas de 
las fuentes de mayor energia del Universo [109-113]. Por ejemplo, provenientes de varios 
nucleos galacticos activos (AGN) cercanos, fueron detectados rayos 7 de alta energia hasta 
alrededor de 10^ GeV [114-117]. Si estos rayos 7 estan originados en el decaimiento de piones 
neutros, entonces deben espcrarsc flujos similares de neutrinos originados en el decaimiento 
dc piones cargados. La produccion dc piones resulta natural cn modelos en los que un flujo de 
protones de altas energias interactiia, en la region de las fuentes, con otros nucleones, o bien 
con fotones en el campo de radiacion ambiente. El detector AMANDA esta comenzando a 
poner a prueba estos modelos a un nivel competitivo con los telescopios de rayos 7 [118-124]. 

Mas alia de las fuentes puntuales, los tclcscopios dc neutrinos tambien pueden medir el 
fondo difuso que proviene de fuentes mas distantcs y de mayor energia, que no son direc- 
tamente visibles mediante rayos 7 debido a la opacidad del fondo cosmico infrarrojo. Sin 
embargo, sera un gran desafio poder distinguir un fondo extragalactico difuso del flujo de 
neutrinos producido por la colision de los rayos cosmicos con la atmosfera terrestre. El 
espectro de neutrinos atmosfericos decrece segiin E^'^, con el indice espectral en el rango 
7 ~ 3 — 3.7, mientras que, para la componente extragalactica, se asume usualmente un 
espectro E'"^. En consecuencia, una componente no atmosferica podrfa ser descubierta a 
partir de un quiebre en el espectro medido. Debajo del quiebre, el espectro estan'a dominado 
por el fondo, y arriba del quiebre, por la serial. 

Los neutrinos atmosfericos ban sido bien identificados a bajas energias por Super - 
Kamiokandc y otros detectores [125, 126], y recientemente tambien fueron medidos a mayores 
energias por AMANDA [127]. En efecto, el estudio de neutrinos atmosfericos dc energias 
sub-GeV y multi-GeV ha sido de enorme importancia; en particular, ha provisto la primera 
evidencia clara en favor de las oscilaciones de neutrinos, demostrando que su masa es no 
nula. Estos neutrinos se producen principalmente en el decaimiento de piones y kaones 
(por ejemplo, 7r~ con yU~ e'VeV^) que son producidos por rayos cosmicos que 

incidcn sobre la atmosfera, generando Uuvias atmosfericas con una componente hadronica 
significativa. 

La vida media de los mesones aumenta con la energia, por el efecto relativista de dilat- 
acion temporal; debido a la interrelacion entre el decaimiento de los mesones y su interaccion 
en el aire, a mayores energias los mesones sufren una atenuacion mayor antes de decaer. De 
un modo similar, los muones producidos en la atmosfera pucdcn propagarse sin decaer y 
alcanzar la superficie terrestre. Los efectos relativistas explican que el indice asociado al 
espectro de neutrinos atmosfericos sea mayor al correspondiente a los rayos cosmicos; un es- 
pectro de rayos cosmicos con indice ^rc ~ 2.7 produce un espectro de neutrinos atmosfericos 
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caracterizado por un mdice en el rango mcncionado, 7 ~ 3—3.7. Dc la longitud dc dccaimien- 
to de los mesones cargados, L = ^ct, donde, en particular, L^^ ~ 5.6 km(£'7r/100 GeV) y 
Lk — 7.5 km(£^;^/TeV), resulta que, por encima de 100 GeV, los piones atraviesan varias 
longitudes de atenuacion en la atmosfera antes de decaer ^; en consecuencia, a estas energias 
los neutrinos son producidos principalmente en el decaimiento de kaones. A la vez, para en- 
ergias por encima del TeV, los kaones tambien se atenuan significativamente antes de decaer, 
y el decaimiento de partfculas charm se convierte en el canal de produccion de neutrinos mas 
relevante. Esta serial, que recibe el nombre de neutrinos atmosfericos prompt charm, no ha 
sido aiin identificada experimentalmente. Por otra parte, existen grandes incertezas desde el 
lado teorico. La produccion de partfculas charm por rayos cosmicos requiere tener en cuenta 
procesos NLO (que incrementan la produccion de partfculas charm por un factor mayor a 2) 
y secciones eficaccs que involucran el uso de funciones dc distribucion partonica para valores 
muy pequenos de la variable dc Bj0rken x, mas alia del rango explorado experimentalmente. 

En el momento en que aparezca, eventualmente, un quiebre en el espectro medido por 
los telescopios de neutrinos, debera analizarse cuidadosamente si la observacion corresponde, 
efectivamcntc, al descubrimicnto dc una nucva componcntc en el flujo de neutrinos de muy 
alta energia, o bien si se trata solo de una fluctuacion. En este sentido, disponer de diferentes 
canales de deteccion sera muy valioso. Si la senal es real, podria constituir un importante 
indicio de nueva fisica, en uno de dos sentidos: (i) como el flujo de neutrinos atmosfericos 
prompt (y, por lo tanto, como un nuevo canal que provee informacion sobre rayos cosmicos 
y QCD), o (ii) como un flujo de neutrinos cxtragalacticos (que daria nucvos indicios sobre 
el Universo de altas energias). Distinguir entre estas posibilidades tambien requerira el uso 
complement ario de diferentes canales de deteccion. 

Las investigaciones sobre los telescopios de neutrinos se han enfocado, hasta aqui, princi- 
palmente hacia el canal de deteccion de interacciones de corriente cargada de V/^, que involucra 
la medicion de muones desplazandose hacia arriba. Dado que, por encima de algunos cicntos 
de GeV, el rango de los muones en hielo excede 1 km (la dimension lineal caracteristica 
de IceCube), el volumcn dc deteccion efectivo no es igual al volumen del detector, sino, en 
cambio, al area del detector multiplicada por el rango de los muones, que es una funcion 
que crece con la energia. Este efecto, combinado con el incremento de la seccion eflcaz de 
interaccion de los neutrinos, tiende a compensar el decrecimiento del flujo de los neutrinos 
incidentes. 

En la Seccion 3.2, analizarcmos en dctallc el canal de trazas Icptonicas, asociado a la 
deteccion de neutrinos muon, y enfocaremos nuestra atencion en las componentes de neu- 
trinos atmosfericos. Si bien este canal de deteccion ha merecido la atencion de numerosas 
investigaciones previas, aqui harcmos notar que el flujo de neutrinos de altas energias (para 

> 10^^ eV) depende significativamente del escenario asumido para explicar la rodilla del 
espectro de rayos cosmicos. En particular, contrastaremos los resultados que se obtienen 
de suponer que los rayos cosmicos consisten solo de protones (una hipotesis simplificadora 

^Una longitud de atenuacion a esta energia corresponde aproximadamente a 120 g/cm^. 
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usualmente asumida en los trabajos previos) con las prcdiccioncs que resultan de escenar- 
ios mas rcalistas, en los que la rodilla es un efecto dependiente de la rigidez magnetica. 
Como veremos, en esta clase de escenarios, el flujo difuso de neutrinos de altas energias 
producido en la incidencia de rayos cosmicos sobre la atmosfera se reduce significativamente; 
en consecuencia, su deteccion en los observatorios de neutrinos se torna mas dificil, pero a 
la vez facilita la biisqueda de fuentes astrofisicas. Asi, ponemos en evidencia la ncccsidad 
de resolver cl origcn dc la rodilla del espcctro dc rayos cosmicos, con el fin dc intcrprctar 
correctamcnte las observaciones de neutrinos de muy alta energia en los nuevos telescopios. 

Luego, en la Seccion 3.3, estudiaremos el canal de las cascadas hadronica y electro- 
magnetica, asociado a la deteccion de neutrinos electron. AUi proponemos un nuevo metodo 
para aislar cl flujo dc neutrinos atmosfericos prompt, cl cual, como ya hemos mencionado, 
es importantc tanto cn si mismo, como en su condicion dc fondo dc la serial cxtragalactica. 
Nuestra propuesta esta enfocada en el canal de eventos de corriente cargada de v^.^ enfatizan- 
do la importancia de considerar el espectro de eventos en terminos de la energia detectable, 
en lugar del producto de flujo por seccion eflcaz como funcion de la energia del neutrino. 
Aunque este canal ha sido, por lo general, ignorado en la literatura, tiene importantes ven- 
tajas. La fraccion de en el flujo atmosferico prompt y en el flujo extragalactico es mayor 
que para el flujo atmosferico convencional de altas energias. Veremos que, en este canal, el 
quiebre espectral ocurre a energias aproximadamente un orden de magnitud menor que para 
el canal de trazas leptonicas; esta es una ventaja importante, porque, a bajas energias, los 
flujos son mayorcs y los cfcctos dc absorcion terrestre son menorcs. Varios autores se ban 
centrado en la deteccion de sin embargo, a energias por debajo de 5 x 10^ GeV, resulta 
dificil separar las interacciones individuales de de aquellas de otros sabores. El canal v^. 
deberia ser particularmente efectivo para los neutrinos atmosfericos prompt, debido a que 
su espectro decrece mas rapidamente que el espectro extragalactico. De este modo, si el 
quiebre espectral ocurre a menor energia, cl flujo atmosferico prompt pucdc dctectarse mas 
facilmente. Por otra parte, la fldelidad espectral cntre la energia del neutrino y la energia 
detectada es mucho mayor que en el caso de interacciones de z/^ de corriente cargada, algo 
importante para la biisqueda de un quiebre espectral. Siguiendo el enfoquc dc la Seccion 
precedente, estudiaremos el efecto que tiene la composicion asumida para los rayos cosmicos, 
comparando nuevamente el caso en el que los rayos cosmicos consisten solo de protones, con 
escenarios en los que la rodilla es un efecto dependiente de la rigidez magnetica. Finalmente, 
tambicn estudiaremos el flujo difuso dc neutrinos galacticos, producidos por la intcraccion dc 
los rayos cosmicos con el medio interestelar, analizando sus componentes y su importancia 
relativa. 
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3.2 El canal de trazas leptonicas: neutrinos muon 
3.2.1 El espectro de rayos cosmicos 

Una manera sencilla de describir el flujo de rayos cosmicos galacticos en un escenario depen- 
diente de la rigidez es asumir que cada componente de carga Z esta dada por 



(3.i: 



donde 0| (0|) es el flujo de rayos cosmicos a bajas (altas) energias ^, dado por 



En esta expresion, 0o es la normalizacion del flujo total y fz la abundancia fraccional, 
ambos correspondientes a la energia de referenda Eq. Eligiendo Eq = 1 TeV, la normalizacion 
y las abundancias son las dadas en la Seccion 2.2; alli tambien estan tabulados los indices 
espectrales az (ver la Tabla ITTj) . Solamente dos parametros quedan libres: Aa, el cambio 
de fndice espectral a traves de la rodilla (que es el mismo para todas las componentes), y 
Ek, que flja la posicion de la rodilla (y se elige convenientemente para ajustar el espectro 
total a las observaciones) . La ec. ()3.1|) provee una interpolacion suave ^ entre el flujo 0| a 
bajas energias {E < ZEk) y 0^ a energias altas {E > ZEk). 

Asumiremos, ademas, una componente extragalactica dada por la ec. ()2.37|) (ver la 
Seccion 2.3), y formada solo por protones ^. La flg l3.1l muestra el ajuste a las observa- 
ciones experiment ales, segiin resulta de tomar Aa = 2/3 (correspondiente, por ejemplo, 
al modelo de la difusion turbulenta y drift desarrollado en el Capitulo 2), tomando el 
parametro Ek = 2.2 x lO^^eV, y Aa = 2 (correspondiente al modelo asumido en [72]), 
con Ek = 3.1 X lO^^eV. Mientras el caso Aa = 2/3 exhibe un acuerdo excelente con las 
observaciones, puede verse que el caso Aa = 2 conduce a la supresion abrupta del flujo total 
por encima de ~ 10^^ eV; en lo que sigue, enfocaremos nuestra atencion en el caso Aa = 2/3. 



3.2.2 El flujo de neutrinos muon atmosfericos 

Los rayos cosmicos que inciden sobre la atmosfera terrestre producen flujos secundarios de 
hadrones y leptones, que pueden describirse por medio de un sistema acoplado de ecuaciones 

^En este escenario, el espectro asociado a cada componente tiene una rodilla, cuya energia es proporcional 
a la carga; en relacion con ella, quedan definidas las regiones de bajas y altas energias que mencionamos 
aqui. 

■^Eventualmente, la ec. H3.1|l puede generalizarse conio = + donde el parametro 

n regula el "ancho" de la rodilla; tomando n ^ 1, resulta un cambio de pendiente muy abrupto. Para el 
caso Aa — 2/3 que consideramos aqui, n — 1 conduce a un acuerdo excelente con los datos experimentales 
(ver fiernmi. 

^De todas formas, la composicion de la componente extragalactica de rayos cosmicos solo afecta las 
predicciones para neutrinos de energia E^, > 10^^ eV. 
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Figura 3.1: Espectros de rayos cosmicos en escenarios dependientes de la rigidez, segiin la 
parametrizacion dada por las ecs. ()3.1|) - ()3.2|1 . Para el cambio de mdice espectral Aa = 2/3, 
tambien se muestran las principales contribuciones al flujo total, correspondientes a niicleos 
de H, He, O y Fe. La Imea recta punteada corresponde al flujo de rayos cosmicos extra- 
galacticos (ec. ()2.37j) ). Para comparar, tambien se muestran las observaciones experimen- 
tales. 

diferenciales que recibe el nombre de "ecuaciones de cascada". Bajo ciertas hipotesis sim- 
plificadoras, las ecuaciones de cascada pueden ser resueltas analfticamente. En esta Seccion, 
investigamos las componentes atmosfericas de los fiujos de neutrinos y antineutrinos muon. 
Debe notarse que, al no distinguir particulas y antiparticulas, aqui nos referimos, en todo 
momento, a la suma de neutrinos y antineutrinos. 

En el rango de interes de este trabajo (10^ GeV< E,, < 10* GeV), las principales contribu- 
ciones al flujo de neutrinos provienen del decaimiento de los mesones vr^, K"^, Kg y K^, que 
producen los neutrinos atmosfericos convencionales que dominan a bajas energias [128-140], 
y del decaimiento de los mesones charm D^.D^.K^ y Dg, que dan lugar al flujo de neu- 
trinos atmosfericos prompt charm que domina a energias altas [141-151]. La contribucion 
de los muones al flujo de neutrinos atmosfericos es completamente despreciable en este con- 
texto; debido a efectos relativistas, por encima de ii^ ~ 10^ GeV los muones Uegan hasta la 
superflcie terrestre sin decaer. Finalmente, debe tenerse en cuenta que, para las energias con- 
sideradas aqui, las oscilaciones de sabor no tienen ningun efecto apreciable en los neutrinos 
atmosfericos. 

La ecuacion de transporte asociada a la componente j— esima de la cascada puede es- 
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cribirse como 

= -^{E,X) - |(E,X) + Y,S,,{E,X) , (3.3) 

donde X es el slant depth definido por la ec. ()1.17p (ver la Seccion 1.5), Xj la longitud de 
interaccion en el aire y la longitud de decaimiento. Estas dos longitudes caracteristicas son 
medidas en unidades de g/ cm^ (es decir, incluyen el factor p{X) que corresponde a la densidad 
local de la atmosfera). El ultimo termino corresponde al proceso de produccion/regeneracion 
dado por 

S,,(E,X)- j^AE^-^—^ , (3.4) 

que acopla las ecuaciones de transporte asociadas a las diferentes componentes de la cascada. 
Puede notarse que akA aqui se refiere a la seccion eficaz total fc+aire, mientras que akA^j 
corresponde al proceso fc+aire — > j+ Como las energfas de ligadura entre los nucleones 
de un nucleo son mucho menores que las energfas tfpicamente involucradas en las colisiones 
nucleares entre rayos cosmicos y nucleos atmosfericos, en este contexto solo es relevante la 
interaccion nucleon-nucleon. En consecuencia, las diferentes componentes de la cascada que 
tomaremos en cuenta son la componente de nucleones (que agrupa a protones y neutrones) 
y las diferentes componentes de mesones M = vr^, K^, K^, D^, D^, Ac y Ds- Como en [133], 
la contribucion de Kg puede incluirse en el flujo de piones. Por simplicidad, consideraremos 
ademas que los linicos procesos relevantes son aquellos que corresponden a la regeneracion 
(es decir, Snn J Smm) y a la produccion de mesones por nucleones (es decir, Snm) [133,142]. 

Consideremos, en particular, flujos de la forma 4>k{E, X) = E~^''g{X) (es decir, planteamos 
la separacion de variables en los flujos, y escribimos el factor dependiente de la energia como 
una ley de potencias). Entonces, resulta que 

5..(i^,^) = ^^^^^..(i^), (3.5) 
donde Zkj es el momento de produccion/regeneracion definido por 

Z,AE) = /'dx ifM^^ , ,3.8) 

Jo (TkAiE) dx 

con X = E/E'. 

Asumiendo que las longitudes de interaccion son independientes de la energia, y que 
la distribucion diferencial de produccion satisface el scaling de Feynman, los propios mo- 
mentos de produccion/regeneracion se tornan independientes de la energia. Asi, para las 
longitudes de interaccion relevantes y los momentos de produccion/regeneracion que involu- 
cran nucleones y mesones no charm, consideramos los valores independientes de la energia 
dados en [2,133]. Ademas, tomamos la dependencia de los momentos con Pk de ajustes a 
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los datos compilados en [2]. Por otra parte, las longitudes de interaccion y los momentos de 
regeneracion asociados a todas las componentes charm pueden tomarse como iguales a las 
correspondientes a los kaones [144]. 

Sin embargo, para la produccion de particulas charm debe tenerse en cuenta el momento 
dependiente de la energia dado por la ec. (j3.6|) . ya que los resultados son muy sensibles a 
la seccion eficaz que se adopta para describir el proceso de produccion charm. En efecto, 
las diferentes estimaciones de los flujos atmosfericos prompt calculados hasta el momento 
difieren en mas de un orden de magnitud [149]. 

Debido a su masa relativamente grande, se considera usualmente que los quarks charm 
son producidos en procesos duros que pueden ser bien descriptos por la cromodinamica 
cuantica perturbativa (pQCD). A orden dominante (LO) en la constante de acoplamiento, 
los procesos que contribuyen a la seccion eficaz de produccion charm son el proceso de fusion 
gluon-gluon gg ^ cc j el proceso de aniquilacion quark-antiquark qq — > cc. Sin embargo, 
en el siguiente orden de aproximacion perturbativa (NLO) aparece el proceso de dispersion 
de gluones gg gg, aportando una contribucion muy significativa a la seccion eficaz total, 
que de esta manera se incrementa en un factor ~ 2 — 2.5 [152-154]. Con respecto a cor- 
recciones de mayor orden, se espera que su contribucion sea relativamente pequena, dado 
que cualitativamente no hay nuevos canales que puedan aparecer. Varias funciones de dis- 
tribucion partonica han sido propuestas, demostrando que proveen predicciones teoricas de 
pQCD que reproducen razonablemente bien los datos disponibles de los aceleradores. Sin 
embargo, para calcular el fiujo atmosferico charm se necesita extrapolar las funciones de 
distribucion de partones hasta fracciones de momento partonico muy pequenas, tfpicamente 
a: ~ 4 GeV/E [148]. Para las energias que debemos considerar aqui, esto queda fuera de la 
region accesible experimentalmente {x > 10~^), y la incerteza en la extrapolacion afecta los 
resultados finales significativamente [144,149]. Con el fin de ilustrar este rango de incertidum- 
bre, mostramos los resultados que se obtienen con dos funciones de estructura diferentes, a 
saber la funcion de distribucion partonica CTEQ3 [155] y el modelo Golec-Biernat, Wiisthoff 
(GBW) [156,157], que incluye efectos de saturacion de gluones. Las secciones eficaces de 
produccion de charm correspondientes fueron obtenidas, ya sea directamente del ajuste dado 
en [148] (para los resultados basados en el modelo de saturacion GBW), o bien ajustando 
los resultados dados en [144] para las funciones de estructura CTEQ3 (evaluadas con los 
parametros M = 2mc = 2/i, con la masa del quark charm rric = 1.3 GeV), e interpolan- 
do los resultados para diferentes energias. Para la seccion eficaz nucleon-aire, usamos la 
parametrizacion dada por [158], 



Una vez que se determinan todas las longitudes de interaccion y los momentos de pro- 
duccion/regeneracion relevantes, los fiujos atmosfericos de nucleones y mesones pueden cal- 
cularse de las ecuaciones acopladas de cascada dadas por ()3.3p y ()3.5p . Recordando que 
estas expresiones fueron obtenidas suponiendo que los nucleones y mesones tienen espectros 




(3.7) 
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de ley de potencias, discutiremos, en primer lugar, los resultados para un flujo inicial de 
nucleones dado por 4>n{E, X = 0) = (po^E"'^, y dejaremos la discusion del caso general para 
mas adelante. 

Bajo las suposiciones que hemos mencionado, el flujo de nucleones se desarroUa indepen- 
dientemente de los flujos de los mesones, y resulta 

ME, X) = e-^/^^0o^E-^ , (3.8) 

donde la longitud de atenuacion de los nucleones se deflne como 

= ■ (3.9) 

Con respecto a las ecuaciones de cascada de los mesones, estas se resuelven usualmente 
considerando por separado la solucion de bajas energias 0^ (que desprecia los terminos de 
interaccion y de regeneracion, dado que \m ^ a bajas energias) y la solucion de altas 
energias 0^ (que, en cambio, desprecia el termino de decaimiento, pues Am ^ A^^ a energias 
muy altas) [2,133,142]. 

La solucion de altas energias esta dada por 



bUE.X) = ^7 ' , M ' (3.10) 

J- — ^NN -L — Aat/ Am 



UE.X) = -^^^e-^/^-0o^E~- . (3.11) 



donde la longitud de atenuacion de los mesones Am se deflne analogamente a la ec. (j3.9p . 

Para el caso de bajas energias, la solucion se obtiene de reemplazar Am — >■ A^ en la ultima 
ecuacion. Debido al factor de Lorentz de dilatacion temporal, la longitud de decaimiento es 
lineal con la energia; luego, vemos que la expresion resultante acopla la dependencia con la 
energia y el slant depth, y no conduce a un simple espectro de ley de potencias. Sin embargo, 
para valores de X no demasiado pequefios, la solucion se reduce a 

mN A AT 

En las ecs. (j3.8|) - (j3.1H) . Z]\rN, Zj^m, J Zmm deben ser evaluados con el indice espectral 
7, correspondiente al flujo de nucleones. Puede notarse tambien que el flujo de mesones 
para altas energias tiene el mismo indice espectral que el flujo de nucleones, mientras que la 
solucion de mesones para bajas energias da lugar a un espectro mas duro, debido a un factor 
adicional (lineal en la energia) implicito en Xf^. 

El flujo de neutrinos producido por el decaimiento debil de mesones puede calcularse 
siguiendo el mismo esquema que ya hemos desarroUado para obtener el flujo atmosferico 
de nucleones y mesones. En efecto, el flujo de neutrinos esta dado por una ecuacion de 
transporte (analoga a la ec. ()3.3p ) que solo contiene los terminos de fuente que resultan del 
decaimiento de mesones, es decir, 

^(E,X) = J2SmAE,X) , (3.12) 

M 
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donde 



La distribucion de decaimiento puede ponerse en terminos de Fmu, el espectro inclusivo de 
neutrinos en el decaimiento del meson M, por medio de la relacion 

1 dT,UE,E')_^^^^^^^^^^^^^^^,^^ ^3^^^^ 



Tm{E) dE 



donde Bpf^, es el branching ratio asociado al decaimiento del meson M en un estado con el 
dado neutrino u. En el Ifmite ultrarrelativista, el espectro inclusivo de neutrinos satisface la 
relacion de escala Fmu^E, E') = Fmu{E / E') / E' [2,133]. Por otro lado, ya hemos derivado 
las soluciones asintoticas para el flujo de los mesones (ecs. ()3.1()j) y 1)3.111) ). que son de la 
forma (Pm{E, X) oc E^^g{X), con /3 = 7 (/? = 7 — 1) para los flujos de altas (bajas) energfas. 
De esta manera, los terminos de fuente dados por la ec. ()3.13|) pueden reescribirse como 



SM.iE,X) = ^^^^^^Z^,l\E), (3.15) 

donde los momentos de decaimiento de los mesones (dependientes de /5) son definidos medi- 
ante ^ ^ 

Zl,^{E) = Bmu [ dx x^-'FMuix) . (3.16) 
Jo 

Los momentos de decaimiento que usamos en este trabajo fueron calculados ajustando la 
dependencia con P de los momentos relevantes tabulados en [142] para 2.7 < /5 < 4, que, a 
la vez, fueron determinados mediante los programas Lund de simulacion Monte Carlo. 

Entonces, los flujos de neutrinos asintoticos para baj as y altas energfas pueden ahora ser 
determinados de las ecs. ()3.12p y ()3.15|) . que involucran a los momentos de decaimiento de 
los mesones a y los correspondientes flujos de mesones para altas y bajas energfas, dados por 
las ecs. fl3.10|) y ()3.11|) . La solucion de bajas energfas resulta 



^AE,X) = Zl^^J^ (1 - e-^/^-) 0o^E-^ . (3.17) 



Como es de esperar, el fiujo de neutrinos se desarroUa rapidamente, en la escala de una 
longitud de interaccion de nucleon, y luego permanece estable. Luego de atravesar toda la 
atmosfera (es decir, en el Ifmite X ^ K^) el fiujo de bajas energfas queda dado por 



= Zl^j^^^<t>oNE-^ . (3.18) 



^Debe aclararse que el branching ratio aqui esta inclui'do en la definicion de los momentos de decaimiento, 
en analogi'a con los momentos de produccion/regeneracion que incluyen la multiplicidad de los estados finales. 
Esto coincide con [142], pero difiere de [133]. 
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De un modo similar, el flujo de altas energfas al nivel de la superficie terrestre es 

^g^F^_v7+l ln(AA//Ajv) ea/ ^ ^-(7+1) (o^q,\ 

V — Zj TV at i — j L\m cos V 

donde Q es el angulo zenital asociado a la direccion de incidencia del rayo cosmico primario, 
y donde 

= ^ , (3.20) 

con itlm y tm la masa y la vida media del meson, respectivamente, y donde = 6.4 km es 
una altura caracterfstica de la escala de variaciones de la densidad en la atmosfera. Dado que 
los momentos de produccion/regeneracion {Zjy^, Zj^m J Zmm) son iguales a los que aparecen 
en los flujos de los mesones, deben ser evaluados usando el mdice espectral 7 asociado al flujo 
de nucleones. La ec. ()3.19p asume un perfil de densidad que corresponde a una atmosfera 
isotermica (es decir, p{h) oc exp(— /i//io)), que es apropiado en la estratosfera {h > 11 km), la 
zona en la que ocurre la mayor parte de las interacciones [133, 142]. Ademas, en la ec. (j3.19|) 
se desprecia la curvatura de la tierra, de modo que, en principio, solo sen'a valida para angulos 
zenitales pequefios. Sin embargo, aun para angulos grandes podemos seguir usando la misma 
expresion, con la simple prescripcion de reemplazar 6 —>■ 6*{6,h = 30 km), donde, para una 
dada linea de vision que corresponde a un angulo 6 en la posicion del observador, el angulo 
6* {6, h) es el angulo zenital que sen'a observado en una posicion en la que esta misma Imea de 
vision esta a una altura h con respecto a la superficie [133]. Por ejemplo, para rayos cosmicos 
que inciden en la direccion horizontal {6 = 90°), resulta 6* ~ arcsin(l + h/R(^)~^ ~ 84.5°. 

Para obtener una expresion general que describa los flujos de neutrinos luego de atravesar 
toda la atmosfera, las soluciones asintoticas ()3.18|) - ()3.19|) pueden ser unificadas por medio 
de la funcion de interpolacion 

= ' ^3.21) 

analoga a la funcion de interpolacion dada en la ec. ()3.H1 para uniflcar los flujos de bajas y 
altas energfas asociados a cada componente de rayos cosmicos galacticos. 

Hasta aquf, hemos descripto un esquema de calculo que provee los flujos de nucle- 
ones, mesones y neutrinos producidos por un flujo inicial de nucleones de ley de potencias, 
(pN^E.X = 0) = 0oAf-E~'^, que incide sobre la parte superior de la atmosfera. Ahora ob- 
tendremos los resultados correspondientes al espectro total de rayos cosmicos que Uega a la 
tierra, tomando en cuenta la contribucion de las diferentes especies nucleares que intervienen 
en su composicion. 

Sea (pz el flujo de rayos cosmicos asociado a la componente nuclear de carga Z y mimero 
masico promedio A. Esta componente nuclear aporta un flujo de nucleones dado por 
4>n,z{En) = A'^(f)z{E = AEn). Si el flujo de cada componente de rayos cosmicos de carga Z 
es (pz{E) = (j)oz{E / Eq)~'^^ , el flujo de nucleones correspondiente esta dado por 

ME) = J2^''''<l^oz[^] . (3.22) 
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Consideremos, en particular, el flujo de rayos cosmicos galacticos parametrizado por las 
ecs. p.ip y ()3.2|) . AUi, cada componente de rayos cosmicos de carga Z esta representada por 
dos flujos de ley de potencias 0| y con indices espectrales az y az+Aa, respectivamente. 
Para cada flujo de ley de potencias, podemos calcular el fiujo de nucleones correspondiente 
(dado por la ec. (j3.22p ) y usarlo para calcular el fiujo de neutrinos asociado. Entonces, 
el flujo final de neutrinos producido por esa dada componente de rayos cosmicos puede 
obtenerse mediante la interpolacion de las dos soluciones, usando nuevamente una funcion 
de interpolacion analoga a la ec. ()3.1|) . Finalmente, el fiujo total de neutrinos resulta de sumar 
la contribucion de todas las componentes de rayos cosmicos galacticos, mas la contribucion 
adicional de la componente extragalactica. 

La fig l3.2l muestra los fiujos de neutrinos atmosfericos que corresponden al escenario de- 
pendiente de la rigidez con Act = 2/3, y con las formulas de parametrizacion y la eleccion 
de parametros dadas en la Seccion 3.2.1. La figura muestra la contribucion de diferentes 
componentes y para diversos casos, a saber los fiujos convencionales (horizontales y verti- 
cales), prompt charm/ GBW y prompt charm/CTEQ3. Para el fiujo prompt charm calculado 
mediante el modelo GBW, la figura tambien muestra los fiujos de neutrino totales en las di- 
recciones horizontal y vertical. Finalmente, en el ultimo caso tambien se indican los fiujos 
de neutrinos convencional, charm y total producidos por la componente extragalactica de 
rayos cosmicos. De la ec. ()3.22p . vemos que el espectro de nucleones que resulta del espectro 
de rayos cosmicos es sensible a la composicion asumida; en consecuencia, tambien el fiujo 
de neutrinos atmosfericos depende de la composicion de los rayos cosmicos. En [142, 159] 
se asume, por simplicidad, que la componente de rayos cosmicos dominante consiste solo de 
protones, con un cambio de mdice espectral Aa = 0.3 que da cuenta de la rodilla del espec- 
tro. Para observar el efecto de la composicion de rayos cosmicos en los fiujos de neutrinos, la 
fig J3.3l compara el fiujo atmosferico horizontal total (con la componente prompt charm cal- 
culada con el modelo GBW) segiin se obtiene para diferentes suposiciones sobre el espectro 
de rayos cosmicos, a saber, el escenario dependiente de la rigidez con Aa = 2/3 (que corre- 
sponde a los resultados dados ya en la fig J3.2p . el mismo espectro total pero asumiendo que 
consiste solo de protones, y el espectro de rayos cosmicos usado en [142,159] (donde se asume 
que el espectro esta constitui'do solo por protones, y que tiene un unico cambio de mdice 
Aa = 0.3 para i^.^ = 5 x lO^^eV). En este ultimo caso, el espectro de rayos cosmicos tiene 
una normalizacion menor, en un factor de ~ 2.6, comparado con los anteriores, y esta algo 
suprimido en comparacion con las observaciones actuales. Comparando los resultados para 
el mismo espectro total de rayos cosmicos, el fiujo de neutrinos correspondiente a la composi- 
cion de diferentes especies nucleares resulta menor que el producido por un espectro de rayos 
cosmicos formado solo por protones; ademas, puede verse que este efecto es mas significativo 
a energfas mas altas. La ec. (j3.22j) muestra que la componente nuclear de carga Z y masa A 
da una contribucion al fiujo de nucleones que esta suprimida por un factor ~ A~^-^ a 

bajas energias, mientras que aparece suprimido por un factor ^2-«z-Aa ^ A~^'^ a energias 
altas. Este efecto no fue tratado adecuadamente en [148], donde la supresion fue considerada 
como dada por un factor por otra parte, el impacto del cambio en la composicion en la 
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Figura 3.2: Flujos atmosfericos {u^ + i/^) para un escenario dependiente de la rigidez con 
Aa = 2/3 (que corresponde al espectro de rayos cosmicos de la fig l3.1|l . Se muestran los flujos 
convencionales (horizontales y verticales), prompt charm/GBW y prompt charm/CTEQ3, 
como tambien los flujos de neutrino totales en las direcciones horizontal y vertical para el 
caso prompt charm/GBW. En este ultimo caso, se indican tambien los flujos de neutrinos 
convencional, charm y total producidos por la componente extragalactica de rayos cosmicos. 

region de la rodilla no fue discutido antes en la literatura. La fuerte supresion del flujo de 
neutrinos producido por las componentes mas pesadas implica que las componentes livianas 
(H y He) son responsables de una gran fraccion de los neutrinos aun por encima de su rodilla, 
de modo que el cambio en la pendiente de las componentes individuales aparece reflejado 
tambien en el cambio de pendiente de los flujos de neutrinos. 

Finalmente, tambien es de interes considerar los flujos de muones atmosfericos producidos 
por rayos cosmicos. Por ejemplo, ha sido sugerido recientemente [149] que la observacion 
(mediante telescopios de neutrinos) de muones atmosfericos propagandose hacia abajo podria 
proveer una medida indirecta del flujo de neutrinos atmosfericos prompt, de modo que esto 
podria ser usado para confrontar las predicciones NLO de pQCD. En efecto, debido a la 
cinematica del decaimiento semileptonico de partfculas charm, el flujo de muones prompt 
coincide con el flujo de neutrinos prompt dentro de un ~ 10%, independientemente de la 
energia y del modelo usado para tratar la produccion de charm atmosferico [149]. Por otro 
lado, el flujo de muones atmosfericos convencionales es alrededor de un factor ~ 5 mayor 
que el flujo de neutrinos muon convencionales dentro del rango de energia de interes en 
este trabajo, y exhibe aproximadamente la misma dependencia con la energia [149]. Efectos 
completamente analogos (a los descriptos aqui para los flujos de neutrinos) deberfan esperarse 
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Figura 3.3: Comparacion del flujo total horizontal atmosferico (z/^ + z/^) para prompt 
charm/GBW producido por rayos cosmicos bajo diferentes suposiciones con respecto al es- 
pectro y a la composicion: el escenario dependiente de la rigidez con Aa = 2/3, que considera 
la contribucion de las componentes con carga 1 < Z < 28; el mismo espectro total de rayos 
cosmicos, pero asumiendo que consiste solo de protones; y el espectro de rayos cosmicos 
usado en [142,159] (TIG), que tambien tiene en cuenta solo a los protones, pero ademas con 
una normalizacion mas baja. 

para los muones atmosfericos producidos por los rayos cosmicos que llegan a la Tierra. 

3.3 El canal de cascadas hadronicas y elect romagneticas: 
neutrinos electron 

3.3.1 El sabor de los neutrinos 

Comencemos por reconsiderar las componentes de neutrinos de altas energfas, incluyendo 
ahora las contribuciones relevantes de todos los sabores de neutrino. De aqui en adelante, 
nos referiremos a los flujos de neutrinos y antineutrinos discriminados por sabor. La fig l3.4l 
muestra las principales componentes, y en la Tabla 3.1 se identifican sus caracteristicas mas 
salientes. 

La componente atmosferica convencional de neutrinos electron es menor, en mas de un 
orden de magnitud, que la correspondiente a los neutrinos muon; su principal contribucion 
proviene del decaimiento de kaones. Por otra parte, una componente de neutrinos tau 
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Figura 3.4: Principales componentes del espectro de neutrinos de altas energias, discrimi- 
nadas por sabor. El flujo atmosferico convencional esta promediado sobre el angulo zenital. 
Se incluye la prediccion de Waxman-Bahcall (WB) para el flujo de GRBs [160], y el li'mite 
mas reciente de AMANDA [123]. 



solo puede aparecer, en el flujo atmosferico convencional, a traves de oscilaciones de sabor. 
Debido a que este mecanismo esta muy suprimido a altas energias, aqui podemos ignorar 
esta contribucion (a modo de ilustracion, ver la flg. 7 de [148]). 

Con respecto a la componente atmosferica prompt, las tasas de produccion de neutrinos 
electron y muon en el decaimiento de los mesones charm son muy similares. En cambio, 
el flujo prompt de neutrinos tau es alrededor de un orden de magnitud menor, y aqui sera 
ignorado. Como en la Seccion precedente, la incerteza en la extrapolacion de funciones 
de distribucion partonica para x pequeiio se ilustra mediante el uso de dos funciones de 
estructura diferentes (a saber, CTEQ3 y el modelo GBW). 

El modelo de rayos cosmicos empleado en el calculo de las componentes atmosfericas es, 
como en la Seccion anterior, un escenario dependiente de la rigidez con Aa = 2/3. Mas ade- 
lante, volveremos a considerar los efectos de variar la composicion de los rayos cosmicos. Debe 
notarse que, mientras que los neutrinos atmosfericos prompt son esencialmente isotropicos, 
los neutrinos atmosfericos convencionales tienen una marcada dependencia zenital, con su 
maximo en la direccion del horizonte; aqui presentamos los flujos convencionales promediados 
sobre el hemisferio superior. 

En la flg J3.4l tambien mostramos el Ifmite mas reciente de AMANDA sobre el flujo de 
neutrinos de altas energias, obtenido de sus analisis de cascadas [122,123]. Algunos de los 
calculos previos del flujo de neutrinos atmosfericos prompt son hasta unos 2 ordenes de 
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Tabla 3.1: Breve resumen de las caracteristicas esenciales de las principales componentes del 
flujo de neutrinos. Los espectros correspondientes aparecen en la fig l3.4[ 



Flujo de neutrinos 


Sabores (z/g : z/^ : Ur) 


Dependencia angular 


atmosferico convencional 
atmosferico prompt 
galactico 
extragalactico 


^:1:0 
1:1:^ 
1:1:1 
1:1:1 


maximo hacia el horizonte 
isotropico 

maximo hacia el centro galactico 
isotropico; fuentes puntuales/transitorias 



magnitud mayores que los flujos que consideramos aqui. Aunque probablemente no sean 
realistas, aiin estos grandes flujos serfan consistentes con el Ifmite actual de AMANDA. 
Aqui, en cambio, asumimos el caso menos favorable (pero mas realista) de un flujo prompt 
pequeiio. For otra parte, tambien consideramos un flujo extragalactico pequefio, dado por 
el modelo de Waxman y Bahcall de Gamma Ray Bursts (GRBs) [124,160]. Fara una fuente 
astroffsica de neutrinos generica, se espera que la razon de los flujos de neutrinos producidos 
sea 1:2:0, transformada por las oscilaciones de neutrinos, durante la propagacion, en 
1:1:1 (aunque no puede descartarse que la existencia de nueva fi'sica en el sector del 
neutrino pueda alterar tanto los flujos como los cocientes de sabor [161-166]). Si los flujos 
reales son mayores a los considerados aqui, entonces la tecnica que proponemos es aiin mas 
facil de implementar. 

3.3.2 Los espectros detectados 

En la flg. 13.41 las componentes prompt y extragalactica emergen sobre el flujo de neutrinos 
atmosfericos convencionales, que constituye el fondo dominante a bajas energias, recien a 
energias por encima de 10^ GeV. Fara ser precisos, esto ocurre para el espectro de z/^ y 
el canal de corriente cargada correspondiente, basado en la deteccion de trazas de muones 
de largo alcance. Si el espectro de pudiera aislarse, entonces el cambio del espectro 
podria ocurrir a energias menores en alrededor de un orden de magnitud, donde los flujos 
son mucho mayores y facilitan la deteccion. Nuestra estrategia, entonces, es la de reducir 
el fondo de neutrinos atmosfericos convencionales excluyendo los eventos de z/^ de corriente 
cargada (caracterizados por las trazas leptonicas), y concentrandonos en los eventos de z/g de 
corriente cargada, que generan lluvias (o cascadas) en el detector. Como muestran la flg. 13.41 
y la Tabla 3.1, las sefiales tienen flujos de z/g y z/^ similares, mientras que el fondo debido a 
los neutrinos atmosfericos convencionales tiene un contenido de z/g comparativamente menor. 
Aiin cuando tambien existe una contribucion de z/^ atmosfericos convencionales al canal de 
las cascadas, debido a interacciones mediadas por la corriente neutra, veremos que su impacto 
es reducido. 

En la interaccion de un neutrino con un nucleon, la energi'a del neutrino E^, se comparte 
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Figura 3.5: Flujos diferenciales de cascadas en funcion de la energia visible E^is, corre- 
spondientes a un detector de 1 km^ de tamafio luego de 10 afios de mediciones, utilizando 
solo neutrinos propagandose hacia abajo. Los flujos incidentes sobre el detector son los que 
aparecen en la fig J3.4l Aqui solo se muestran las componentes del espectro de neutrinos con- 
vencionales atmosfericos, con la contribucion por separado de las interacciones de corriente 
cargada (CC) y corriente neutra (NC). La importancia relativa del canal CC crece con 
respecto al canal z/^ NC a mayores energias, debido al decrecimiento de {y). 



entre el quark saliente, al que le corresponde una fraccion y, y el lepton saliente, que se Ueva 
una fraccion l — y. La secciones eficaces diferenciales para interacciones de corriente cargada 
y neutra tienen su maximo en y = 0. En una interaccion Uf, de corriente cargada, el quark 
inicia una Uuvia hadronica de energia ~ yE^, y el electron una Uuvia electro magnetica de 
energia ~ {l—y)E^, de modo que la energia visible total resulta Eyis ~ Ei, (asumiendo que las 
lluvias hadronica y electromagnetica son indistinguibles en el detector). En una interaccion 
de corriente neutra, E^g es menor en un factor (y) ^ 0.4 — 0.3 (decreciente con el aumento de 
la energia) [167]. Ademas, las secciones eficaces totales para la corriente neutra son menores 
que las correspondientes a la corriente cargada, crNc/<^cc — 0.4 [167]. Tomando en cuenta 
que el espectro de neutrinos atmosfericos convencionales decrece muy rapidamente (con un 
indice espectral 7 ~ 3 — 3.7), el fiujo de lluvias originadas en interacciones de corriente neutra 
de z/^ queda suprimido por un factor ~ (y)*-^"^^ x aNc/cfcc, es decir, aproximadamente un 
orden de magnitud en relacion con el fiujo de lluvias originadas en interacciones de corriente 
cargada de u^. 

En consecuencia, en el espectro detectable de lluvias producidas por neutrinos atmosfericos 
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convencionales, las contribuciones de v^. y i/^ son comparables: la diferencia en el flujo (ver 
la fig l3.4|) se compensa con la diferencia en la detectabilidad. Nuestros resultados para los 
neutrinos atmosfericos convencionales se muestran en la fig ]3.5l Como hemos mencionado, 
estamos excluyendo eventos de z/^ de corriente cargada, que pueden ser reconocidos mediante 
trazas de muones de largo alcance. Los espectros que muestra la figura se calcularon convolu- 
cionando los espectros de neutrinos asumidos con la seccion eficaz diferencial (promediada 
entre neutrinos y antineutrinos) [167]. La fig l3.5l muestra que las tecnicas que describimos 
aqui pueden reducir significativamente el fondo de neutrinos atmosfericos convencionales. 

Dado que la contribucion de z/g y a los flujos de neutrinos atmosfericos prompt y 
extragalacticos es comparable, el flujo de las cascadas correspondientes quedara dominado 
por las interacciones de corriente cargada de z/g. Aunque aqui incluimos las interacciones de 
corriente neutra de todos los sabores relevantes, estas contribuciones podn'an ser ignoradas; 
por ejemplo, esto puede verse en la fig ]3.51 comparando las contribuciones de las corrientes 
neutras y cargadas de v^,. 

Hasta aqui, no hemos mencionado las interacciones de v^^ en el caso en que su contribucion 
sea relevante (ver la tabla 3.1). A energias por debajo de -Ej, ~ 5 x 10^ GeV, su interaccion 
de corriente cargada produce solamente Uuvias (con E^is ~ debido a la corta vida media 
del lepton tau. Por encima de esa energia, la longitud de la traza leptonica es suficiente para 
ser identificada y separada de la Uuvia; sin embargo, a energias tan altas, los fiujos son muy 
pequenos, de modo que la separacion de esos eventos es esencialmente irrelevante. Cuando 
Vr esta presente en el flujo, incluimos sus contribuciones de corriente cargada y neutra en 
el espectro de lluvias detectadas. De todos modos, dado que la fraccion de en el flujo de 
neutrinos atmosfericos prompt es muy pequeiia, la identificacion directa de cualquier evento 
de Vr serfa un firme indicador de una serial extragalactica. Una descripcion mas detallada 
de las caracten'sticas de los diferentes tipos de eventos en un detector de neutrinos, y de su 
detectabilidad relativa, se brinda en [161]. 

La fig ]3. 61 muestra los espectros de lluvias detectables correspondientes a los flujos de la 
fig ]3.4l La energia a la que la serial prompt atmosferica, o bien la extragalactica, podria 
emerger del fondo de neutrinos atmosfericos convencionales, es ahora un orden de magnitud 
menor. Esto no hubiera sido evidente si, en cambio, hubieramos presentado el flujo (o el pro- 
ducto del flujo por la seccion eficaz) en funcion de la energia del neutrino. Hooper et al. [150] 
tambien propusieron la medicion de los neutrinos atmosfericos prompt y extragalacticos por 
medio de la deteccion de lluvias. Sin embargo, existen ciertas diferencias importantes entre 
su propuesta y la nuestra. Aqui asumimos que todos los eventos de de corriente cargada 
pueden ser exclmdos, identificando las trazas muonicas de largo alcance; por el contrario, 
en [150] se incluyen todos los eventos en que la interaccion ocurre en el volumen del detector. 
Por otra parte, nosotros calculamos los espectros en terminos de la energia visible, en lugar 
de la energia del neutrino; esto tiene un gran impacto en la reduccion del fondo de neutrinos 
atmosfericos convencionales. Asi, tomar en cuenta estos efectos nos permite considerar un 
caso mas realista, en el que los flujos asumidos para los neutrinos atmosfericos prompt y 
extragalacticos son sustancialmente menores a los estudiados en [150]. 
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Figura 3.6: Espectros diferenciales de cascadas en funcion de la energfa visible E^is, para un 
detector de tamafio km'^ luego de 10 afios de operacion, usando solo neutrinos propagandose 
hacia abajo, y los espectros de la fig l3.41 La definicion de las componentes galacticas, y la 
normalizacion de los flujos, se explican en el texto. 



El canal de las Uuvias no provee informacion especi'fica acerca del sabor de los neutrinos, 
y no puede distinguir entre eventos de corriente cargada y neutra. De todos modos, esto no 
representa una gran desventaja, y esta compensado por la mayor fidelidad que se consigue 
entre la energfa del neutrino y la visible, un aspecto esencial para resolver un posible quiebre 
en el espectro. La resolucion angular es moderada (~ 20°, comparada con la resolucion ^ 1° 
para el canal de z/^ de corriente cargada), pero es perfectamente adecuada para un flujo 
isotropico. 

Una de las desventajas del canal de deteccion de lluvias es que los muones atmosfericos 
pueden producir un fondo significativo cuando, pasando cerca del detector, inician una Uuvia 
a partir de un evento de bremsstrahlung. En efecto, el limite actual de lluvias de AMANDA 
esta dado, precisamente, por el fondo producido por este tipo de eventos [122,123]. Sin 
embargo, dado que este es un efecto de superficie, el tamafio (mucho mayor) de IceCube 
deberfa permftir una reduccion sustancial de este fondo, conservando a la vez un volumen 
fiducial suficientemente grande. Una exclusion similar de la region externa del detector sera 
tambien necesaria, para eliminar los eventos de z/^ de corriente cargada en los que la traza 
del muon escapa a la deteccion (de modo que solo queda registrada la lluvia liadronica). 
Naturalmente, estos cortes reduciran la exposicion que liemos asumido aquf. Para evitar 
los efectos de absorcion en la Tierra, liemos considerado solamente neutrinos propagandose 
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hacia abajo. Los neutrinos que atraviesan un diametro terrestre complete son absorbidos 
para energias alrededor de 4 x 10^ GeV; para distancias menores, la energia de absorcion 
es significativamente mayor (por ejemplo, ver la fig. 2 de [168]). Luego, la exposicion en el 
rango relevante de energias podria incrementarse tomando tambien en cuenta una fraccion 
considerable de los eventos de propagacion hacia arriba. 

Aunque el flujo de neutrinos atmosfcricos convencionales tiene su maximo hacia el hor- 
izontc, aqui ha sido promediado sobrc todo cl hcmisferio superior. Eliminando los eventos 
cercanos al horizonte, podriamos obtener sefiales (no convencionales) identificables a menores 
energias. Un estudio mas complete de la sensibilidad y la capacidad de separar las compo- 
nentes del flujo de neutrinos, usando, por ejemplo, las tecnicas Monte Carlo desarroUadas 
en [122], seria de gran interes. 

3.3.3 El flujo de neutrinos galacticos 

Hasta aquf, hemos discutido los flujos de neutrinos producidos en la atmosfera o en fuentes 
extragalacticas, omitiendo el flujo galactico difuso que resulta del decaimiento de piones y 
muones, producidos en la interaccion de los rayos cosmicos con el medio interestelar (ISM) 
[159, 169-172] (ignoraremos un posible flujo de neutrinos originado en una fuente puntual en 
el centre galactico). 

El flujo galactico difuso puede ser determinade siguiende el mismo precedimiente que ya 
hemos descripto, en la Seccion precedente, para la produccion de los neutrinos atmosfericos. 
Por simplicidad, asumiremos que los rayos cosmicos estan distribuidos homogeneamente 
dentro del disco galactico. Las partfculas ultrarrelativistas interactuan con el gas ambiente 
del ISM, que constituye un plasma no relativista de muy baja densidad, formado funda- 
mentalmente por hidrogeno atomico y molecular. Debemos tener en cuenta que, debido a la 
densidad extremadamente baja del ISM, la interaccion de las partfculas secundarias produci- 
das en las interacciones nucleon-nucleon son completamente despreciables; en consecuencia, 
todas las partfculas secundarias decaen mucho antes de interactuar, independientemente de 
la energfa considerada. Asf, la principal contribucion al flujo de neutrinos difusos producido 
en la galaxia proviene del decaimiento de piones y muones; las contribuciones adicionales 
que provienen de la produccion y el decaimiento de mesones mas pesados pueden ser des- 
preciadas. 

Consideremos un flujo inicial de nucleones dado por 4>n{E,X = 0) = 4>qnE~"' (el flujo 
de neutrinos inducidos por el espectro total de rayos cosmicos puede obtenerse siguiendo 
el procedimiento que se uso en el calculo de los flujos atmosfericos). En este contexto, el 
slant depth atmosferico debe ser reemplazado por la columna de densidad atravesada en el 
ISM a lo largo de la Ifnea de vision; debido a la densidad muy baja del ISM, es ^ X ^ Ajv- 

^En cstc caso, A^r corrcspondc a la longitud dc interaccion para miclconcs propagandosc a travcs del ISM. 
Notese que \n, expresado en g/cm^, resulta muy similar al correspondiente a la propagacion de nucleones 
en el aire. En cambio, la longitud de decaimiento de los mesones A^, siendo lineal en la densidad de masa 
del medio de propagacion, resulta mucho menor en el ISM que en la atmosfera. 
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Entonces, de la ec. ()lll7j] . el flujo de neutrinos muon producido por el decaimiento de piones 
en el ISM es 

X) = ZI^Znm^<PqnE-^ . (3.23) 
An 

Para calcular la contribucion (al flujo de neutrinos muon) que resulta del decaimiento de 
muones, es suficiente con una sencilla estimacion a partir de los resultados ya obtenidos para 
el decaimiento de piones. De la cinematica del decaimiento, la fraccion media de energfa 
(relativa a la partfcula primaria) en el decaimiento vr — >■ /x + i/^ es K^^ = 0.21 para el y 
0.79 para el muon [128]. Analogamente, en el decaimiento /x — e + i^e + ^'^t la fraccion efectiva 
de energia para el z/^ resultante es 0.35 [128], es decir, una fraccion = 0.28 relativa al 
pion original. En consecuencia, se espera que la contribucion del decaimiento de muones 
sea aproximadamente un factor [K^/K-,^)"' mayor que el flujo de neutrinos proveniente del 
decaimiento de piones (por ejemplo, un factor de 2.1 para un mdice 7 = 2.7 en el espectro 
de nucleones, en buen acuerdo con calculos mas detallados [159]). 

Un calculo similar nos permite estimar el flujo de neutrinos electron. En el decaimiento 
ji e+z/e + z/^, la fraccion de energfa efectiva para el z/e es 0.3 [128]. Entonces, el flujo de z/g es 
aproximadamente un factor (0.3/0.35)''' del flujo de producido en el decaimiento de muones 
(es decir, un factor ~ 0.7 para un indice 7 = 2.7 en el espectro de nucleones, nuevamente 
en buen acuerdo con [159]). Finalmente, las oscilaciones de sabor de los neutrinos galacticos 
redistribuyen los flujos producidos de z/^ y z/g uniformemente entre los tres sabores. 

Dado que el flujo galactico es lineal en la columna de densidad atravesada en el ISM 
a lo largo de la linea de vision, su maximo esta en la direccion del centro galactico, y 
es, por lo tanto, anisotropico en coordenadas galacticas. En [170], la densidad del ISM 
esta dada en funcion de coordenadas galacticas, y se asume una densidad minima n = 
0.087 cm"^. Considerando que el halo galactico contiene esta densidad de materia minima, 
las anisotropi'as en la columna de densidad resultan despreciables, excepto en la direccion del 
centro galactico. Asumiendo un halo de 20 kpc de radio y una escala vertical caracteri'stica 
de 2 kpc, la columna de densidad en una direccion tfpica es x-^o cg — 10^^ cm~^, mientras 
que, en la direccion del centro galactico, es xqg — 10^^ cm~^. Consideraremos por separado 
la contribucion de la region del centro galactico {\b\ < 10° y |/| < 10°, que corresponde a un 
angulo solido AficG = 0.12 sr) y la contribucion promedio de todas las demas direcciones en 
el hemisferio superior (que corresponde a AQ^o cc = 6.16 sr). Esta separacion es consistente 
con la resolucion angular esperada para la deteccion de cascadas en IceCube. 

La fig ]3.6l muestra los flujos diferenciales de las cascadas correspondientes a las compo- 
nentes galacticas, comparados con los flujos de las componentes que ya hemos discutido. 
Excluyendo la direccion del centro galactico, el flujo galactico puede despreciarse, en com- 
paracion con las otras componentes del espectro. En la direccion del centro galactico, en 
cambio, el flujo galactico es comparable con las otras componentes. Sin embargo, en la 
fig ]3.61 el flujo en la direccion del centro galactico aparece como si fuera isotropico. Para 
calcular la tasa real de eventos en la direccion del centro galactico, este y los otros flujos 
deben renormalizarse, a fin de tener en cuenta la reduccion en la aceptancia angular; al 



83 



multiplicar todos los flujos por Ar^cc/^^no cg — 0.02, los flujos resultantes son demasiado 
pequeiios para ser observados. 

3.3.4 Los efectos de la composicion de rayos cosmicos 

Las figs l3.7r a)-(b) muestran las tasas integrales de eventos de cascada, correspondientes a las 
componentes relevantes del espectro total de neutrinos de altas energias. En la fig l3.7r a). se 
asume que todas las especies nucleares en el rango 1 < Z < 28 contribuyen a la composicion 
en un escenario dependiente de la rigidez; en la figEKb), se asume que el mismo espectro de 
rayos cosmicos esta compuesto solo por protones. Debe notarse que la composicion asumida 
afecta a los flujos de neutrinos atmosfericos, tanto convencionales como prompt. Como 
en figuras anteriores, se utilizan dos curvas para las predicciones asociadas a los neutrinos 
atmosfericos prompt, con el fin de indicar el rango de incerteza que resulta de usar distintas 
prescripciones para la region de x pequeiio. 

Como ya hemos observado en la Seccion anterior, al considerar una composicion de 
rayos cosmicos en la que interviene una mezcla de especies nucleares, el flujo de neutrinos 
atmosfericos tiene un impacto relativamente menor, y resulta en un fondo significativamente 
reducido para la deteccion de un flujo extragalactico difuso originado en fuentes no resueltas. 

Si bien las tasas integrales que obtenemos son pequenas, el numero de eventos obtenido 
en estas predicciones es realista y razonable. Por otra parte, el Ifmite de AMANDA permite 
la existencia de flujos mayores en alrededor de 2 ordenes de magnitud, que increment arfan 
las tasas de eventos detectados en un factor similar. Aiin antes de descubrir los flujos de 
neutrinos atmosfericos prompt o extragalacticos, la tecnica que describimos aqui permitirfa 
una medicion del flujo atmosferico convencional con gran estadi'stica, algo esencial para 
determinar confiablemente su extrapolacion a mayores energias. 

El quiebre espectral entre el fondo de neutrinos atmosfericos convencionales, que domina 
en la region de bajas energias, y la nueva serial, que domina a energias altas, ocurre para 
una energia menor, en alrededor de un orden de magnitud, que al considerar el canal usual 
de deteccion de trazas leptonicas (comparar las figs J3.7r aVfb) con la fig l3.4p . Una vez que el 
quiebre en el espectro sea observado, diversas caracteristicas permitiran distinguir el flujo de 
neutrinos atmosfericos prompt de un flujo extragalactico; por ejemplo, el mdice del espectro 
observado, la deteccion de una componente de neutrinos tau, o la identificacion de fuentes 
puntuales o transitorias. 
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Figura 3.7: Tasas integrales de eventos de cascada en funcion de la energi'a visible 
para un detector de tamaiio km^ luego de 10 afios de operacion, asumiendo diferentes com- 
posiciones de rayos cosmicos: (a) el escenario dependiente de la rigidez con Aa = 2/3, que 
considera la contribucion de componentes con 1 < Z < 28; (b) el mismo espectro total, pero 
asumiendo que consiste solo de protones. La Ifnea marcada "AMANDA" indica el espectro 
integral resultante asumiendo una ley de potencias E'"^, sin una energia de corte superior, 
normalizado por el li'mite diferencial de AMANDA. 



Capitulo 4 



Fotodesintegracion de rayos cosmicos: 
una explicacion alternativa para la 
rodilla del espectro 

Como una propuesta alternativa para explicar la rodilla del espectro de rayos 
cosmicos, en este Capitulo exploramos el escenario en el que la rodilla resulta 
de la interaccion de los rayos cosmicos con un fondo de fotones opticos y UV 
hlandos que puede estar presente en el entorno de las fuentes. En este esce- 
nario, la rodilla del espectro reBeja la energia umbral asociada a los procesos de 
fotodesintegracion nuclear y de produccion de fotomesones por protones, y la 
composicion exhibe una tendencia hacia componentes mas livianas por encima 
de la rodilla. 

4.1 El escenario de la fotodesintegracion nuclear 

En el sistema en reposo de un rayo cosmico de energia E y masa ArUp, la energia de un 
foton de energia e aparece incrementada por el factor relativista 7 = E/ArripC^. Teniendo en 
cuenta que la fotodesintegracion nuclear tiene tfpicamente una energia umbral en torno de 
algunos MeV (en el sistema en reposo del micleo), resulta que un foton en el rango optico 
(e = 1 — 10 eV) puede desintegrar nucleos con energias E > A x 10^^ eV. En consecuencia, 
la presencia de un fondo de fotones opticos y UV blandos en el entorno de las fuentes podria 
explicar el origen de la rodilla del espectro de rayos cosmicos. Este escenario fue propuesto 
por Hillas [173] en 1979, y luego desarrollado por Karakula y Tkaczyk [174]. 

Aqui presentamos una reconsideracion de este escenario, reevaluando en detalle la propa- 
gacion de nucleos y protones en la region de las fuentes, e incorporando elementos adicionales 
que modifican las predicciones de [174]. En particular, en la reevaluacion de las tasas de 
fotodesintegracion nuclear incorporamos los procesos de emision de mas de un nucleon (que, 
en efecto, constituyen el canal dominante para nucleos pesados), y mejoramos el tratamiento 
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de la energia umbral en la resonancia dipolar gigante asociada a cada nilcleo estable (adop- 
tando los recientes resultados de Stecker y Salamon [175]). Tambicn discutimos el impacto 
del Uamado "mecanismo del neutron" [176,177], que permite que los neutrones escapen de la 
fuente sin perdidas de energia adicionales. Finalmente, introduciendo una energia de corte 
(cutojf) inferior en la distribucion de fotones de ley de potencias, eliminamos la supresion 
abrupta que se observa en el espectro calculado en [174] por cncima dc lO^'' cV. Con la imple- 
mcntacion de estas consideracioncs adicionales, la columna dc dcnsidad dc cncrgia rcqucrida 
para reproducir la rodilla del espectro dc rayos cosmicos en cstc escenario resulta menor, en 
un orden de magnitud, que las estimadas previamente en [174]. 



4.2 La propagacion de rayos cosmicos 

4.2.1 La fuente y el espectro del fondo de radiacion 

En este trabajo, asumiremos que la fuente emite nucleos en el rango de niimero masico 
1 < A < 56. Siendo que, por lo general, existe un linico isotopo estable para cada niicleo 
de masa Arrip en todo el rango que va desde hasta ^^Fe, consideraremos una linica carga 
Z asociada a cada valor de A. Asumiremos que los flujos diferenciales de rayos cosmicos 
emitidos por la fuente son espectros de ley de potencias, 

0O(£;) = $0 E-/^' (0 < i < 55) , (4.1) 

donde i = 56 — A. Las intensidades $^ e indices cspcctralcs j3i se obtienen a partir de las 
observaciones a bajas energias, compiladas en [9,72] (para mas detalles, ver la Seccion 2.2). 

Para el fondo de fotones que rodea a la fuente, estudiaremos el caso de una distribucion 
de Planck, 

P 

""^'^ " 2C(3)(A;BT)3cxp(e/A;BT)-l ' ^^'^^ 

donde es la constante de Boltzmann, T la tcmpcratura absoluta j C la funcion zeta de 
Riemann (tal que C(3) ~ 1.202), y una distribucion de ley de potencias, 

^(e) = P(« - 1) - 7^1) (4-3) 

donde y cm son los cutoffs inferior y superior, respectivamente, y donde a es el mdice 
espectral. En ambas distribuciones, n{e) representa el espectro diferencial correspondiente a 
la densidad de fotones de energia e, mientras que p = f de n(e) es el mimero total de fotones 
por unidad de volumen. Llamando pe a la densidad total de energia, resulta 

Pe ~ 2.631kBTp (4.4) 
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para el espectro de Planck, mientras que 
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(4.5) 



corresponde a la distribucion de ley de potencias. Siguiendo a [174], asumiremos a = 1.3. Sin 
embargo, a diferencia de [174], aqul adoptamos un cutoff inferior que excluye la presencia 
de un abundante numero de fotones de baja energi'a (infrarrojos), que de otro modo sen'an 
dominantes en la fotodesintegracion por encima de 10^'' eV. 

Con respecto a la distribucion espacial del fondo de radiacion, consideraremos que es 
homogenea en toda una region (de dimension lineal L) circundante a la fuente. De cualquier 
modo, los resultados dependeran solo de la columna de densidad de fotones total, integrada 
a lo largo de la trayectoria de los rayos cosmicos desde su produccion en la fuente hasta su 
escape de la region del fondo de fotones. Despreciando efectos de confinamiento en la fuente, 
la propagacion de los rayos cosmicos es esencialmente rectilinea; la columna de densidad 
total atravesada, expresada en terminos de la densidad de energi'a, es PeL. 

4.2.2 Propagacion de nucleos: fotodesintegracion 

El principal mecanismo de perdida de energi'a para nucleos con E < 10^^ eV propagandose 
en un fondo de radiacion optico es el proceso de fotodesintegracion. Un micleo (de numero 
masico A = 56 — i j factor de Lorentz 7) que se propaga a traves de una distribucion de 
fotones n{e) (dada, por ejemplo, por las ecs. ()4.2|1 6 ()4.3j] ) tiene una probabilidad de fision 
con emision de j nucleones dada por 



donde aij es la seccion eficaz de fotodesintegracion correspondiente al proceso, e la energfa 
del foton en el sistema del observador y e' su energi'a en el sistema en reposo del micleo. Para 
calcular Rij, ajustamos aij a los parametros dados en las Tablas 1 y 2 de [178], mientras que 
los umbrales e^^^ j^- fueron tomados de la Tabla 1 de [175]. Para tener en cuenta todos los 
canales de la reaccion, resulta util definir la tasa de emision efectiva, 



A bajas energias {e[f^r,ij < ^' < 30 MeV), la seccion eficaz (Tij{e') esta dominada por la 
resonancia dipolar gigante, y la fotodesintegracion ocurre fundamentalmente mediante la 
emision de uno o dos nucleones. En [174, 178] se asume una linica energfa umbral e^^^^ j^- = 
2 MeV para todos los canales de la reaccion y todos los nucleos. Sin embargo, de acuerdo a 
recientes consideraciones [175], la energfa umbral de las reacciones depende tanto del nucleo 




(4.6) 




(4.7) 
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Figura 4.1: Tasas de emision efectiva versus energi'a por nucleon, para niicleos de ^^Fe 
propagandose en un espectro de fotones de Planck (con ksT = 10 eV). -Rqc/ calculado 
usando la energia umbral fija e^^^^ = 2 MeV, mientras que Ro,ef corresponde a los nuevos 
resultados [175] que involucran energias umbrales especificas. 

como del numero de nucleones emitidos. Incorporando los nuevos resultados, la energia 
umbral se desplaza a valores mayores; en la emision de un linico nucleon, la energia umbral 
tipica resulta ~ 10 MeV, mientras que la correspondiente a la emision de dos nucleones es 
ahora ~ 20 MeV. Aunque la aparicion de la rodilla en este escenario es, en ultima instancia, 
un efecto de umbral, el nuevo calculo, mejorado con la inclusion del cambio en los umbrales 
de reaccion, no modifica significativamente las tasas de emision resultantes en comparacion 
con aquellas obtenidas con la energia umbral fija ej^r.ij — 2 MeV. Esto puede apreciarse en 
la fig J4.11 donde se muestran graficos de Ri^^f/PE en funcion de la energia por nucleon E* 
para un micleo de ^^Fe {i = 0) que se propaga a traves de un espectro de fotones de Planck 
(con ksT = 10 eV). La tasa de emision efectiva Ro,ef (calculada con las energias umbrales 
especificas) exliibe solo un pequeno desplazamiento hacia mayores energias con respecto a 
aquella obtenida por medio de la energia umbral fija (denotada Rq l^). La razon de un efecto 
tan poco significativo es que la fotodesintegracion esta dominada por la region del pico de 
la resonancia gigante, situado tipicamente alrededor de ~ 20 MeV para la emision de un 
nucleon, y en torno de ~ 26 MeV para la emision de dos nucleones, es decir, bien por encima 
de las energias umbrales. 

A mayores energias, crij (e') es aproximadamente plana y la emision de varios nucleones 
adquiere una mayor probabilidad, tornandose dominante para el caso de los niicleos pesados, 
para los cuales la probabilidad de emision de un linico nucleon es de solo 10%. Para estimar 
la relevancia de los procesos de emision de varios nucleones, comparamos la tasa efectiva 
que resulta de las ecs. ()4.6|1 y ()4.7|) con la tasa de emision que desprecia los procesos que 
involucran mas de un nucleon (es decir, la que se emplea en [174]). La fig J4.2r a) muestra 
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Figura 4.2: Tasas de emision efectiva versus energia por nucleoli, para niicleos de ^^Fe, (a) 
para un espectro de Planck (con /c^T = 1.8 eV), y (b) para una distribucion de ley de 
potencias con mdice espectral a = 1.3 y un cutoff superior e*/ = 20 eV. La tasa de emision 
de varios nucleones corresponde a un cutoff inferior = 0.8 eV, mientras que las tasas de 
emision de un linico nucleon son para = 10~^ y 0.8 eV, como se indica. 

Ri,ef/ Pe en funcion de la energia por nucleon E* para un niicleo de ^^Fe {i = 0) que se propaga 
a traves de un espectro de fotones de Planck (con fc^T = 1.8 eV). Resulta evidente que la 
emision de varios nucleones es significativa, incrementando la tasa de emision a energias altas 
por un factor ~ 4 con respecto a los resultados que involucran la emision de un unico nucleon. 
Analogamente, la fig l4.2r b) muestra los resultados correspondientes a un espectro de fotones 
de ley de potencias con un cutoff superior ejv/ = 20 eV. En esta figura tambien puede verse el 
efecto de introducir un cutoff inferior no despreciable en la distribucion. Se observa que, con 
un cutoff despreciable (por ejemplo, = 10~^ eV) la tasa de emision crece monotonamente, 
debido a que la contribucion dominante a la fotodesintegracion viene dada, en este caso, por 
los fotones de baja energia (infrarrojos), que son mas abundantes para mayores valores de 
E*. Estos fotones mantienen el predominio de la resonancia dipolar gigante aiin a energias 
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nucleares altas, y la emision de varies nucleones resulta poco significativa. Sin embargo, si 
consideramos cutoffs inferiores no despreciables (por ejemplo, = 0.8 eV), encontramos 
que, a partir de un cierto valor de E*, la resonancia dipolar gigante queda excluida, y la 
contribucion dominante a las tasas de emision es la del plateau de altas energi'as de cxj. 
En emision satura y adquiere tambien un plateau. 

Comparando las dos curvas correspondientes a la tasa de emision con cutoffs no despreciables, 
observamos aqui nuevamente que los procesos de emision de varios nucleones conducen a un 
incremento significativo (por el mismo factor ~ 4 que ya habfamos encontrado antes) en la 
tasa de emision. De esta manera, vemos que asumir un cutoff inferior en el espectro de fotones 
de ley de potencias debe conducir a resultados similares a los que resultan de considerar un 
espectro de Planck. En este caso, siendo que las tasas de emision de varios nucleones son 
considerablemente mayores (que las correspondientes a la emision de un linico nucleon), es 
de esperar que se requiera una densidad de fotones menor a la estimada previamente en [174] 
para producir efectos similares. 

Debe notarse que, en las figs l4.1l v l4.2t las tasas de emision aparecen graficadas en funcion 
de la energfa por nucleon E*, en lugar de la energia por nucleo E. Esto es debido a que 
E* permanece constante durante la fotodesintegracion y, por otra parte, porque la forma 
de las secciones eficaces de los diferentes micleos son similares, con sus respectivos maximos 
ubicados aproximadamente en los mismos valores de E*. 

Dado que el proceso dominante que afecta el transporte de los micleos a traves del fondo 
de radiacion es la fotodesintegracion nuclear, asumiremos un modelo simplificado en el que 
los micleos se propagan en forma rectilmea, despreciando otros efectos (tales como difusion, 
conveccion, fragmentacion nuclear (spallation), reaceleracion, etc.). Entonces, las ecuaciones 
de transporte que describen la propagacion de los micleos resultan 



donde (f)i{E*,x) es el fiujo diferencial correspondiente a un nucleo de masa A = 56 — i 
{0 < i < 54) con una energfa por nucleon E* a la distancia de propagacion x (medida desde la 
fuente), Rij{E*) es la tasa de emision, y 6ij es la delta de Kronecker. Con las simplificaciones 
asumidas, este sistema de 55 ecuaciones diferenciales acopladas admite soluciones analiticas 
exactas; la solucion (es decir, los flujos diferenciales de todas las especies nucleares a la 
distancia de propagacion L) viene dada por 




ME*,X)J2R^AE*) + il-So^)J2Ri^-J)AE*)<P^~AE*,x) , (4.8) 



dx 





(4.9) 



donde 




(4.10) 



i=o 
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% = TTTT^^ „ rm.^ ' P^^a z > J • (4.11) 



YlUi b(i^k)j{E*)R(^i_k)k{E* 

Puede notarse que, dado que las tasas de emision Rij son lineales en la densidad de 
energi'a de fotones pe, los flujos resultantes dependen solo de la columna de densidad pcL, 
como era de esperar. Notese ademas que, dado que las ecuaciones de transporte dependen 
de la energi'a por nucleon E* , los espectros finales de las diferentes componentes nucleares 
de rayos cosmicos requieren recalcular las soluciones (ecs. fl4.9|) - ()4.1H) ) en terminos de la 
energi'a por niicleo E. 



4.2.3 Propagacion de protones: produccion de fotomesones 

El principal mecanismo de perdida de energi'a que debe tenerse en cuenta en la propagacion 
de los protones es el proceso de produccion de fotomesones. Deberemos distinguir entre el 
flujo de protones primarios producidos en la fuente (0^=55) y los nucleones liberados como 
producto de la fotodesintegracion nuclear {(^^ y para los protones secundarios y neutrones, 
respectivamente). Sea (f)f{E,x) el flujo diferencial de nucleones de energfa E producidos 
dentro de una distancia x desde la fuente por la fotodesintegracion de nucleos de masa 
A = 56 — z y energi'a por nucleon E* = E/A. Entonces, la ecuacion que gobierna la evolucion 
de los fiujos es 

^^l^l^ = <P,{E*,x)Y,JR^AEl , para0<^<54, (4.12) 

con la condicion inicial (j)f^{E,x = 0) = 0. La solucion es 

/ J2k>i kRik{E 



r- 



V ^k>\Rjk{E*) 



1-exp -^i?,fc(E*)L 



k>\ 



(4.13) 



Los fiujos totales de protones secundarios y neutrones se obtienen mediante las sumas 
pesadas sobre el fndice i: 



54 

56 



ni?,L) = 5^-^0f(i?,L) (4.14) 



i=0 

y 

^E, ^) = E (1 - 0f (i?, ^) , (4.15) 

i=0 ^ ^ 

donde Zj es la carga del micleo i— esimo. 

En presencia de un campo magnetico en la region de la fuente, las partfculas cargadas 
quedan retenidas allf por algiin tiempo, mientras que las no cargadas pueden escapar mas 
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rapidamente. Este fenomeno, el "mecanismo del neutron" , ha sido propuesto en el contexto 
del estudio de mecanismos de produccion de rayos cosmicos ultra-energeticos [176,177]; aqui 
discutiremos la relevancia de este mecanismo en el escenario que estamos estudiando. Bajo 
estas consideraciones, asumiremos que el flujo de neutrones 0"" se suma al flujo saliente total 
de rayos cosmicos, despreciando cualquier perdida de energi'a adicional en su propagacion a 
traves del fondo de fotones que rodea a la fuente. En cambio, el flujo de protones secundarios 
(p^^, combinado con el flujo inicial de protones primarios 0^=55, experimenta perdidas de 
energia debido al proceso de produccion de fotopiones, a traves de las reacciones 

p + 7-^p + 7r° (4.16) 

y 

p + 7 ^ n + 7r+ . (4.17) 
La probabilidad de interaccion viene dada por 

1 l"^ n(e) f'^"'^ 
g{E) = — de^ de'e'a,,{e')K{e') , (4.18) 

donde, analogamente a la ec. ()4.6|1 . 7 = E/rripC^ es el factor de Lorentz del proton, a^p 
es la seccion eflcaz correspondiente a la interaccion 7^ (obtenida de ajustes a los datos 
experiment ales compilados en [179]), K es el coeflciente de inelasticidad (deflnido como la 
perdida de energia relativa promedio del proton), e es la energia del foton en el sistema del 
observador, y e' su energia en el sistema en reposo del proton. La energia umbral para la 
produccion de fotopiones es e^^^^ = 145 MeV. 

Entonces, Uamando 0^(-E, x) y 0*"(i?, x) a los flujos diferenciales de protones y neutrones 
de energia E (producidos por las reacciones ()4.16|) y ()4.17p . respectivamente) luego de atrav- 
esar una distancia x desde la fuente, las ecuaciones de propagacion correspondientes quedan 
dadas por 

''^'^•^) -^,iEmE,.)^l^JJL^]^(Ji^^ ,4.19) 



dx V ' / 2(1- Ky \1- K J ^ \l- K 

y (despreciando perdidas de energia en el flujo de neutrones) 



dx 2{1-Ky\l-Kj^ \1-K 

donde los factores ^ resultan de la probabilidad relativa de ocurrencia de las reacciones (j4.16|) 
y (j4.17p . Al coeflciente de inelasticidad se le dio el valor tipico K = 0.3, independiente de 
la energia del proton [174]. Las condiciones iniciales correspondientes a las ecuaciones de 
propagacion, ()4.19p y ()4.20p . son 

cf)P(E,x = 0) = <PU,iE) + r%E,L) (4.21) 
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y 



\E,x = 0) = 



(4.22) 



respectivamente. En la primera se asume, por simplicidad, que los protones secundarios 
producidos por la fotodesintegracion nuclear se suman directamente al espectro inicial de 
protones. Esta suposicion esta justificada debido a que, en los casos en que la produccion de 
fotopiones es relevante, la fotodesintegracion es muy eficiente y los nucleos son desintegrados 
rapidamente. La solucion a la ec. ()4.19p es 



\E,L) = J2<P^{E,L) 



1=0 



donde 

mientras que, para / > 

E 



bl{E,L) = r{E,0)exp{-g{E)L) , 



{E,L) 



(l-K) 



{2(1 -K)) 



E 



[i-Ky 



exp -g 



L 



j=l ^ ^ ' n=0 1 Im-o g y 

Para la componente de neutrones, la solucion a la ec. ()4.2()|1 resulta 



9 



1=0 



donde 



r-{E,x) = Y,<priE 



X) 



2(1 -K) 



1 - exp -g 



E 



l-K 



L 



mientras que, para que I > 0, 



r{E,L) 



£^ n I 

11^ 



(2(1 -i^)) 



i+i 



j=0 



E 



I 

E- 

n=o g 



1 - exp ( ( ) L 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



9 I (i_/^)m+i J 9 y(i-Kyi+i J 

La linealidad de la tasa de produccion de fotopiones 9{E) con respecto a la densidad de 
energia de fotones pe (ver la ec. (j4.18p ) hace que los flujos de protones y neutrones dependan 
solo de la columna de densidad PeL, tal como ocurre en el caso de los nucleos. 
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4.3 Result ados y comparacion con las observaciones 

Hasta aqui, hemos planteado las ecuaciones que describen la propagacion de los rayos 
cosmicos en el entorno de las fuentes, y encontramos su solucion exacta: ecs. ()4.9j) - ()4.1H) 
para los nucleos, ecs. ()4.13j) . ()4.15|) y ()4.2(ij) - ()4.28j) para los neutrones (que, eventualmente, 
decaen en protones fuera de la region de la fuente), y ecs. ()4.23j) - ()4.25j) para los protones (que 
sufren perdidas de energia debido a la produccion de fotopiones), siendo la suma de todas 
estas contribuciones el flujo total de rayos cosmicos (ptotai- Ahora podemos examinar como 
depende 4>totai de la distribucion de fotones asumida para el entorno de la fuente, es decir, 
como es su dependencia en terminos de la columna de densidad peL y de los parametros 
involucrados en las distribuciones espectrales de fotones (ecs. (j4.2|) y (|4.3p ). 

La fig J4.3r a) muestra el flujo diferencial total de rayos cosmicos (en funcion de la energfa 
por particula) para un espectro de Planck con fc^T = 1.8 eV y diferentes valores de la colum- 
na de densidad peL. Analogamente, la fig l4.3r b) muestra los resultados correspondientes a 
un espectro de Planck con fc^T = 10 eV. En las figuras tambien aparecen los espectros 
observados mediante diferentes experimentos (CASABLANCA, DICE, Tibet, PROTON y 
AKENO) [73,74,180]. Segun se observa claramente en las figuras, al aumentar la temperatu- 
ra la rodilla se desplaza liacia menores energfas, mientras que el efecto de incrementar peL es 
el de intensificar la supresion del flujo por encima de la rodilla. Este comportamiento resulta 
del hecho de que la fotodesintegracion nuclear es el proceso responsable de la formacion de 
la rodilla (debido a que los procesos de produccion de fotopiones aparecen recien a mayores 
energias). Entonces, si la temperatura de la distribucion (y, por lo tanto, la energfa media 
de los fotones) es mayor, los nucleos comienzan a desintegrarse a energfas mas bajas. Por 
otra parte, la propagacion de rayos cosmicos no puede verse afectada por peL debajo de la 
rodilla, a energias en las que no hay fotodesintegracion; en cambio, por encima de la rodilla, 
las tasas de desintegracion son lineales en pe- En consecuencia, debe esperarse que peL sea 
la cantidad que regule la supresion del espectro en la region de la rodilla, tal como se observa 
en las figuras. 

El mejor ajuste a las observaciones para el espectro de fotones de Planck, reportado 
en [174], ocurre para fc^T = 1.8 eV y peL = 2.25 x 10^^ eV/cm^; para comparar, la 
figES^a) muestra nuestros resultados para la misma temperatura. Como puede observarse, 
los nuevos resultados son compatibles con los datos experiment ales para valores de columna 
de densidad alrededor de peL ~ 5 x 10^^ eV/cm^, es decir, menores en un factor ~ 50 que el 
resultado previo. La figES^b) muestra que los resultados correspondientes a ksT = 10 eV 
estan en buen acuerdo con los datos experimentales para columnas de densidad en el rango 
PeL = 5 X 10^7 - 2 X lO^^ eV/cm^. 

Las figs l4.4f a) y (b) muestran la composicion masica media, (InA), correspondiente a 
espectros de Planck con fc^T = 1.8 eV y ksT = 10 eV, respectivamente, en comparacion con 
los datos de diferentes experimentos [73,74,181]. Las soluciones con los mayores valores de la 
columna de densidad peL predicen que los rayos cosmicos arriba de la rodilla estan mayor- 
mente constitufdos por protones, en desacuerdo con los datos experimentales. La contribu- 
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Figura 4.3: Flujo total de rayos cosmicos para espectros de fotones dc Planck con (a) ksT — 
1.8 eV y (b) ksT = 10 eV, y diferentes valores para la columna de densidad PeL, como se 
indica. Los resultados se comparan con los datos observados por diferentes experimentos. 

cion no despreciable de componentes mas pesadas a energias por encima de la rodilla impone 
una cota superior sobre la columna de densidad alrededor de {pEL)max ~ 2 x 10^^ eV/cm^. 
Este resultado hace ver nuevamente que, en este escenario, los procesos de emision de varios 
nucleones juegan un rol muy significativo en la propagacion de los rayos cosmicos en la region 
de la fuente. 

Podemos observar tambien que, aunque las soluciones con columnas de densidad rela- 
tivamente bajas {peL < {pEL)max) estan en un acuerdo razonable con algunos conjuntos 
de datos experimentales dentro de las barras de error, los resultados que se obtienen para 
la composicion arriba de la rodilla parecen estar demasiado suprimidos. Sin embargo, debe 
tenerse en cuenta que aqui hemos considerado solo una fuente con una columna de densidad 
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Figura 4.4: Composicion masica de rayos cosmicos para cspcctros dc fotoncs dc Planck con 
(a) ksT = 1.8 eV y (b) ksT = 10 eV, y diferentes valores para la columna de densidad PeL, 
como se indica. Los resultados se comparan con las observaciones de diferentes experimentos. 

fija, mientras que en una situacion realista se tendrian muchas fuentes, cada una, en prin- 
cipio, con un entorno de radiacion diferente. Evitando calculos detallados, consideremos el 
efecto que resultan'a de agregar una fuente adicional. Asi, si esta segunda fuente tuviera una 
columna de densidad peL algo menor, daria lugar a un espectro menos suprimido; entonces, 
mientras la primera fuente seria responsable de la formacion de la rodilla, la segunda man- 
tendria una componente nuclear mas alia de la rodilla. Por otra parte, en nuestros calculos 
hemos asumido que todos los niicleos rccorren la misma distancia L en la region que con- 
tiene la radiacion; es decir, no hemos tcnido en cuenta que, para cualquier dada cncrgia, 
los campos magneticos presentes en el entorno de la fuente retendrian mas efectivamente a 
los rayos cosmicos mas pesados. Tomando esto en cuenta (por ejemplo, asumiendo que un 
micleo recorre una distancia proporcional a su carga) encontrariamos que, mientras la desin- 
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Figura 4.5: Flujo total de rayos cosmicos para espectros de fotones de ley de potencias 
con a = 1.3, em = 20 eV, peL = 10^* eV y diferentes valores para el cutoff inferior de la 
distribucion de fotones, como se indica. Los resultados se comparan con las observaciones 
de diferentes experimentos. 

tegracion podria ser muy efectiva para niicleos pesados (y seria, de heclio, la responsable 
de la formacion de la rodilla), seria mucho menos efectiva para los nucleos mas livianos (ya 
que la columna de densidad para ellos seria mucho menor) y entonces podrian sobrevivir sin 
desintegrarse hasta energias mas alia de la rodilla. 

De este modo, vemos que este escenario puede reproducir el comportamiento obser- 
vado por algunos experimentos, en el que la composicion cambia liacia el predominio de 
componentes mas livianas por encima de la rodilla. Sin embargo, como hemos comenta- 
do anteriormente, existe una marcada discrepancia en el conjunto de datos experimentales 
disponibles. Los resultados de los experimentos mas recientes (KASCADE [77] y EAS- 
TOP/MACRO [78]) parecen concordar en que la composicion se torna mas pesada por 
encima de la rodilla, favoreciendo asi a los escenarios dependientes de la rigidez (como, por 
ejemplo, el escenario de la difusion turbulenta y drift estudiado en el Capitulo 2), aunque 
otros experimentos muestran la tendencia contraria (por ejemplo, DICE [73]). 

Los resultados para un espectro de fotones de ley de potencias (con un cutoff inferior 
en el rango infrarrojo) no presentan mayores diferencias con los obtenidos para espectros 
de Planck. En lo que sigue, fijaremos el cutoff superior, tomando el valor cm = 20 eV. La 
fig ]4.5l muestra el fiujo total de rayos cosmicos (en funcion de la energfa por partfcula) para 
distribuciones de ley de potencias con una columna de densidad peL = 10^^ eV/cm^ y varios 
valores para el cutoff inferior del espectro; los espectros obtenidos reproducen adecuadamente 
la rodilla. Analizando la dependencia de los resultados con el cutoff inferior del espectro de 
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fotones, vemos que, mas alia de pequenas diferencias en la normalizacion, los valores mas 
bajos de producen flujos mas suprimidos en la region de energias altas. De la expresion 
para la tasa de emision de la fotodesintegracion (ec. ()4.6|) ). puede verse que el rol de 
es el de determinar un valor para la energia de los rayos cosmicos E, tal que, debajo de E, 
domina la resonancia dipolar gigante, mientras que, por encima de E, prevalece el plateau 
de a propio del regimen de altas energias. A medida que bajamos e^, disponemos de un 
numero mayor de fotones de baja energia, de modo que la supresion del flujo impuesta por 
la resonancia gigante se extiende a energias de particula mayores. 

Este analisis permite explicar apropiadamente los resultados de la fig J4.51 como asi tam- 
bien los que muestra la fig l4.2f b). Los resultados obtenidos previamente para espectros de 
fotones de ley de potencias [174] exhiben una supresion del flujo extremadamente abrupta 
a altas energias (por encima de E = 10^^ eV) que parecen estar en desacuerdo con las ob- 
servaciones del espectro y, especialmente, de la composicion de los rayos cosmicos. Como 
vemos aqui, esto se atribuye al valor despreciable que se adopto alii para e^; el uso de un 
cutoff no despreciable evita la supresion abrupta del flujo de rayos cosmicos y provee un 
mejor acuerdo con las observaciones en esta region de energias. Por otra parte, debe no- 
tarse que una solucion que exhiba un flujo de rayos cosmicos excesivo a energias mayores 
que 10^^ eV no es problematica, ya que es de esperar que a energias tan altas el escape de 
rayos cosmicos de la galaxia comience a ser un mecanismo efectivo de supresion del flujo. 
Adicionalmente, la aceleracion de los rayos cosmicos en la fuente podria ser menos efectiva 
a altas energias, contribuyendo tambien a suprimir el flujo. En consecuencia, teniendo en 
cuenta la existencia de estos mecanismos de supresion adicionales, no seria razonable esperar 
una explicacion del espectro de rayos cosmicos solo mediante procesos de fotodesintegracion 
mas alia de ~ 10^"^ eV. 

La fig ]4. 61 muestra las componentes de protones y neutrones {(ff + y 0" + 0*", respec- 
tivamente) y el flujo total de rayos cosmicos (ptotai, para el espectro de fotones de Planck 
con ksT = 1.8 eV y peL = 2 x 10^^ eV/cm^. Con el fin de estimar el posible impacto del 
mecanismo del neutron, tambien aparece graficado el flujo total de rayos cosmicos, 4>totaiy 
se obtiene ignorando el escape de neutrones. En la figura se observa que el mecanismo del 
neutron incrementa el flujo total de rayos cosmicos en ~ 100% para ii^ > 5 x 10^^ eV. Si, en 
cambio, hubieramos adoptado una columna de densidad peL = 10^^ eV/cm^, el incremento 
de flujo correspondiente hubiera sido de solo ~ 10%. Sin embargo, debe notarse que estos re- 
sultados corresponden a columnas de densidad muy grandes (de hecho, incompatibles con la 
composicion de rayos cosmicos observada); los efectos que resultan del mecanismo de escape 
del neutron resultan, en cambio, despreciables para menores valores de PeL, como tambien 
las perdidas de energia de los protones por produccion de fotopiones. 

Finalmente, estimaremos las condiciones que deben imponerse sobre las fuentes para 
que puedan producir las densidades de fotones que requiere este escenario. Una fuente 
isotropica con una densidad de energia de fotones pe,o en su superflcie (con radio vq) tendra 
a su alrededor una densidad de energia de fotones dada por p_B(r) = PEflijo/rY . Esto 
signiflca que los rayos cosmicos que abandonan la fuente radialmente atraviesan tipicamente 
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Figura 4.6: Flujo de protones (0^ + 0^^), neutrones + y total {(ptotai) para un espectro 
dc fotones de Planck con ksT = 1.8 eV y peL = 2 x 10^^ eV/cm^. Tambien se muestra el 
flujo total de rayos cosmicos {(p^otai) obtiene ignorando el escape de neutrones. 

una columna de densidad de fotones (per) ~ Pe,qTq. For otra parte, si los rayos cosmicos 
estan confinados magneticamente cerca de la fuente, el camino que recorreran a traves de 
la region con la radiacion sera mayor, es decir, (peL) — X{pEr), con A ^ 1. Si, con el 
fin de estimar las maximas densidades que es posible producir, asumimos que la fuente 
radi'a con una luminosidad L = Aur^pEflC menor que el Kmitc de Eddington (es decir, 
L < 

LEdd = ^'KGMrripc/ gt — 1-3 x 10^^ M/M© erg/s, donde M es la masa de la fuente 
y aT la seccion eficaz de Thomson) , entonces el radio de la fuente debe satisfacer 

ro < ^ ~ — — 2 X lO'^m , 4.29) 

de modo que la fuente debe ser compacta (como un pulsar) y, en general, mas pequena que el 

tamaiio ti'pico dc una cnana blanca. Sin embargo, debe tenerse en cuenta que, en situaciones 
no estacionarias (como, por cjcmplo, cn explosiones de supernova), se pueden alcanzar lumi- 
nosidades mayores que LEdd- Otra consideracion de interes es que, si la fuente produce un 
espectro de cuerpo negro con temperatura T y esta radiando en el Kmite de Eddington, la 
relacion ^TrrlaT"^ = LEdd implica que T ~ 1.8 x 10^ eV(M/M0)i/4(lO km/ro)^/^ de modo 
que, para que la energi'a de los fotones este en la region 6ptica/UV, los radios algo mayores 
que aquellos de las estrellas de neutrones estarian favorecidos (excepto que la luminosidad 
sea menor que LEdd)- 



Capitulo 5 
Conclusiones 



En esta Tesis Doctoral estudiamos diversos aspectos relacionados con la propagacion de rayos 
cosmicos galacticos de muy alta cnergia, con particular enfasis cn la region que comprende 
las dos rodillas del espectro (es decir, en el rango E ~ 10^^ — 10^*^ eV). 

En el Capitulo 2, consideramos en detalle la difusion y los efectos de drift en la propa- 
gacion de rayos cosmicos en configuraciones plausibles de campos magneticos galacticos, 
adoptando coeficientes de difusion apropiados para los niveles de alta turbulencia presentes 
en la galaxia, y asumiendo un espectro de Kolmogorov de fiuctuaciones. En cstc csccnario, 
cada componente nuclear de rayos cosmicos adquiere una rodilla (con un cambio de mdicc 
espectral Aa ~ 2/3) que ocurre a una energia proporcional a la carga correspondiente; de 
este modo, la rodilla del espectro resulta del escape de las componentes livianas (esencial- 
mente, protones y particulas alfa), mientras que el escape de las componentes mas pesadas 
(basicamente, los micleos de hierro) forma la segunda rodilla. Encontramos que el escenario 
de la difusion turbulenta y drift describe adecuadamcnte el espectro observado desde energias 
por debajo de la primer rodilla hasta las correspondientes a la region del tobillo, en donde 
tiene lugar un cambio hacia el predominio de la componente extragalactica de rayos cosmicos. 
La composicion que predice este escenario tiende a ser mas pesada a mayores energias a traves 
de la rodilla, y esta en excelente acucrdo con los datos de algunos experimentos, como, por 
ejemplo, las mas recientes observaciones de KASCADE y EAS-TOP/MACRO. 

Por otro lado, la amplitud de anisotropia experimenta un crossover desde el regimen 
dominado por la difusion transversal al regimen dominado por los drifts, que aparece alrede- 
dor de la region de la rodilla. En efecto, esta es una caracteristica particular que diferencia 
claramente las prcdiccioncs de cstc escenario de las de otras propuestas. Ademas, exam- 
inando la fase de anisotropia del primer armonico en ascension recta, encontramos que el 
maximo de la anisotropia podria estar oricntado hacia la direccion del centro galactico para 
energias cercanas a 10^^ eV, tal como fuc observado por AGASA y SUGAR. Las predicciones 
de este escenario resultan tambien de potencial interes para las futuras observaciones en el 
hemisferio sur con el observatorio AUGER. 

Finalmente, complementamos este estudio con la determinacion numerica, mediante sim- 
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ulaciones extensivas Monte Carlo, de los coeficientes de difusion paralela, transversal y anti- 
simetrica que describen la propagacion de rayos cosmicos en condiciones de alta turbulencia. 
Extcndicndo investigaciones previas, analizamos diferentes tipos de espectro de turbulencia 
y proponemos parametrizaciones sencillas que describen los resultados obtenidos en cada 
caso. 

En el Capitulo 3, calculamos el flujo difuso de altas energias de los neutrinos producidos 

per rayos cosmicos que inciden sobre la atmosfcra, o bien en su propagacion en el medio 
interestelar de la galaxia. Los resultados dependen significativamente del espectro y de la 
composicion asumidos. En particular, considerando que la composicion se torna mas pesada 
por encima de la rodilla (es decir, como resulta de considerar escenarios que explican la rodilla 
como debida a un cfecto dependiente de la rigidez) el fondo de neutrinos atmosfcricos se 
reduce significativamente, facilitando la biisqueda de fuentes astrofisicas. Lucgo dc invcstigar 
este fcnomcno en el canal usual de trazas leptonicas, asociado a la deteccion de neutrinos 
muon, proponemos un nuevo metodo para aislar el flujo de neutrinos atmosfericos prompt 
a traves del canal de las cascadas hadronica y electromagnetica, asociado a la deteccion 
de neutrinos electron. Este cstudio mucstra la importancia de resolver la composicion para 
estimar apropiadamente el fondo difuso de neutrinos de altas energias producido por los rayos 
cosmicos, y es particularmente relevante en la posible deteccion de neutrinos extraterrestres 
en experimentos tales como AMANDA y ICECUBE. 

En el Capitulo 4, consideramos en detalle la fotodesintegracion nuclear de rayos cosmicos, 
en el csccnario en que su intcraccion con fotones opticos y UV blandos en el entorno de la 
fuente da origen a la supresion del flujo y al cambio en la composicion en la region de la 
rodilla. El tratamiento mas preciso de los procesos nucleares, particularmente la inclusion 
de procesos que involucran la emision de varios nucleones, implica que las columnas de den- 
sidad de fotones requeridas sean menores (por un orden de magnitud) a las calculadas en 
investigaciones previas. Como una propucsta alternativa para explicar la rodilla del espec- 
tro, este modelo predice que la composicion se torna mas liviana a partir de ~ 10^^ eV; 
esto constituye una diferencia sustancial con los escenarios dependientes de la rigidez, que 
predicen la tendencia opuesta, y podra ponersc a prucba con las mediciones mas precisas 
que se obtendran proximamente en KASCADE-Grande. 

En conclusion, enfocando nuestra atencion hacia la propagacion de rayos cosmicos galacticos 
de energia, hemos elaborado propuestas que explican las rodillas del espectro, la 

composicion y las anisotropias; encontramos, ademas, que una solucion a cstos problcmas 
podria tambicn contribuir signiflcativamente al conocimiento de la distribucion de las fuentes 
de rayos cosmicos en la galaxia, de la region en que los rayos cosmicos se aceleran, de la es- 
tructura del campo magnetico galactico, y de los flujos de neutrinos de muy alta energia que 
podran detectarse en los nuevos observatorios. 



Apendice A 

El tensor de difusion 



El caracter difusivo del transporte de rayos cosmicos en la galaxia resulta de considerar 
al campo random como un mecanismo efectivo de colision. Las componentes diagonales 
del tensor de difusion corresponden a la "difusion browniana", es decir, a la aparicion de 
un flujo de parti'culas en respuesta a un gradiente de concentracion, ambos en la misma 
direccion (J = —DVN). Como el campo regular rompe la isotropia, debe distinguirse la 
direccion transversal al campo regular de la direccion paralela al mismo; considerando que 
el campo magnetico esta localmente dirigido en la direccion z, las componentes diagonales 
quedan definidas como D^x = Dyy = D±, D^z = D\\- 

El origen fisico del drift macroscopico queda ilustrado por la fig lA.ll que muestra las 
orbitas helicoidales que corresponden a parti'culas de carga positiva, proyectadas sobre un 
piano perpendicular al campo magnetico regular (por simplicidad, aqui ignoramos la distor- 
sion que introduce la componente random del campo magnetico). El resultado de promediar 
el flujo sobre una distribucion de parti'culas es la aparicion de una corriente macroscopica 
perpendicular al gradiente de concentracion; debe notarse que este efecto resulta del com- 
portamiento colectivo de la distribucion de parti'culas, y esta ausente a nivel microscopico 
(es decir, considerando trayectorias de particulas individuales) . Por otra parte, es inmediato 
observar que la componente de drift es antisimetrica: D^y = —Dy^ = D^. Finalmente, 
siendo que la fuerza de Lorentz es perpendicular al campo magnetico que la produce, este 
efecto esta ausente en los pianos x — z, y — z. 

En terminos de las componentes de b = Bo/|Bo|, el vector unitario en la direccion del 
campo magnetico regular, el tensor de difusion queda dado por la ec. ()2.9|) . 

Dij = (Dii - D±) bibj + D±5ij + DAeijkbk ■ 
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Figura A.l: Esquema que ilustra el origen fisico del drift macroscopico (de [13]). 



Apendice B 

El formalismo de Kubo 



La difusion de particulas cargadas en campos magneticos random puede calcularse a partir 
de la descorrclacion entre difcrcntes componentes de la velocidad de una particula (una 
vez promediada sobre una distribucion isotropica de muchas particulas, y sobre muchas 
configuracioncs diferentes de campo random), en lo que se conoce como el formalismo de 
Kubo [66,69,99]. 

La funcion de correlacion de velocidades puede definirse como 

Rij{t) = {viito)vj{to + t)) (t > 0) , (B.l) 

donde (...) denota el promedio sobre un ensemble de configuraciones de campo random. 
Esta definicion asume la ergodicidad del sistema (es decir, que Rij es indcpcndiente de to). 
Asumiremos que Rij{t) suficientemente rapido en el Kmite i — > oo, de modo que 

poo 

/ dt Rij{t) < oo , (B.2) 
Jo 

y definiremos Tef como algiin valor de t suficientemente grande, tal que 

/ ' dt Rij{t) > / dt Rij{t) . (B.3) 



Calculemos el Kmite 



Siendo que 



(Axi(T)AxJT)) 



{Axi{T)Axj{T)) = [ dt [ dt'{vi{t)vj{t')) 
Jo Jo 

dt'Rji{t')+ / dt'Rij{t') 
Jo 



(B.5) 
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resulta que, si i 7^ j y Rijit) = —Rji(t) = R(t), entonces 

r-T rt 



(Ax,(T)Ax,(T)) 



T/2 i>T-t 

dt / dt'R{t') -I dt I dt'R{t') 

Jt Jt/2 JT-t 



Si Rij{t) = Rjiit) = R{t), entonces 

(Ax.(T)Aa;,(T)) 



(B.6) 



dt+ dt+ dt 

JT.f JT-T,f 

Asumiendo T ^ Tg/, el primer termino queda 
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dtR{t) + / dtR{t) 
Jo 



(B.7) 



~ 2 Te/ / dt R{t) 




(B.8) 



analogamente, se llega al mismo resultado para el tercer termino. La contribucion del segundo 
termino es 



T-T, 



dt 



ef 



dt'R{t' 



T-t 



dt'R{t' 



2T I ^ dt R{t) 



(B.9) 



en el limite T ^ 00, esta resulta la unica contribucion no nula al coeficiente de difusion. En 
este caso, Uegamos finalmente a la formula de Kubo ()2.2ip . 



Dij = dt Rij{t) 
Jo 



(B.IO) 



Apendice C 

Soluciones exactas de la ecuacion de 
difusion 



Adoptando un modelo simplificado en el que los campos, los coeficientes de difusion y las 
fuentes tienen una dependencia espacial sencilla, la ecuacion de difusion puede resolverse 
anali'ticamente [59]. Las soluciones exactas que se obtienen en este modelo simple muestran 
el origen de la rodilla en el escenario de la difusion turbulenta y drift. 

Teniendo en cuenta una linica componente de campo magnetico regular, asumiremos 
que su amplitud es |Bo(r, 2)| ~ 1/r. De las ecs. fl2.26|) y ()2.28|) . resulta esencialmente 

1/3 

que D± (X r/ y Da oc r^, entonces, podemos despreciar la dependencia espacial de la 
difusion transversal {D±{r, z) = cte.) y asumir DA{r, z) oc r. Para las fuentes, consideramos 
Q{r,z)=2hsq{r)6{z). 

Si el campo magnetico regular no tiene inversiones espaciales (es decir, h^[r,z) = — 1 en 
toda la galaxia), las velocidades de drift (j2.14j) - (j2.15|) quedan reducidas a 

= , u, = -IDaIt = cte. , (C.l) 

y la ecuacion de difusion ()2.13p resulta 



d'^ d"^ Id Uz d 
+ 



2h 

N{r,z) = --^q{r)6{z) . (C.2) 



Como una simplificacion adicional, asumiremos que la altura del halo es muclio menor que 
el radio de la galaxia {H <C R), de modo que los gradientes radiales puedan despreciarse 
frente a los gradientes en z ^. Para z 0, la ecuacion se reduce a 

dz' ~ dz ' ^^-^^ 



^Otra manera de hacer nulos los gradientes radiales es restringirse al caso en el que las fuentes estan 
uniformemente distribuidas en todo el disco galactico (es decir, tomando q{r) — etc.). 
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z/H 



Figura C.l: Perfil vertical de la densidad de rayos cosmicos, en el caso de un campo magnetico 
regular simetrico con respecto al piano galactico, para diferentes valores dew = 2HDA/rD^. 

la solucion es 

^ ^ r (e"'(i-^/^)-l)(l-e-"') (z>0) ,^ 

donde w = 2HDA/rD±. De integrar la ec. ()C.2|) en el entorno del disco galactico (es decir, 
calculando lim^^Q He^^ •••) ^ reemplazar la solucion ()(y.4|l . resulta 

Hr) = ^[^\ q{r) • (C.5) 
w senh[w) 

Debido a la ausencia de drifts radiales, la dependencia radial de la densidad de rayos 
cosmicos es identica a la distribucion radial de fuentes; en particular, si g(ro) = en un dado 
ro, alii es tambien N{rQ,z) = 0. Para vq tal que g(ro) 7^ 0, el perfil vertical de la densidad 
de rayos cosmicos refleja la influencia de la corriente macroscopica Uz < 0, tal como muestra 

la fjgirm 

En el piano galactico, la solucion es 

N{r,z = 0) = -^q{r)tanh(^^^] . (C.6) 

Teniendo en cuenta que D± oc E^^^ y Da ex. E, y tomando q oc E~^ para el espectro 
diferencial de produccion en las fuentes, vemos que, a bajas energias (para HDa/tD^ <^ 1) 
es oc E~f^~^^^, mientras que a energias altas (para HDa/tDx^ ^ 1) resulta oc E^^^^. 
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Figura C.2: Perfil vertical de la densidad de rayos cosmicos, en el caso de un campo 
magnetico regular antisimetrico con respecto al piano galactico, para diferentes valores de 
w = 2HDA/rD±. 



En consecuencia, la solucion exliibe una "rodilla" (correspondiente a un cambio de mdice 
espectral igual a 2/3) que ocurre en torno de HDaItD^ ^ 1. 

En el caso de un campo regular antisimetrico respecto del piano galactico (es decir, tal 
que b^{r,z) = —sgn{z)), los drifts resultan 

2D A 

Ur = 2DaS{z) , = sgn{z) ; (C.7) 

r 

es decir, los drifts verticales estan ahora dirigidos hacia el piano galactico, y sobre el piano 
aparece un drift radial singular, dirigido hacia el exterior de la galaxia. La solucion (en la 
aproximacion H <^ R) es 

dy q{yr)y-'+^'i'-^~'") . (C.8) 

Debido al drift radial, la densidad en un dado ro recibe la contribucion de todas las fuentes 
ubicadas en r < ro, mientras que los drifts verticales tienden a concentrar los rayos cosmicos 
en la region del piano galactico; la fig l(I2l muestra un perfil vertical caracteristico. 
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